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A Notice to Readers 


Because the cost of publishing the proceedings of the Electromagnetic Wave Theory Sym- 
posium (vol. AP-4, no. 3; July, 1956) exceeded expectations, the funds of the PGAP were temporar- 
ily depleted. As a result, six papers originally scheduled to appear in the present issue had to be 
withdrawn in order to maintain the year-end treasury balance in the black. The withdrawn papers 
will be added to the January, 1957 issue. We apologize to the authors for this unavoidable and 


unexpected delay. 
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—The Editor 


news and views 


MESSAGE FROM THE SECRETARY 


I would like to report several items pertaining to our 
present activity. Each year there are, of course, several 
changes in the membership of the PGAP Administrative 
Committee. The following is a list of members with 
their responsibilities and the dates at which the various 
terms of office expire. 


Honorary Member 


Lester C. Van Atta, Hughes Research Laboratory, Cul- 
ver City, Calif. (Indefinite term: an honor accorded 
the PGAP founder) 

Members Ex-Officio—(Past chairmen continue as mem- 
bers ex-officio for three years after their term.) 

Philip S. Carter, Radio Corporation of America, Rocky 
Point, N.Y.—Development of PGAP (1957). 


_ Delmer C. Ports, Jansky and Bailey, Washington, D.C. 


—Development of PGAP (1959). 


Regular Members 


John I. Bohnert, Naval Research Laboratory, Wash- 
ington D.C.—Vice-Chairman (1958). 


Henry C. Booker, Cornell University, Ithaca, N. Y.— 


Chairman (1958). 

Arthur Dorne, Dorne and Margolin, Inc., 30 Sylvester 
Street, Westbury N. Y.—News and Views Editor and 
PGAP Representative to National Convention Program 
Committee. 

Frederick T. Haddock, University of Michigan, Ann 
Arbor, Mich.—WNot yet assigned (1959). 

Robert A. Helliwell, Stanford University, Stanford, 
Calif—Membership (1957). 

Jack W. Herbstreit, National Bureau of Standards, 
Boulder, Colo.—Not yet assigned (1959). 


Donald D. King, Radiation Laboratory, Johns Hopkins 
University, Baltimore, Md.—Promotion of Local 
Chapter Organization and Activity (1958). 

Robert L. Mattingly, Bell Telephone Laboratory, 
Whippany, N. J.—~Secretary-Treasurer and Public 
Relations Man Responsible for Publication of Promo- 
tional Information in Other Technical Journals (1959). 

Victor H. Rumsey, University of Illinois, Urbana, I].— 
Propagation Papers Procurement (1957). 

John B. Smyth, Smyth Research Associates, 3930 
Fourth Avenue, San Diego, Calif—TRANSACTIONS 
Editor (1958). 

Roy C. Spencer, Sylvania Electric Products, 100 First 
Street, Waltham, Mass.—TRANSACTIONS-PROCEED- 
INGS Liaison (1957). 

Archie W. Straiton, University of Texas, Austin, Tex.— 
Antenna Papers Procurement (1957). 

A recent ruling of the Administrative Committee will 
make available the roster of the PGAP membership 
(2211, of which 276 are students as of September 30) to 
any organization which is currently subscribing toa 
year’s institutional listing in our TRANSACTIONS. The 
charge for this information will be nominal. 

Our PGAP local chapters are continuing to grow in 
number. At present there are eight in existence and two 
more sections, Syracuse, N. Y. and Akron, Ohio, have 
expressed an interest in this affiliation. Items from six 
of these chapters will be found in the next section, 
“Chapter News.” The other two sections are Albuquer- 
que, N. M. and Orange Belt, Calif. The chairmen of 
these chapters from whom information regarding their 
activity may be obtained are respectively, G. Oltman, 
Sandia Corporation, Albuquerque, N. M. and T. J. Gib- 
ban, 204 Princeton Avenue, Claremont, Calif. 
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We have received many acknowledgments from 
foreign libraries attesting to the effectiveness of our 
policy of introducing our TRANSACTIONS into these 
libraries on a complementary basis. 

RoBERT L. MATTINGLY 


CHAPTER NEWS 
Chicago Section 


The first meeting for the 1956-1957 season of the 
Chicago Chapter of the PGAP was held on Friday, 
October 19, at the Western Society of Engineers Build- 
ing, 84 East Randolph Street. This was a field trip 
through WGN-TV’s new television station on top of 
the Prudential Building. They are operating at maxi- 
mum power of 316kw erp or as they say “the highest 
power from the highest tower.” 

The next meeting of the chapter will be November 
16 at the Western Society of Engineers Building at 
which time Dr. R. F. H. Yang of Andrew Corporation 
will give a talk on “The Development of Microwave 
Antennas.” 


Denver-Boulder Chapter 


On May 23, 1956, the Denver-Boulder Chapter of the 
PGAP celebrated the completion of its first operating 
year. 

The activities marking the anniversary began with a 
technical meeting, open to the public, held at the Radio 
Building of the NBS, Boulder Laboratories, and con- 
cluded with an annual meeting and evening banquet for 
members at the Memorial Building of the University of 
Colorado. A reception was held between the two events 
for members of the PGAP and their friends, at the 
home of Dr. and Mrs. E. K. Smith. 

The technical session featured Dr. A. D. Wheelon of 
the Ramo-Wooldridge Corporation of Los Angeles. Dr. 
Wheelon spoke on the subject of certain theoretical 
aspects of scatter propagation. The meeting concluded 
with the showing of slides taken during visits to England 
and other European countries by Drs. Wait and Cot- 
tony. 

The chapter has held meetings frequently since then 
at the National Bureau of Standards with Dr. James 
Wait of NBS, Chapter Chairman, presiding. The most 
recent of these were as follows: 

July 19, 1956—A paper entitled, “A Survey of Current 
Ionospheric Research at the Cavendish Laboratory,” 
was presented by Dr. Glenn Keitel, physicist, recently 
Fulbright Fellow at Cavendish Laboratory, Cambridge 
University, England. 

August 8, 1956—“Recent RDF Research at the Uni- 
versity of Illinois,” a paper by Dr. Harold D. Webb, 
Associate Professor and Co-Supervisor of RDF Re- 
search, University of Illinois, was presented. 

August 28, 1956—A paper entitled, “Turbulence in 
the Ionosphere” was presented by Dr. H. G. Booker, 
Professor of Electrical Engineering, at Cornell Uni- 
versity. 
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Los Angeles Chapter 


The Los Angeles Chapter of the PGAP held its 
September meeting jointly with the General IRE Sec- 
tion meeting for which it provided the speaker, Dr. 
L. C. Van Atta of Hughes Aircraft Company, who 
spoke on “Microwave Antenna Problems.” 

The newly elected officers of the group, Robert Elliott 
(Rantec Corp.), Chairman; Robert Spegen, Vice-Chatr- 
man; and Thomas Kinaga, Secretary-Treasurer, met 
with committee chairmen at an organizational meeting 
on October 11 to plan activities for the coming year. 


Philadelphia Chapter 


The first meeting of the Philadelphia Chapter will be 
held November 28, 1956. C. T. McCoy, of Philco, Re- 
search Division, will present a paper, entitled “Princi- 
pals of Low Noise Microwave Receiver Design” with 
emphasis upon the microwave properties of crystals and 
their significance to antenna problems. 


San Diego Chapter 


An election was held for the next year’s officers of the 
newly formed San Diego combined Professional Group 
of Antennas and Propagation and Microwave Theory 
and Techniques. Those elected were the following: Dr. 
John B. Smyth, Chairman; R. E. Honer, Vice-Chair- 
man; W. L. Teeter, Secretary. Dave Proctor was tem- 
porary chairman. 


Washington Chapter 


The Washington Chapter of the Professional Group 
on Antennas and Propagation announces the election of 
the following officers for the 1956-1957 year: William 
W. Balwanz, of Naval Research Laboratory, Chairman; 
Harry Fine, of Federal Communications Commission, 
Vice-Chairman; Herschel T. Ward, of Melpar, Inc., 
Secretary. Terms of office will run concurrently with 
those of the Washington Section Officers. 


RECENT MEETINGS 
URSI Y 


The Fall Meeting of the U.S.A. National Committee, 
URSI, was held October 11 and 12, 1956 at the Uni- 
versity of California, Berkeley, Calif. The participating 
URSI National Commissions were: Commission 2, 
Radio and Troposphere—Chairman, Dr. J. B. Smyth; 
Commission 3, Ionospheric Radio—Chairman, Dr. 
M. G. Morgan; Commission 4, Radio Noise of Ter- 
restrial Origin—Chairman, Prof. A. W. Sullivan. 

There were three sessions held each day with a dinner 
and social evening on October 11 and a U.S.A. National 
Committee meeting on the evening of October 12. 


Symposium on Optics and Microwaves 


In cooperation with George Washington University, 
the Optical Society of America, and the Office of Naval 


1956 


Research, the IRE Professional Group on Antennas and 
Propagation again presented an outstanding technical 
symposium in Lisner Auditorium, George Washington 
University, Washington, D. C., on November 14-16, 
1956. Devoted to the subject of “Optics and Micro- 
waves” the program comprised 24 papers and was ar- 
ranged to emphasize many significant aspects of the two 
fields so that it might be clearly seen where they differ 
and where they are alike. 


CCIR 


The Eighth Plenary Assembly of the International 
Radio Consultative Committee (CCIR) was held in 
Warsaw, Poland, between August 9 and September 13, 
1956. 

It was well-attended and much good work was ac- 
complished. There were in attendance some 400 dele- 
gates and supporting personnel from 41 countries which 
participate in this work. 

The members of the United States delegation were 
designated from the representatives of government 
agencies and private companies and organizations who 
had participated in the preparatory work in the United 
States for the Eighth Plenary Assembly and are as fol- 
lows: J. S. Cross, Chairman; E. W. Allen, Vice-Chair- 
man; D. K. Bailey; E. W. Bemis; J. B. Coleman; 
W. Q. Crichlow; H. E. Dinger; Mrs. F. T. Dowling; 
L. W. Garvey; W. D. George; Miss Gertrude Grabow- 
ska; Mrs. J. G. Harris; J. W. Herbstreit; A. G. Jensen; 
Miss E. K. Kraft; W. Mason; A. G. Skrivseth; G. S. 
Turner; A. P. Walker; N. White, and F. H. Willis. 

The scene of the conference was the Palace of Culture 
and Science, the 750-foot skyscraper which dominates 
the Warsaw skyline, in which excellent facilities were 
provided for the plenary sessions, for the meetings of 
the Study Groups which dealt with the specific tech- 
nical questions considered by the Conference, and for 
the personal comfort and convenience of those in at- 
tendance at the Conference. 

No attempt will be made here to discuss in detail the 
many technical questions which were considered, but 
some idea of their scope and significance can be ob- 
tained by referring to the titles of the fourteen Study 
Groups which have been organized to deal with them. 
In order they are: I—Transmitters; [I—Receivers; II] 
—Fixed Systems; IV—Ground-Wave Propagation; V— 
Tropospheric Propagation; VI—Ionospheric Propaga- 
tion; ViI—Standard Frequency and Time Signals; VIII 
—Monitoring; IX—Radio Relay Systems; X—Broad- 
casting; XI—Television; XII—Tropical Broadcasting; 
XIII—Mobile Systems, and XIV—Vocabulary. 

The work of these committees has been steadily in- 
creasing over the last few years to the extent that it was 
becoming difficult for the Study Groups to resolve the 
large number of problems confronting them during the 
short period of time available to them during the 
Plenary Assemblies. Consequently, it was agreed to re- 
organize the work so as to spread a large portion of the 
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load to the three-year intervals between assemblies. 
This, together with some rearrangement of work assign- 
ments to the various Study Groups, should result in 
more effective operation of the committee. 

The present Director of the CCIR, Dr. Balth van der 
Pol, is scheduled to retire on December 31, 1956. Six 
applications for the position were accepted and Dr. 
E. Metzler of Switzerland was elected as the new 
Director. 

An invitation was extended to the CCIR to hold its 
Ninth Plenary Assembly in the United States and was 
accepted by acclamation. Although the exact time and 
place was not fixed, it will probably be held in either 
1959 or 1960. 

On the whole the contributions of the United States 
were well-received by the Assembly and served as the 
basis for many of the final documents adopted. The 
continued participation of the United States in future 
CCIR work will have a profound influence on telecom- 
munication developments in other countries and will 
thus benefit the United States in the final analysis. The 
value of the CCIR as a forum for the exchange of views 
and information among the most qualified engineers and 
scientists of this and other countries should be recog- 
nized, and to this end, the CCIR should be supported 
by all those in the United States—both government and 
industry—concerned with the over-all improvement and 
advancement of communications. 

E. W. ALLEN, Ch. Eng. 
Federal Communication Commission 
Washington D. C, 


FuTURE MEETINGS 
Symposium on Very-Low-Frequency Propagation 


This symposium, sponsored by the National Bureau 
of Standards and the Professional Group on Antennas 
and Propagation of the IRE, at NBS, Boulder, Colo. 
is to be held January 23-25, 1957. The following com- 
mittees and their respective chairmen are responsible 
for its presentation: 


Steering Committee—J. R. Wait 

Local Arrangements—R. Silberstein 

Finance—J. R. Johler 

Publicity—C. H. Bragaw 

Panel Discussions—T. N. Gautier and R. A. Helliwell 

Regular Technical Papers—J. M. Watts and J. R. 
Wait. 


F. W. Brown, K. A. Norton, and R. J. Slutz will be act- 
ing in an advisory capacity. 

Persons wishing to attend this meeting should notify 
Mrs. M. Halter, NBS, Boulder, Colo., as soon as possi- 
ble. Although papers will still be accepted until Novem- 
ber 30, the program is essentially frozen now. The 
papers which have been accepted so far are listed on the 
next page. No more than one or two changes or addi- 
tions are anticipated. 
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“Some Physical Problems in the Generation and 
Propagation of VLF Radiation” by E. L. Hill, 
Department of Physics, University of Minnesota. 
“Studies of High Power VLF Antennas” by W. Gus- 
tafson and E. Devaney, U. S. Navy Electronics 
Laboratory, San Diego. 

“Some Properties and Applications of Magneto- 
Ionic Theory at VLF” by R. A. Helliwell, Radio 
Propagation Laboratory, Stanford University. 
“The Relation Between Group Delay of a Whistler 
and the Distribution of Ionization Along the Ray 
Path” by R. L. Smith, Radio Propagation Labora- 
tory, Stanford University. 

“Measurement and Interpretation of the Polariza- 
tion and Angle of Arrival of Whistlers’ by J. H. 
Crary, Radio Propagation Laboratory, Stanford 
University. 

“The Effect of the Earth’s Magnetic Field on the 
Transmission and Reflection of VLF Waves at the 
Lower Edge of the Ionosphere” by I. Yabroff, 
Radio Propagation Laboratory, Stanford Uni- 
versity. 

“Records of VLF Hiss at Boulder, Colorado, 
During 1956” by J. M. Watts, National Bureau of 
Standards, Boulder. 

“Extra-Terrestrial Origins of VLF Signals” by 
R. Gallet, National Bureau of Standards, Boulder. 
“Extensions to the Geometrical Optics of Sky Wave 
Propagation at VLF” by J. R. Wait and A. Mur- 
phy, National Bureau of Standards, Boulder. 
“Wave Guide Mode Calculations for VLF Iono- 
spheric Propagation Including the Influence of 
Ground Conductivity” by J. R. Wait and H. H. 
Howe, National Bureau of Standards, Boulder. 
“A Study of Signal-Versus-Distance Data at VLF” 
by J. L. Heritage and S. Weisbrod of Smyth Re- 
search Associates and J. E. Bickel of U. S. Navy 
Electronics Laboratory. 

“Basic Experimental Studies of the Magnetic Field 
of Electromagnetic Sources Immersed in a Semi- 
Infinite Conducting Medium” by M. B. Kraich- 
man, U.S. Navy Ordnance Laboratory, White Oak, 
Silver Spring, Md. 

“A Technique for the Rapid Analysis of Whistlers” 
by J. K. Grierson and L. R. O. Storey, Radio 
Physics Laboratory, Ottawa, Canada. 

“A Method to Interpret the Dispersion Curves of 
Whistlers” by L. R. O. Storey, Radio Physics 
Laboratory, Ottawa, Canada. 

“Relations Between the Character of Atmospherics 
and Their Place of Origin” by J. Chapman and 
E. T. Pierce, Cavendish Laboratory, Cambridge, 
England. 

“Survey of Investigations of VLF Propagation at 
Cambridge” by K. G. Budden, Cavendish Labora- 
tory, Cambridge, England. 
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October 


“A Study of VLF Ground Wave Propagation in 
Alaska” by G. M. Stanley, Geophysical Institute, 
College, Alaska. : 

“The Phase and Group Velocity of the VLF 
Ground Wave” by J. R. Johler, National Bureau of 
Standards, Boulder. 

“Polarization of the Ground Wave of a Radio 
Atmospheric” by A. W. Sullivan, University of 
Florida. 

“Noise Investigation at VLF by the National 
Bureau of Standards” by W. O. Crichlow, National 
Bureau of Standards, Boulder. 

“Spectrum Analysis of Spherics” by W. Taylor, 
National Bureau of Standards, Boulder. 
“Statistical Descriptions of Atmospheric Radio 
Noise” by A. D. Watt, National Bureau of Stand- 
ards, Boulder. 

“On the Polarization of Spherics” by A. G. Jean, 
National Bureau of Standards, Boulder. 

“The Effect of Receiver Band Width on the Amph- 
tude Distribution of VLF Atmospheric Noise” by 
F. F. Fulton, Jr., National Bureau of Standards, 
Boulder. 

“Our Present State of Knowledge of the Lower 
Ionosphere” by A. H. Waynick, Ionospheric Re- 
search Laboratory, State College, Pa. 

“Heavy Ion Effects in Audio-Frequency Propaga- 
tion” by C. O. Hines, Radio Physics Laboratory, 
Ottawa, Canada. 

“Some Recent Measurements of Atmospheric 
Noise in Canada” by C. A. McKerrow, Radio 
Physics Laboratory, Ottawa, Canada. 
“Performance and Design Criteria for High Power 
VLF Antennas” by W. M. Brown, Bureau of 
Ships, Washington, D. C. 

“Tonospheric Irregularities Causing Random Fad- 
ing of VLF ” by S. A. Bowhill, Pennsylvania State 
University, University Park, Pa. 

“75 Ke/Sec Pulsed Normal Incidence Ionospheric 
Sounding” by A. J. Ferraro, J. J. Gibbons, R. E. 
Houston, W. J. Ross, and E. R. Schmerling, 
Pennsylvania State University, University Park, 
Pa, 

“A VLF Antenna for Generating a Horizontally 
Polarized Radiation Field” by R. S. Macmillan, 
R. M. Golden, and W. V. T. Rusch, California 
Institute of Technology, Pasadena, Calif. 
“Artificial Generation of Low-Frequency Atmos- 
pheric Pulses” by M. M. Newman, Lightning 
and Transients Research Institute, Minneapolis, 
Minn. 

“The Numerical Solution of the Differential 
Equations Governing the Reflexion of Long Radio 
Waves from the Ionosphere” by K. G. Budden, 
Cavendish Laboratory, University of Cambridge, 
Cambridge, England. 
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Symposium on Communication Theory 
and Antenna Design 


The Air Force Cambridge Research Laboratories are 
co-sponsoring a Symposium on Communication Theory 
and Antenna Design with the Boston University, 
Physical Research Laboratories, to be held at Hillel 
House, Boston University, Boston, Mass., on January 
9-11, 1957. 

For further information write to: Miss Alice Cahill, 
CRRD, Hq. Air Force Cambridge Research Center, Air 
Research and Development Command, Laurence G. 
Hanscom Field, Bedford, Mass. 

Following is the tentative program. There will proba- 
bly be some additions, and there are two changes possi- 
bly to be made. 


January 9, 1957 
Morning 


Registration: 9:00-9:30 A.M. 

Opening Address: 9:30-9:45 a.m.—Dr. Duncan Mac- 
Donald, Boston University, Dr. L. M. Hollings- 
worth, AFCRC. 

Keynote Talk—Dr. Roy. C. Gunter, Jr., ARCRC- 
Clark University. 

Mathematical Introduction J—Dr. Charles Bumer, 
Clark University. 

Mathematical Introduction II]—Dr. 
Holt, AFCRC-Tufts University. 


F. Sheppard 


Afternoon 


Application to Electronics I—Dr. Arthur Kohlenberg, 
Melpar, Inc. 

Application to Electronics [I—Dr. Peter Elias, Mas- 
sachusetts Institute of Technology. 


Evening 


Symposium Banquet at Laurence G. Hanscom Field, 
Bedford, Mass.: 7:30 P.M. 

Speaker: Col. Gordon T. Gould, Jr., Director, Com- 
munications and Electronics Division, Hq. Air Re- 
search and Development Command. 


January 10, 1957 
Morning 


Application to Optics I—Dr. Edward O’ Neill, Boston 
University. 
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Application to Optics II—Lt. George B. Parrent, Jr., 
ATER. 


Afternoon 
Electro-Optical Systems in Cascade—Dr. Otto 
Schade, Radio Corp. of America. 
Application to Radio Astronomy J—Dr. R. N. 


Bracewell, Div. of Radio Physics, the Common- 
wealth Scientific and Industrial Research Organiza- 
tion, Australia. 

Application to Radio Astronomy JI—To be an- 
nounced. 


Evening 


7:00 p.m. Discussion Group on the Application of 
Antenna Problems: Group Leader: Dr. Roy C. 
Spencer, Sylvania Electric Cor Products Corp. 


January 11, 1957 
Morning 


Antennas I—Dr. John Ruze, AFCRC-Radiation En- 
gineering Laboratory. 

Antennas I]—Dr. W. H. Stell, Div. of Physics, the 
Commonwealth Scientific and Industrial Research 
Organization, Australia. 

Antennas III—To be announced. 


Afternoon 


Classified—Contributed Papers on Special Antenna 
Problems. 


OTHER NOTES 


There is still time to select topics for a PGAP sym- 
posium at the IRE National Convention in March. Dr. 
J. T. Bolljahn, of Stanford Research Institute, who con- 
tinues as PGAP Meetings Chairman, although he is no 
longer a member of the Administrative Committee, 
urges anyone with suggestions for suitable topics to com- 
municate them to him as soon as possible. 

The membership is reminded that material for publi- 
cation in this department is always welcome. News 
about PGAP members, special events, activities, or 
University courses, as well as comments or opinions on 
PGAP affairs are all solicited. Such items should be sent 
to Arthur Dorne, News and Views Editor, Dorne and 
Margolin, 30 Sylvester St., Westbury, N. Y. 


Lo Sea 
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contributions 


Circularly-Polarized Biconical Horns* 
C. GOATLEY? anp F. D. GREEN} 


Summary—A new technique is described for obtaining circular 
polarization from a biconical horn. This polarization is obtained 
through use of a spatial array of thin conducting elements between 
the cone faces; orientation of the elements varies with distance 
from the feed. The method has been used successfully at S band and 
X band. 

This technique was applied to horns which have maximum radia- 
tion intensity in the plane normal to the cone axes as well as to other 
horns with beams tilted out of this plane. The beam tilt was obtained 
by an appropriate choice of cone configurations. 


GENERAL DESCRIPTION 
Feed System 


HE BICONICAL horn is among the most versa- 

tile of omnidirectional antennas because of its 

simplicity of design and, with the appropriate 
type of feed, its excellent broadband characteristics. 
Early work on the antenna! showed that it could be de- 
signed for either horizontal or vertical polarization, the 
choice of which largely determines the type of feed 
required. 

Excitation of the antenna to produce vertical polari- 
zation (TEM mode) requires the least complex feed 
arrangement. This feed is easily formed by flaring the 
inner and outer conductors of a coaxial transmission line 


* Manuscript received by the PGAP, October 27, 1955; revised 
manuscript received, July 16, 1956. Presented at the National Elec- 
tronics Conference, Chicago, IIl., October, 1955. 

+ Melpar, Inc., Falls Church, Va. 

1W. L. Barrow, L. J. Chu, and J. J. Jansen, “Biconical electro- 
magnetic horns,” Proc. IRE, vol. 27, pp. 769-779; December, 1939. 


so that the inner conductor joins the upper cone and the 
outer conductor joins the lower cone. 

Horizontally polarized (TEo: mode) versions of the 
antenna can be obtained by utilizing feeds such as the 
loop, slotted waveguide, or tridipole. The remaining 
possibility to be discussed is that of producing circular 
or elliptical polarization. 

The most straightforward approach to producing cir- 
cular polarization would suggest the design of a circu- 
larly-polarized feed. The complexity of adaptable feed 
types made it desirable to find a simpler method. Ac- 
cordingly, use was made of the flared coaxial feed, which 
is mechanically simple and inherently broadband, and 
the antenna was altered between the feed and the aper- 
ture to obtain circular polarization. 


Depolarizing Elements 


The method of obtaining elliptical and circular 
polarization is somewhat unconventional. It consists 
basically of interposing a dense array of thin, conducting 
elements on cylindrical plastic forms between the faces 
of the cones. Fig. 1 shows an S-band model of the 
antenna disassembled to view the array configuration. 
This model is 12 inches in diameter and 6 inches high. 
The cylinders are formed from thin polyethylene sheet 
material; the conducting elements are of silver paint 
applied with a Speedball pen. 

Conducting elements on the innermost cylinder of the 
array are in a plane normal to the common axis of the 
cones, and those on succeeding cylinders are inclined in 
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Fig. 1—S-band horn with depolarizing elements. 
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EXPERIMENTAL DATA 
Polarization 


The alteration of the electric field caused by insertion 
of the elements is such that elliptical polarization of 
voltage axial ratio less than 5:1 is obtained over the 
frequency band from 2000-5200 me and the polarization 
is very nearly circular over a limited frequency band 
centering at 3500 mc. 

Polarization patterns of the antenna were measured 
by rotating a pyramidal waveguide horn in a plane nor- 
mal to the direction of maximum radiation. This ar- 
rangement permits only the measurement of the polari- 
zation pattern and did not provide information on the 
relative phase of the field components to determine the 
sense of rotation of the resultant field vector. 

Polarization patterns of the S-band horn measured at 
3000 and 5000 mc are shown in Fig. 2(a) and 2(b). Char- 
acteristic of the patterns is a clockwise rotation of the 


oS 
\\be 


180 
(b) 


Fig. 2—(a) Polarization pattern—3000 mc. (b) Polarization pattern—5000 me. 


equal increments at progressively greater angles so that 
the elements on the outer cylinder are at 45°. The ele- 
ments are one inch long (one-quarter wavelength at 
3000 mc) and are evenly spaced about each row. The 
rows are symmetrically disposed above and below a 
centerline about the circumference of each cylinder. 
Center-to-center spacing between adjacent elements in 
a row is 2 inches on the outer three cylinders and, to 
accommodate the desired element length, is extended to 
approximately twice that on the inner two cylinders. 
Radial spacing between each of the five cylinders is 
one inch. 


axis with increasing frequency and a smooth variation 
in the axial ratio as a function of frequency. A graph of 
the axial ratio variation over a limited frequency band 
is shown in Fig. 3 and indicates the close approach to 
circular polarization at 3500 mc as measured in the 
plane of maximum radiation. 


Radiation Patterns 


Radiation patterns measured in the vertical and 
azimuth planes at 3500 mc are shown in Fig. 4(a) and 
4(b) for both vertical and horizontal polarization. Pat- 
terns in the vertical plane showed slightly narrower 
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Asymmetry in the patterns is believed to be largely 
due to imperfections in the experimental model such as 
joints in the plastic cylinders, irregularities in the metal 
cones, and slight eccentricities in the coaxial feed at the 
flare joints. These irregularities also show up as slight 
variations in the polarization patterns as measured at 
different azimuth positions around the horn. 


Gain Measurement 


Gain measurements were made on the antenna for 
vertical and horizontal polarization at three azimuth 
positions spaced 120° apart about the horn. The results 
for vertical polarization are given in Table I as measured 
at one azimuth position. Differences in the measured 
gain occur at different azimuth positions, as would be 
expected from an examination of the radiation patterns. 
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Fig. 4—(a) Vertical plane radiation pattern. (b) Azimuth plane radiation pattern. 


beamwidths at all frequencies for vertically-polarized 
radiation than for horizontally-polarized radiation. 
Similarly, the sidelobe level is consistently lower for 
horizontal polarization over the entire frequency range. 
The narrower beamwidths, and to some extent the lower 
sidelobe level, can be attributed to the difference in 
aperture illumination for the two fields. The field pro- 
ducing vertical polarization illuminates the aperture of 
the horn uniformly while the field producing horizontal 
polarization is distributed sinusoidally across the aper- 
ture. 


The values given in Table I are representative of values 
obtained at all three azimuth positions. 


TABLE I 
GAIN OF ANTENNA FOR VERTICALLY POLARIZED SIGNAL 


Frequency (mc) Gain Over Isotropic (db) 
3000 1.46 
3500 0.70 
4000 —1.03 
4500 1.20 
5000 2.80 
5500 DeOs 
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Match characteristics of the S-band antenna are 
indicated by the vswr which remains less than 2:1 over 
a 2:1 frequency band. 


SCALING AND BEAM TILTING 


In order to achieve full utilization of the antenna, it 
was desirable to investigate the feasibility of tilting the 
beam out of the horizontal plane and scaling for use at 
other frequencies. Both endeavors were successful in 
scaling the antenna for use at X band, and the experi- 
mental model is shown in Fig. 5. When assembled, it is 
10 inches in diameter and 23 inches high. The lower cone 
is replaced with a disk to achieve beam tilting. 


Fig. 5—X-band horn with depolarizing elements. 


This model also contains 5 cylinders constructed in 
the same manner as those used in the S-band model. The 
cylinders are separated radially by 0.3 inches. There are 
5 rows of elements on the outer cylinder, 4 rows on each 
of the next 2, and 3 rows on each of the last 2 cylinders. 
Each element is 0.35 inches long, and center-to-center 

spacing between the elements is approximately 0.4 
inches on all cylinders except the inner one where the 
spacing is increased to 0.5 inches. Vertical plane radia- 
tion pattern of the scaled model is shown in Fig. 6 as 
measured at 10,000 mc. Use of a single cone with a 
circular plate in place of the lower cone was sufficient 
to tilt the beam approximately 15° out of the horizontal 
plane. Polarization patterns show the same general 
characteristics as those measured with the S-band 
model, but it was found that element densities apprecia- 
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Fig. 6—Vertical plane radiation patterns. 


bly lower than used in this model tend to introduce un- 
desirable variations in the polarization characteristics 
at different azimuth positions. 


ANALYSIS OF PHASE DIFFERENCE BETWEEN VERTICAL 
AND HoRIZONTAL FIELD COMPONENTS 


All of the present designs of the antenna have been 
developed experimentally. Information which would 
permit design of a horn with predictable characteristics: 
depends on knowing considerably more about the func- 
tion of the array of elements. It is evident that the ele- 
ments induce horizontal components of the field, but the 
complete array must also be considered as an artificial 
dielectric medium filling the space between the cone 
faces. The properties of certain types of artificial di- 
electrics can be determined analytically with a high de- 
gree of exactness. Such is not the case with an arti- 
ficial dielectric consisting of thin conducting pins.’ 
Empirical data on such a medium is available, but it 
does not appear to be directly applicable to the case 
where the pins vary in orientation through the medium. 
For this reason, an exact analysis of the mechanisms 
producing elliptical polarization from the antenna 
appears impractical. However, a qualitative analysis is 
useful to demonstrate why a phase difference is intro- 
duced between horizontal and vertical components of 
the field, and the concept of optical path length is con- 
venient for this purpose. Fig. 7 shows a cross section of 


2S. A. Schelkunoff and H. T. Friis, “Antenna Theory and Prac- 
tice,” John Wiley and Sons, Inc., New York, N. Y., pp. 577-579; 
1952. 
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Fig. 7—Cross section of biconical horn. 


the horn. R, is the radius at which the cross section 
height, a, becomes too small to propagate the horizontal 
field components. Let mq be the effective index of re- 
fraction’ for the horizontally-polarized field component. 
For this component it may be shown that my=Aj,/X- 
where ), is in effect the guide wavelength, and X, is the 
wavelength in a medium of index m,. The value of m, is 
assumed to vary as a function of radius because of the 
changing orientation of the elements forming the di- 
electric medium. From these considerations a general 
expression for 2; can be given in the form 


tn = fr(R). (1) 
It follows immediately that 
ee es (2) 
f,(R) 


where J is the free-space wavelength. 

The cutoff wavelength \, may also be obtained in 
terms of R from Fig. 7 and (1). Where a=2R tan a, 
dh. = [4R tan a] [f,(R)], and 


Ly 1 
ae SSS (3) 


peed val 
Xe 
the expression for the total index of refraction becomes 


1 
(deh . (4) 


3 a Gee tan : [f.(R) a 


The propagation velocity of the vertical field com- 
ponent depends only on the index of refraction of the 
dielectric medium. The index for that polarization can 
be immediately written in a general form as a function of 
R so that 


hy = al: (5) 


3 The index of refraction equals the square root of the diclectric 
constant. 
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Since the optical path length is defined as Jn dR, the 
phase difference between the components can be ob- 
tained in the form 


as te 
9 = =(f /(R)dR 
\ 0 


ik jo ieee 7 ‘\ (6) 
? y/ Ds ers tan ae 


Although the functions describing the effective index 
of refraction for either polarization are not explicitly 
defined, the form of (6) suggests that a phase difference 
should exist and that it will vary with frequency. 

Some question must exist regarding the choice of R, 
as the lower limit of the second integral since it involves 
the assumption that the horizontal field component is 
generated at R,. Such an assumption must be considered 
to have limited validity, but investigation has not been 
extensive enough to determine an apparent phase- 
center for the horizontal component. 

The relative amplitude of the field components is also 
dependent on the configuration of the array of elements, 
and it seems reasonable to expect that both orientation 
and density of the elements play a significant part. No 
investigation of this has been made experimentally or 
otherwise. 


APPLICATIONS OF THE ANTENNA 


One of the most interesting possible uses for the de- — 
polarized horn was discussed in a recent paper by 
Kelleher and Morrow.‘ Their particular application con- . 
sisted of employing the horn to feed a symmetrical 
toroidal lens which can be designed to shape the beam 
in the vertical plane. The probable appearance of such a 
system is shown in Fig. 8. At present its practicability 
seems limited only by the weight of the lens material. 


—— FEED 


LENS~ 


Fig. 8—Toroidal lens with horn feed. 
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Phase Centers of Microwave Antennas* 
DAVID CARTER} 


Summary—This paper is concerned with the location of the 
phase centers of microwave antennas. The inadequacy of conven- 
tional aperture theory for the accurate description of phase centers 
is discussed. Formulas are developed and, for numerical indications, 
calculations are made for paraboloidal reflectors of different f/D 
ratios and a class of primary patterns which provide an approximate 
representation of a great many common feeds. The results are pre- 
sented in graphical form to provide useful design information and 
show the dependence of principal E- and H-plane phase center 
location on feed and dish parameters. Contrary to the prediction of 
aperture theory, it is shown that the phase centers of axially sym- 
metric antennas are not in the aperture plane, but they are dispersed 
about it. 


HERE [Sa good deal of ambiguity on the location 
“De the phase centers of microwave antennas. Dif- 

ferent answers are obtained depending on the 
approach used. For some purposes, where antennas are 
rotated or used in interferometric systems, knowledge of 
phase center location is desirable and often necessary. 
For small antennas this information is rather easily 
measured in the laboratory. However, the large ranges 
required for the measurement of large aperture antennas 
at small wavelengths make such measurements ex- 
ceedingly difficult. Reflections in the range, mechanical 
inaccuracies, etc., may well make accurate measure- 
ments impossible. 

Theoretical calculations have been made to obtain 
phase distribution over the main lobe.! These calcula- 
tions are based on the aperture distribution method of 
calculating far-field patterns?~* and show the phase and 
amplitude distributions in the far field for different as- 
sumed aperture-field distributions. While this theory 
satisfactorily predicts the far-field amplitude distribu- 
tions for small angles off-axis, it cannot accurately locate 
the phase centers of real antennas. 

The major difficulty is the fact that the total aperture- 
field distribution is generally not known and not used 
in the formulations referred to above. To obtain these 
formulations from information that might usually be 
known in microwave antennas, such as feed patterns 
(amplitude, phase and polarization) and reflector geom- 
etry, approximations are made which seriously affect 
the location of the phase center. In the first place the 
location of the “aperture” plane is somewhat arbitrary.® 


* Manuscript received by the PGAP, November 10, 1955; re- 
vised manuscript received, July 30, 1956. 

+ San Jose State College, San Jose, Calif. 

1C. C. Allen, “Radiation patterns for aperture antennas with 
nonlinear phase distributions,” 1953 IRE CoNVENTION REcorD, Part 

nO=12" 

ft Te. C. Spencer, “Fourier Integral Methods of Pattern Analysis,” 
R. L. Rep. No. 762-1; January 21, 1956. 

8H. T. Friis and W. D. Lewis, “Radar antennas,” Bell Sys. Tech. 
J., vol. 26, pp. 232-246; April, 1947. 

4 Silver, S., “Microwave Antenna Theory and Design,” McGraw- 
Hill Book Co., Inc., New York, N. Y., 1949. 

5 Tbid., p. 158. 


If one employs the usual methods of geometrical optics, 
for example, to obtain the aperture-phase and amplitude 
distributions from a point source feed pattern and a 
paraboloid reflector,® one obtains the same distributions 
on any two parallel planes which are perpendicular to 
the axis of the reflector. Calculations of the far-field 
patterns from these aperture distributions by means of 
the Fraunhofer field approximation,’ 


\e~ TkR 


AR 


Up ~—— | P(E, nei trhy dtn, 

A 
then would be identical and one would obtain phase cen- 
ters differing in location by the axial separation of these 


two arbitrarily chosen planes. The notation used here is 
the same as that of S. Silver.* (See Fig. 1.) 


ia 


Fig. 1—Coordinate system. The origin is at the focal point of 
the reflector. 


A second significant discrepancy occurs in the fact 
that there is no counterpart of the longitudinal com- 
ponent of reflector current in the aperture distribution.® 
This longitudinal component of current makes no sig- 
nificant contribution to the amplitude patterns at 
angles near the axis, and it has generally been neglected. 
However, it does contribute to the curvature of the far- 
field phase pattern and is, therefore, significant in the 
calculation of phase centers. For this reason, phase 
centers derived from the far-field patterns of repre- 
sentative aperture distributions,® will not coincide with 
the actual phase centers of the antennas represented by 
these aperture distributions. In fact it will be shown 
that paraboloids illuminated by point source feeds with 
axial symmetric patterns have principal plane phase 
centers which are not coincident and whose location de- 
pends on the feed and dish parameters. This is in con- 
trast to the results of aperture theory which would pre- 


6 [bid., p. 419. 
Plbtd pels. 
8 [bid., p. 420. 
9 Allen, op. cit., p. 9. 
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dict coincident phase centers located in the aperture 
plane in such cases. 

This latter statement follows directly from the fact 
that the Fourier transform of an even function is real. 
Thus the Fraunhofer field of an aperture distribution, 
F(é, 7), is given by! 


i] : foo} i> 2) : 
PES oul if w(E, er batt kun) ded 
Pp XR ee ( q ) 


where u(&, 1) vanishes everywhere except inside the 
aperture area, wherein it coincides with the aperture 
distribution. Since the integral is real if w(&, 7) is a real 
(except for an arbitrary complex constant) function 
even with respect to both of its arguments, the phase of 
U, will be independent of direction, (6, @), except for 
discrete points when U, passes through zero making the 
phase change by z radians. In the case of the paraboloid 
illuminated by a point source feed with axial symmetric 
patterns, having negligible cross-polarization component 
of reflected field, the aperture distribution of the prin- 
cipal polarization component is axially symmetric and 
phase constant.'! Therefore the corresponding aperture 
distribution representation, “(£, 7) possesses more than 
enough symmetry to make its Fourier transform real. 
Hence the equiphase surface in the region of the main 
lobe in the Fraunhofer field is given by kR=constant, 
which is a sphere with center at the origin of coordi- 
nates. Thus, for this example, the aperture distribution 
approximation predicts that there will be a single point 
phase center located at the origin which is taken at the 
center of the aperture plane. 

To obtain a more accurate description of the phase 
center of such antennas, the phase distribution in the 
far field will be calculated from the current distribution 
over the reflector. Inclusion of the longitudinal com- 
ponent of reflector current will be seen to separate the 
principal H- and H-plane centers of phase. In addition 
this approach will locate these phase centers with re- 
spect to the vertex of the paraboloid. This will be done 
in terms of a primary feed gain parameter, and the 
paraboloids angular aperture to show the variation of 
these phase center locations with feed and dish pa- 
rameters. 

Choosing a class of primary patterns which provide 
an approximate representation of a great many common 
feeds, the axially symmetric point source primary feed 
gain function will be taken as (see Fig. 2) 


Ga(¥) = G(0) cos" v = G00) (- = =)". 
1+ x 


Then assuming as before that for most feeds of in- 
terest the cross polarization component of reflected 
field may be neglected in calculating the principal 
polarization diffraction patterns, it has been shown that 


10 Silver, op. cit., p. 174. 
1 Toid., p. 419. 
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x 
Fig. 2—Coordinate system. The aperture is in the x=y plane. 


principal Z- and H-plane patterns of the radiation field 
of such a paraboloid antenna is given by” 
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" D. Carter, “Wide angle radiation in pencil beam antennas,” 
J. Appl. Phys., vol. 26, pp. 645-652; June, 1955. 
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The notation here is the same as Carter’s,” and the 
origin of coordinates is taken at the focus of the pa- 
raboloid. From this expression for the principal E- and 
H-plane patterns, it can be seen that the equiphase 
contours are given by 


f 0 T.,sin6 + [2 cos @ 
const.— 27.cos¢— = tana! 
2 k T.2.sin@ — I, cosé@ 
R=, in the E-plane 
6 Ix 
const. — 27 eis alin = thea - plane 


te] 4 x1 


Here one can see the effect of the longitudinal com- 
ponent of reflector current in separating the principal 
&- and H-plane phase centers. Thus if J and J,. were 
negligible the equiphase contour in the £ plane would 
reduce to the same expression as obtained for the H 
plane. The fact that the longitudinal component of re- 
flector current contributes to the principal £-plane 
equiphase contour and does not contribute in the 
H plane is not surprising. Consideration of the distri- 
bution of this component on the paraboloid as repre- 
sented in Fig. 3 shows that it has odd symmetry about 
the H plane and the contributions from both sides of 
this plane cancel everywhere in the plane. This does not 
happen in the principal £ plane because the longitudinal 
component of reflector current has even symmetry 
about that plane. 


4-PLANE 


H-PLAME 


Fig. 3—Longitudinal component of reflector current. The distribution 
of this current component is represented by the horizontal arrows 
and the dots and crosses to illustrate the symmetry about the 
principal E and H planes. The vertical arrow at the focus indi- 
cates the polarization of the feed. 


To obtain the phase centers, the radius of curvature 
of the equiphase contours will be evaluated on the axis 
of the main lobe. If R(@) is the equiphase contour, the 
radius of curvature, p, is given by 
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The phase patterns of axially symmetric antennas are 
axially symmetric and have continuous first derivatives. 
Therefore, on the axis, the first derivative of RK with re- 
spect to @ vanishes and 


R(0) 


j € =) 
R dé? Jo 


ree =| 
p(0) = aR +| = ts 


Carter: Phase Centers of Microwave Antennas 


599 


in the far zone. Hence the phase center lies behind the 
origin a distance equal to 


PR] 
46? ty 


(0) — RO) =| 


It will be convenient to express the phase center loca- 
tions as fractions of the focal length in front of the vertex 
of the paraboloid. Calling this normalized distance ¢, it 
can be seen that 


i= ORO) 1-- [= 
dO? Ngo 


if di 


Performing the somewhat laborious operations indicated 
it can be shown that for 6=0, 
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This last formula locates the phase centers as a func- 
tion of the primary feed gain parameter, 2, and the 
angular aperture of the paraboloid, ¥=2 tan! x. 
Choosing a range of the parameters ” and WV to obtain 
representations of the most common configurations, the 
phase center locations were evaluated and plotted in 
Figs. 4 and 5. For purposes of comparison a plot of 
aperture plane location vs angular aperture was added 
to Fig. 5. 

Contrary to the prediction of aperture theory, Figs. 
4 and 5 indicate that the phase centers are not in the 
aperture plane, but they are dispersed about it. For the 
most part, they lie between the aperture plane and the 
vertex which is the region to the left of the aperture 
plane locus in Fig. 5. The H-plane phase centers always 
lie behind the aperture plane, near the vertex. This is 
also true for the E-plane phase centers in most cases. 
However, the axial component of reflector current, 
which does contribute in the & plane, pushes the phase 
centers further ahead of the vertex than their cor- 
responding location in the principal H plane. The reason 
for this is that the reflector’s curvature forces the axial 
component of current distribution to lie ahead of the 
vertex, producing a maximum of this component dis- 
tribution near the aperture plane. In some cases of low 
primary feed gain and small angular aperture, Fig. 5 
indicates that the principal E-plane phase centers are 
located in front of the aperture plane. 
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Fig. 4—Principal plane phase centers. These curves show the varia- 
tion of the phase center locations with primary feed gain in the 
principal E and H planes of a paraboloid reflector fed by a pri- 
mary having a gain function. 


2(n+1)cos*¥ for OSyS 
Gp = 


Finally it should be noted from Figs. 4 and 5 that, 
in general, decreasing the paraboloid angular aperture 
and increasing primary feed gain moves both principal 
E- and H-plane phase centers towards the vertex of the 
paraboloid. This was to be expected since either of these 
two variations tends toward the limiting case of reflec- 
tion by a plane mirror of a point source, a case in which 
image theory tells us that the phase center is located in 
the reflecting plane. Thus decreasing the angular aper- 
ture means either decreasing the paraboloid diameter 
for a fixed focal length or increasing the focal length for 
a fixed diameter. In the former case the aperture plane 
is moved back toward the vertex and in the latter case 
the vertex is moved forward toward the aperture plane, 
both variations effectively producing a flatter reflector. 
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Fig. 5—Principal plane phase centers. These curves show the varia- 
tion of the phase center locations with angular aperture in the 
principal # and H planes of a paraboloid reflector fed by a pri- 
mary having a gain function, 


2(n + 1) cos"y for 0s¥s> 
G = 


Similarly increasing primary feed gain increases the cur- 
rent in the zero curvature vertex region at the expense 
of the currents near the aperture plane, effectively pro- 
ducing a flatter current distribution. Since either of 
these variations produces a flatter current distribution 
located near the vertex, the phase centers should move 
together toward the vertex. 

In conclusion, it should be stressed that these results 
apply strictly only for feeds which are axially sym- 
metric point sources located at the focus of the parabo- 
loid. This provides an approximation for a great many 
common feeds. However, the phase centers of the re- 
flected radiation field can be varied and even pushed be- 
hind the dish by defocusing and separating the primary 
feed phase centers. 
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A New Method for the Measurement of the Average Dielectric 
Constant of the Underground Medium on Site* 
MOHAMED A. H. EL-SAIDt 


Summary—An electromagnetic interference pattern! in the far- 
distance condition is utilized to determine the average dielectric 
constant of the underground propagational medium. The method 
depends upon the determination of the surface wave velocity by 
means of measuring the first self-resonance frequency of a dipole 
wire laid on the earth’s surface. 


INTRODUCTION 


HE TECHNIQUE of sending radio waves to 

large distances into the interior of the earth has 

recently become important to the art of geo- 
physical prospection.! The dielectric constant of the 
underground material over the propagational path is one 
of the important parameters of this technique and a 
knowledge of its magnitude is essential both from the 
viewpoint of quantitative understanding of the propaga- 
tion of em waves into the earth’s interior, and also for 
prospecting purposes. Direct methods depending upon 
laboratory measurements on rock samples from the site 
are not indicative of the conditions prevailing in the 
underground propagational medium. This paper pro- 
vides a new method for the determination of the dielec- 
tric constant of a site with good accuracy. 


AVERAGE DIELECTRIC CONSTANT 


The inhomogeneity of the underground material 
causes the dielectric constant to vary from point to 
point along a propagational path. This will cause the 
phase velocity and the wavelength to vary as the wave 
progresses into the material. However, an average 
value of the wave velocity and wavelength may be en- 
visaged if we consider an average value for the dielectric 
constant. This is true as long as the wavelength is many 
times greater than the variational distances, which is 
the case in the operating frequency range utilized for 
prospecting. 

When a dipole is laid on the earth’s surface and is 
excited by a voltage, the wave which is set up will have 
in the far-zone three so-called fundamental components: 
one traveling and expanding into the upper space, 
another into the interior of the earth, and a third along 
the earth’s surface so that the boundary condition is 
satisfied. In the presence of an underground reflecting 
surface such as a water layer or a metallic ore, the 
underground wave will be reflected back to the earth’s 
surface thus interfering with the surface wave. It is 
known! that the resulting em interference pattern at 


* Manuscript received by the PGAP, December 28, 1955; revised 
manuscript received, April 3, 1956. 

t Dept. of Energ., Cairo University, Giza, Egypt. 

1M, A. H. EI- Said, “Geophysical prospection of underground 
water in the desert by means of electromagnetic interference fringes, ” 
Proc. IRE, vol. 44, pp. 24-30; January, 1956. 


distances on the earth’s surface which are much larger 
than the depth of the reflecting surface is characterized 
by fringes occurring at regular intervals given by: 


1a 
alt(ver- —) 


in the variable frequency and variable distance inter- 
ference methods respectively, where: 


(1a) 


I 


Delay frequency Fy 


(1b) 


Delay distance do 


c=velocity of light, 

€,=average dielectric constant of the underground 
material, 

a=ratio of the surface wave velocity to the velocity 
of light, 

f=operating frequency of the transmitter, 

d=distance between the transmitting and receiving 
dipoles. 


The values of Fy and do are independent of the depth 
h of the reflecting surface provided that d>>h. Also, the 
product Fod =fdo is substantially constant, a criterion 
which identifies the so-called “far-distance condition.” 
Eqs. (1a) and (1b) may also be re-written as: 


sree 
va= (+5 rioN ey ef pata 


This far-distance relationship comprises therefore a 
method for the determination of ¢, if a is known. The 
measurement of a on site can be effected by means of 
the following method. 


(1c) 


SURFACE WAVE VELOCITY 


A dipole laid on the earth’s surface will exhibit a self- 
resonance frequency which is less than that of the same 
dipole in free space by virtue of the fact that the di- 
electric constant of the earth’s surface material is 
greater than that of air. Fig. 1 shows the observed be-’ 
havior of the driving point impedance of a dipole an- 
tenna laid on the ground in the desert land during 
summertime; 7) and xo are the resistive and reactive 
series components, respectively. These components 
were measured by means of a balanced radio frequency 
bridge placed at the dipole center. It will be seen that 
these curves have the same general character as those 
of the same antenna in free space.? Further, the first 
self-resonance frequency of such laid dipole is approxi- 
mately half that of the first self-antiresonance. This 


2 See Ronald W. King, “Electromagnetic Engineering,” McGraw- 
Hill Book Co., Inc., New York, N. Y. 
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Fig. 1—Observed driving point impedance components of a di- 
pole wire laid on the earth’s surface in the desert. 


denotes that the voltage wave is reflected at the open 
end with negligible phase-shift. Moreover, the measured 
current distribution along this dipole wire when excited 
with a voltage at its first self-resonance frequency in- 
dicated a practically sine curve. 

In general, the driving point impedance of an excited 
dipole is a complex quantity, and when the wave set up 
by such antenna is reflected at the open end and reaches 
the driving point with a phase retardation of 7, the 
impedance becomes minimum and is resistive. Under 
this condition, the antenna exhibits its first self- 
resonance. For a dipole laid on the earth’s surface, the 
phase-shift depends upon the velocity of the wave 
which is set up by the antenna on the earth’s surface. 
The velocity of propagation may thus be determined by 
a measurement of the antenna self-resonance frequency. 
To a first approximation, this is the same velocity with 
which the surface wave propagates in the far zone, thus 
L=y,/4 where X, is the surface wavelength, from which: 


a = AfoL/¢ = I/We (2) 


where Z is the halflength of the dipole, fo is its first self- 
resonance frequency, and e, is a so-called “combined 
dielectric constant” which is between that of air and 
that of the earth’s surface material. 


EXPERIMENTAL WORK 


Fig. 2 shows the observed values of the length L 
against the first self-resonance frequency fo of the dipole; 
calculated values of the product fol and a are also 
shown. These curves are taken for a dipole wire having 
a diameter of 1 millimeter covered with rubber insula- 
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Fig. 2—Parameters of a dipole wire laid on the earth’s surface in the 
desert and excited from the center at its first self-resonance 
frequency; the deduced surface wave velocity is shown. 


tion whose diameter is 4 millimeters. It will be seen that 
the product foL, and hence also a, are substantially con- 
stant over a wide range of frequencies. An investigation 
of the parameters of the earth’s material over this fre- 
quency range was made by measurements on sample 
materials taken at various depths from the surface. 
These measurements were effected by filling a cylindri- 
cal air condenser with the material and measuring the 
susceptance and conductance components of the sample 
by means of the radio frequency bridge. Values of the 
dielectric constant ¢ and the resistivity p were calcu- 
lated accordingly, and the results are shown in Fig. 3. 
It will be seen that p decreases rapidly with depth 
owing to the increased moisture content which also re- 
sults in an increased value of e«. When the samples were 
left to dry in the sun, the remeasured constants assumed 
the values at zero depth. Further, the observed effect of 
the increase of frequency from 100 kc to 1000 ke on the 
surface sample was to reduce ¢€ from 8.2 to 6.3, that is, 
a ratio of 1.3. This ratio increases considerably with 
depth. The observed increase of a by a factor of 1.08 
over the frequency range from 200 kc to 1000 ke may be 
attributed mainly to the above mentioned changes of «€ 
with frequency. It should be considered, however, that 
the relative effects of the large changes in € with depth 
are very much reduced by virtue of the fact that the 
energy storage in the material decreases with depth be- 
cause of increased distances from the antenna wire and 
increased losses. Observed changes in resistivity will 
have only a second order effect on the actual value of a. 
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Fig. 3—Dielectric constant and resistivity of sample materials of a 
desert site against depth at frequencies of 100 kc and 1000 kc. 


Fig. 4 shows a sample of the observed far-distance 
interference patterns obtained at some place in Senai, 
Egypt, where the water table is at a depth of about 500 
meters. It is interesting to observe that the depth of the 
fringes is large. The distance over which this pattern 
was taken is d=2100 meters, that is, a ratio of 
m=d/h=4.2. Since it is not usually convenient in the 
field to use distances of the order of 10 times the depth, 
it can be shown that (see Appendix) for m=>2 we have: 


a liasian ia) 
(bat) [or 2), 


In the experiment of Fig. 4, observed Fy) =110 kc, meas- 
ured \/e,=I/a=2.68, and the calculated value of Ve, 
is 3.58 for that particular site. 

The results established here were primarily intended 
for use in the determination of the water table in the 
desert by means of the electromagnetic interference 
fringe method. For this purpose the proper procedure 
for prospecting would be as follows: 


(1d) 
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Fig. 4—Sample of the “far-distance” interference patterns of a desert 


site: Abu-Aweigla-Senai, point 2. d=2100 m, F=110 ke, /€ 
= 268 9/5) 5.58. 
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Fig. 5—Parameters for “far-distance” interference patterns, 
Egyptian desert, summertime. 


1) Determine a of the site. 

2) Obtain two interference patterns, one at an ap- 
propriate distance and the other at a far-distance. 

3) Estimate Ve, from (1c). 

4) Determine the approximate value of h from the 
information contained in the appropriate distance 
pattern and from this estimate m. 

5) From (1d) and the information contained in the 
far-distance pattern, obtain the correct value of 
er, 

6) Repeat 4) to determine the correct value of h if 
necessary. 


Typical values of the product /Pod=fdy and a@ are 
shown in Fig. 5 against \e,, for the Egyptian desert 
during summertime. 
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APPENDIX 
Let d=mh, where m>>1 at far distances. We then 


have 
Y 
(DD = a(1 + =), 
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and it can be shown that in the variable frequency 
interference method we get 


48.) 1 C 2 
€, = (— am 1 5 
vex(—+5)/( +5] 


In the variable distance interference method, at far 
distances we have: dnii—dn2d, and if m>>1 we can say 
that mMri2m, =m. With these approximations, it can 
be shown that 


» 
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The Image Method of Beam Shaping* 
PAUL T. HUTCHISONt 


Summary—The image scheme of beam shaping for use at micro- 
wave frequencies is described. Cosecant-squared radiation patterns 
are obtained with parabolic reflectors no larger than those required 
to give pencil beams of commensurate beam widths. Radiation pat- 
terns calculated using diffraction theory are compared with meas- 
ured patterns of a paraboloidal dish fed by a horn feed and one im- 
age horn. Experimental patterns are included to show the effects of 
variation of all parameters. A qualitative analysis of a paraboloid 
fed by a horn and several images is shown to agree with measured 
results. 


INTRODUCTION 
ee Sar antennas having a so-called “co- 


secant-squared” radiation pattern! in the verti- 
cal plane and a uniformly narrow beam in the 
other plane are in widespread use for airborne mapping 
and other radar applications. One scheme that has been 
used to obtain such a radiation pattern is the multiple 
feed type antenna.? This paper describes the image 
scheme of beam shaping which employs a section of a 
paraboloidal reflector fed by an electromagnetic horn 
and an image of the horn provided by a flat conducting 
ground plane below the axis of the reflector. Because of 
the similarity of these systems, a brief description of the 
former will be included. 
A paraboloidal reflector or dish has the property of 
transforming a spherical wave into a plane wave if the 
spherical wave source is placed at the focal point of the 


* Manuscript received by the PGAP, July 7, 1955; revised manu- 
script received January 31, 1956. 
Georgia Inst. Tech., Atlanta, Ga. Formerly with Raytheon 
Mfg. Co., Newton 58, Mass. 
1S. Silver, “Microwave Antenna Theory and Design,” McGraw- 
Hill Publishing Co., Inc., New York, N. Y., pp. 465-471; 1949. 
2 Tbid., p. 473. 


dish. If a small horn feed which closely approximates a 
spherical wave source is displaced from the axis of the 
dish, the resulting secondary beam will be displaced on 
the opposite side of the axis by an angle proportional to 
the feed displacement. This property of the paraboloid 
is the basic principle in the design of the image and the 
multiple feed antennas. 

If a paraboloidal reflector is fed by two horns, one of 
which is displaced from the axis, there will be two 
secondary beams with an angular displacement deter- 
mined by the angular separation of the horns. The 
amplitude of each beam is approximately proportional 
to the square root of the power fed to each horn. If the 
horns are spaced and fed properly, the combination of 
the beams will give a smooth cosecant-squared pattern 
over limited angles. 

A side view of an image type antenna employing a 
flat ground plane and a cut paraboloid of focal length F 
is shown in Fig. 1. The rays that are reflected from the 
ground plane seemingly come from a point d/2 below 
the ground plane. The feed below the ground plane is 
an image of the main horn, and its directivity is a func- 
tion of the main horn radiation pattern. The waves ap- 
parently emanating from the image are the waves that 
are reflected from the ground plane, so the power di- 
vision between the main and image horns is determined 
by the tilt angle ¢. If the wave is horizontally polarized, 
the rays undergo a 180° phase-reversal at the ground 
plane, so effectively the horn and its image are fed 180° 
out of phase. The image horn causes a secondary beam 
displaced from the axis of the paraboloid by an angle 
proportional to d/F and having an amplitude deter- 
mined by the tilt of the horn feed. 
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The problem is to find the radiation pattern of this 
image type antenna and to show how to control it. The 
discussion will be divided into the following categories: 
1) theoretical analysis, 2) the effects of variation of 
parameters and a comparison of calculated and meas- 
ured results, and 3) modifications of the ground plane 
to give cosecant squared coverage over wider angles. 


PARABOLA 


FOCAL 
a eg gan 
GROUND PLANE SES ay 


Fig, 1—Cross section of image type antenna with flat 
ground plane. 


THEORETICAL BEAM SHAPE 


The antenna to be analyzed isacut paraboloid of rev- 
olution with a flat ground plane displaced from the 
axis of the dish by a distance d/2. The small horn feed, 
horizontally polarized, is at the focal point and tilted 
an angle @ above the axis of the dish. The problem of 
finding the resultant radiation pattern of this antenna 
will be solved by finding first the secondary beams pro- 
duced by the main and image horns individually and 
combining the beams using the principle of superpo- 
sition. An intermediate step will be an analysis of the 
amplitude and phase variations of the field intensity 
across the mouth of the dish. 

The field distribution over the aperture of the dish 
due to the main feed, F(y), is determined by the radia- 
tion pattern and tilt of the horn. In Fig. 2(b) the ampli- 
tude and phase of the field intensity across the aperture 
due to illumination by the main feed are plotted as a 
function of y. Since the horn is at the focal point, the 
phase of the intensity is constant over the aperture 
y=-—B to y=A. The origin is arbitrarily placed at the 
point of maximum intensity to simplify the integrals 
encountered later. 


Hutchison: The Image Method of Beam Shaping 


605 


A 


> 


Fy, / 


Pe Sy 
ee / UNIFORM 
/ PHASE 


@ POSITION Y IN APERTURE 


MAGNITUDE AND PHASE 
OF INTENSITY ACROSS 
APERTURE 


Fig. 2—Horn and reflector relationship and illumination of 
reflector by main horn, 


The field distribution over the aperture due to illumi- 
nation of the dish by the image, G(y), is a function of 
the main horn radiation pattern as well as the displace- 
ment of the ground plane. G(y) will be maximum at the 
intersection of the ground plane and the parabola, and 
G(—B) is approximately equal to F(—B) in magnitude. 
The image is at a distance d below the focal point, and 
since an off-axis fed paraboloid is not a perfect beam 
forming element, the phase of G(y) varies with y. The 
phase is determined by the method of ray tracing as 
shown in Fig. 3. The total path of a ray from the image 


Fig. 3—Geometry for calculating pathlengths of 
rays from image. 


to the aperture is a+. All rays emanating from the 
main feed travel F+T before reaching the aperture of 
the dish, so 27/[(a+6) —(F+T) ] +7 will be the phase 
angle of G(y) relative to F(y). The expressions for a and 
b are: 
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1 
sry: [(4F2 + y.2)? + 16F2d(d + 2y:)]}2 (A) 


lee 7 2 


(4F? + y.2)? + 16F2dy, 


Fig. 4 shows the pathlength variation for one particular 
setting of the ground plane. For small values of d/F, this 
curve can be closely approximated by the straight 
dashed lines shown in Fig. 4. This approximation simpli- 
fies the problem of finding the shape of the secondary 
beam due to illumination of the aperture by the image. 
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Fig. 4—Pathlength variation as a function of position across 
mouth of reflector. 


In Fig. 5 the amplitude and phase of the field intensity 
due to illumination by the image are plotted as a func- 
tion of y. The phase difference of F(y) and G(y) at the 
origin is Bc +7, where ¢ is the difference in pathlengths in 
centimeters if 8 is in radians per centimeter. 

The resultant field intensity at any point is the 
combination of F(y) and G(y), taking into consideration 
their time phase difference. The total intensity at 
y=-—B is small since F(—B) and G(—B) are almost 
equal in magnitude and nearly 180° out of phase. As 
the tilt of the horn is changed, the magnitudes of F(y) 
and G(y) change, but the phases are unaltered. As d is 
increased, the slopes m; and my are changed, so the 
phase of G(y) is changed, but the other functions are 
constant. 

The next step in the solution of the problem is to 
calculate the secondary radiation pattern. By using 
Huygen’s Principle,? the secondary pattern in the yz or 
vertical plane can be expressed in terms of F(y) and 
G(y) over the rectangular aperture. 


5 Bell Laboratories Staff, “Radar Systems & Components,” D. 
Van Nostrand Co., Inc., New York, N. Y., p. 476; 1949. 
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Fig. 5—Tllumination of the reflector by the image horn. 
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y=A 


G(y) €i8u(sin 6— “ay (3) 


y=k 
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E is the electric field intensity at any point in the yz 
plane. 

K is a constant of proportionality. 

GB =27/d where Xd is the wavelength. 

F(y) is the illumination of the aperture by the main 
horn. 

G(y) is the illumination of the aperture by the image 
horn. 

m, and mz are slopes of the phase variation curve as 
shown in Fig. 4. 

c, s, A, B, and k are points or distances shown in 
Fig. 5: 


The first term in (3) represents the radiation pattern 
of the antenna when no ground plane is used. This 
pattern due to illumination of the reflector by the main 
feed will be called the main beam and will be maximum 
along the axis of the dish regardless of F(y). The re- 
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maining terms give the radiation pattern of the dish 
illuminated by the image. This secondary beam, called 
the image beam, will have maximum radiation at some 
angle above the axis of the paraboloid. Since the value 
of G(y) is comparatively small in the region of the 
aperture y=k to y=A, the image beam peak will be 
determined by the second term in (3) which gives 
maximum radiation at an angle @=arcsin m,. The factor 
BC +7 gives the phase difference of the real components 
of the field intensities in the main and image beams. 
Because of the positioning of the origin, the quadrature 
components of field intensity are small compared to the 
real components, so BC +7 gives a good indication of 
the total phase difference of the intensities. This analy- 
sis of a three-dimensional problem by a two-dimensional 
approach is obviously in error unless the variation of 
the field intensity across the aperture has the same 
variation in y for all parts of the aperture. This assump- 
tion was made and a comparison of the calculated and 
measured results shows that this is justifiable in most 
cases. Probably the main source of error in this analysis 
is the assumption that the horn feed is a point source. 


EXPERIMENTAL RESULTS 


The reflector used in this experiment was a rectangu- 
lar section of a spun aluminum paraboloid having a 
21-inch focal length and vertical and horizontal dimen- 
sions of 36 and 72 inches. The horizontal half-power 
beam width at the peak of the beam was 1.4° and the 
absolute gain 40 db operating at a wavelength of 3.2 cm. 
In calculating the radiation patterns of the antenna 
employing a flat ground plane, both F(y) and G(y) were 
assumed sinusoidal since the feed produced almost this 
type of illumination. The calculated patterns of the 
main and image beams and the phase difference of the 
field intensities in the two beams are shown in Fig. 6 
for d/F=0.072 and ¢ =29°. The phase difference of the 
field intensities at the origin due to the main and image 
horns is 36°. However, the phase difference of the radia- 
tion fields will not be exactly 36° since the relative 
phase of the intensity in each beam varies with @. This 
phase difference is never greater than about 30° in the 
region where a deep null might occur, that is the region 
of approximately equal intensities. If a cosecant- 
squared pattern is desired, the two beams must add so 
the power at any angle @ is varied as cosecant-squared 
theta, so any dip or irregularity in the pattern would 
ruin it. The theoretical and measured radiation pat- 
terns for this setting of the ground plane are shown in 
Fig. 7(a). As will be shown later, the phase relation- 
ships will not be so favorable for other values of d/F. 


VARIATION OF PARAMETERS 


To get aclear picture of the operation of this antenna, 
it will be helpful to observe the changes in beam shape 
as the following parameters are changed: 1) ground 


Hutchison: The Image Method of Beam Shaping 


607 


A= 3.2cm 
d/F=0.072 


FIELD INTENSITY IN UNITS 


PHASE DIFFERENCE OF FIELD INTENSITY 
INMAIN AND IMAGE BEAMS IN DEGREES 


54). Sons Gen oe be Awan 
ANGLE @ IN DEGREES 


g. 6—Relative amplitude and phase of field intensities in 
main and image secondary beams. 
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plane level, 2) tilt of the feed, and 3) frequency. 

As the ratio d/F is increased, the position of the 
image and hence the image beam is displaced further 
from the axis of the paraboloid. If the main and image 
beams are far enough apart, 7 or 8° in this case, there 
will be a distinct null about midway between them 
since neither radiation pattern shows any appreciable 
energy in this region. Experiments show the presence 
of a deep null (25 decibels) when d/F is greater than 
0.143. When d/F is between 0.072 and 0.12, the beams 
are positioned so the radiation pattern could be smooth 
similar to Fig. 7(a) if the phase relationships between 
the beams are correct. Calculations show a deep null is 
likely when d/F is 0.12 since BC +7 for this setting is 
180°. As discussed previously, the intensities are nearly 
in phase when ¢=29° and d/F is 0.072. For the inter- 
mediate setting of d/F equal to 0.095, the phase differ- 
ence is 108°, so a slight dip is observed between the two 
beams. Figs. 7(a) and 7(b) (next page) show the be- 
havior of the radiation patterns as d/F is changed. 

As the tilt @ of the horn feed is changed, effectively 
the power division between the main and image horns 
is changed. The relative magnitudes of the two second- 
ary beams will change, but their positions will not. As 
the horn is tilted, 8C 47 changes from 36° when ¢@ =29° 
to 80° when ¢ =40°, but the relative magnitudes of the 
real and quadrature components in each beam also 
change. Calculations show that changing the tilt of the 
feed has some small effect on the phase relationships, 
but the change in the relative magnitudes of the second- 
ary beams is the predominant effect. Fig. 8 shows the 
measured patterns for several values of ¢ when d/F is 
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Fig. 7—Calculated and measured radiation patterns of image 
type antenna for different values of d/F. 


0.072. Changing the magnitude of the image beam by a 
large per cent usually requires changing the radiation 
pattern of the horn feed. 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


October 


DECIBELS BELOW PEAK 
To) 


d=3.2cm,d/F=0072 


ei Gi aa O 2 4 6 
ANGLE @ IN DEGREES 


Fig. 8—Measured radiation patterns of image type antenna for 
several tilt angles of the horn feed. 


As the frequency is changed, say +5 per cent, the 
only change in the beam shape that should be observed 
is the usual change in beamwidths. If the ground plane 
is 2 inch, d/F=0.072, below the axis of the dish, 
GC+ varies from 25 to 47° at the two extreme fre- 
quencies, which may cause the beam shape to change a 
few per cent, but not enough to cause a noticeable dip 
or null. This antenna was tested over a 10-per-cent 
band, and the radiation patterns were almost identical 
except for side lobe levels. 


SELECTION OF FocAL LENGTH 


The phase difference of the field intensities in the 
main and image beams is determined by d, F, and the 
wavelength A, and the angular position of the image 
beam depends on the ratio d/F. 

The position of the image beam must be close enough 
to the main beam to make the radiation pattern a 
smooth curve similar to Fig. 7(a). From the calculated 
beamwidths of the main and image beams, the desired 
angular position of the image beam can be determined. 
The ratio d/F is then determined from the empirical 
expression: 


O(d/F) = tan 0. (4) 


O is the beam deviation factor which is defined as the 
ratio of the angular displacement of the secondary 
beam from the axis to the angular displacement of the 


4 Silver, op. cit., p. 488. 


1956 


feed from the axis. 9» is the desired position of the lobe 
or beam due to the offset horn or image. 

In most cases the antenna will have the horn tilted so 
the peak illumination from the main feed is about D/3 
above the bottom of the dish. D is the vertical height of 
the dish. If some value of F is assumed and the values 
of d and y,=D/3 are substituted into (1) and (2), the 
value of a+0 for this tilt angle can be calculated. In 
order to satisfy the phase relationships for a smooth 
pattern, @+6 should be about one-half wavelength 
greater than F+T. 


MULTIPLE IMAGES 


The addition of the main beam and the beam due to a 
single image gives very limited cosecant-squared cover- 
age especially when the main beam is narrow. By shap- 
ing the ground plane as shown in Fig. 9, two images are 


~~ *IMAGE | 
- IMAGE 2 


Fig. 9—Image antenna using ground plane to provide 
two images. 


provided. In the two-image system, each image effec- 
tively illuminates only a portion of the reflector, so the 
vertical beamwidths of the image beams are wider than 
beamwidths of the single-image system. By use of the 
dual image ground plane, the coverage was increased 
from about 5° as shown in Fig. 7(a) to about 16°. Chang- 
ing the shape of the ground plane to a curved surface in 
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Fig. 10—Radiation pattern of an image type antenna using a ground 
plane curved in one plane to provide multiple images. 
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the vertical plane feeds the paraboloid with a series of 
images which gives a smooth pattern as shown in Fig. 10. 

One advantage of this antenna is the high gain ob- 
tained because the presence of the ground plane allows 
more nearly constant illumination of the aperture in 
the vertical plane, without spillover. Also, reflections 
from the reflector back into the horn feed are a mini- 
mum since the horn is at a point of low intensity due to 
cancellation of the fields from the main and image 
horns. 


CONCLUSION 


The well-known image principle used in optics is 
applied to microwave antennas for the purpose of beam 
shaping. A method of calculating the radiation pattern 
of a paraboloidal dish fed by a horn and its image and 
the ways of controlling this pattern have been presented. 
Measured radiation patterns of the antenna fed by a 
horn and one image are shown to agree with the theo- 
retical patterns. A qualitative analysis of a ground 
plane which provides several images is substantiated by 
actual patterns which show this type antenna provides 
good cosecant-squared coverage over limited angles. 

The flat ground plane was first used by Paul Hines, 
of the Raytheon Manufacturing Company, to provide 
shaped beams from a horn feed of pillbox construction. 
The author is indebted to Dr. L. J. Chu, of the Massa- 
chusetts Institute of Technology, for his helpful discus- 
sions during the preparation of this paper. 
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Summary—Experimental measurements on three loop antenna 
configurations are presented. The technique for obtaining impedance 
and current distributions using a single-wire transmission line over 
an image plane is described with particular attention given to the 
difficulties encountered. The results are reproduced in graphical 
form, and for loops where theoretical results are available, curves com- 
paring theory and experiment are presented. 


INTRODUCTION 
4 | YHE CLOSED thin-wire transmitting loop in its 


various configurations has been useful for some 

time because of its valuable directional properties 

and broad-band operational characteristics. Several 
theoretical studies have been made on the circular 
loop,’ but due to the difficulties encountered in ob- 
taining solutions for the current distribution and input 
impedance, and the complicated numerical evaluation 
involved, only approximate theoretical results which 
agree qualitatively with experiment have been obtained. 
Recent solutions for the circular and square loops 
have been found respectively by Storer? and King.® 
Storer has used a variational method in which Hallén’s 
Fourier-series solution for the current has been modified 
to avoid the difficulties of convergence previously en- 
countered and hence facilitate the numerical computa- 
tions. As a result, a series of tables and curves has been 
compiled of impedance and current distribution on 
circular loops of circumference up to 2.5 wavelengths for 
several thicknesses of wire. Theoretical data obtained 
by Storer provide a quantitative comparison in this 
paper. King has applied an iteration method to the 
problem of the square loop driven from one corner or 
driven from the center of one side. By superimposing 
the solutions of four independent integral equations 
representing four-phase sequences when the loop is 
driven at each corner by arbitrary voltages, a general 


* Manuscript received by the PGAP, October 5, 1955; revised 
manuscript received, January 31, 1956. Paper given at the 1955 
National Conference on Aeronautical Electronics, Dayton, Ohio. 
Published as Cruft Lab. Tech. Rep. No. 213, Harvard Univ., Cam- 
bridge, Mass., May, 1955. This research was ‘made possible through 
the joint support of the Navy Department (Office of Naval Re- 
search), the Signal Corps of the U. S. eine and the U. S. Air Force, 
under ONR Contract NSori- HAS), ABA OY 

+ Gordon McKay Lab., Harvard Univ. Cambridge, Mass. 

1K. Hallén, “Theoretical inv vestigations into the transmitting and 
receiving qualities of antennae,’ Nova Acta Upsala, ser. IV, vol. 11, 


pp. 1- 44: 1938. 

2G. Glinski, “Note on circular loop antennas with non-uniform 
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1947. 
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IRE Trans., vol. AP-4, pp. 393-407; July, 1956. 


Proc. Camb. 


ON ANTENNAS 


AND PROPAGATION October 


Measurements* 


KENNEDY} 


solution for the current distribution is obtained. A 
similar solution will result for the square loop driven at 
the center of a side. Numerical evaluation of. the results 
is in process. 


IMAGE-PLANE MEASURING TECHNIQUE 


By the theorem of images a single conductor over an 
infinite, perfectly conducting ground plane produces 
fields identical to those of a single conductor and its 
image. If the conditions of a plane surface, perfect con- 
ductivity, and infinite extent are satisfactorily approxi- 
mated for the ground plane, very good experimental re- 
sults can be obtained by the image-plane technique. In 
addition to simplicity of construction, other advantages 
of the system are 1) that it is inherently balanced and 
thus a high degree of electrical symmetry obtains, and 
2) the presence of the image plane serves as an electrical 
shield to isolate operator and measuring equipment from 
the antenna. 

The square and circular loops of interest are used to 
terminate a single wire transmission line mounted over 
an image plane. The input impedance of the loop is 
measured by sampling the rf signal along the driven 
transmission line with a small movable probe. The 
standing-wave-ratio or width-of-the-minimum method 
is used to determine impedance. The amplitude and 
phase of the current are obtained with a similar probe 
moved along the loop itself. The phase is determined by 
mixing the signal picked up on the loop with that from 
a matched transmission line and varying the position of 
the probe on the matched line to give a minimum signal. 


Transmission-Line Description 


The image plane used in this investigation (see Fig. 1) 
is of $-inch duralium and measures 12 feet by 10 feet. 
The center section, comprised of a two-foot strip, is 
highly polished and divided into four removable panels 
for easy adaptation to other measuring techniques. In 
the near vicinity of the transmission line and antenna it 
was found necessary to provide a smooth, highly con- 
ducting transition from one panel to another with the 
aid of silver paint and aluminum foil. 

A single-wire transmission line of }-inch brass tubing 
is mounted over part of the polished section of the image 
plane. A coaxial driving unit incorporating an ordinary 
type-V connector is soldered directly to the line. An 
inner line conductor terminated in a small search probe 
serves as a means of moving the probe along a 1/16-inch 
slot in the outer conductor as well as transmitting the 
rf signal to the detecting system. A rack-and-pinion de- 
vice of fine pitch provides an accurate method of meas- 
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Fig. 1—Vertical image plane. 


uring the motion of the probe. Adjustments on the 
carriage allow for a wide range of line spacings and for 
rotation of the slot through 180°. For the measurements 
described in this report a full line spacing of 1.246 cm 
was maintained, establishing a transmission-line charac- 
teristic resistance of 


b 
Uv O Cosa (—) = 155.4 ohms, 


2a 


where 0 is the distance between centers of the line and 
its image, and a is the radius of the conductor. 

The transmission line is supported at a position \/4 
behind the driving point by a short-circuiting plate. 
This quarter-wave section presents a high impedance to 
the generator. However, the distance from driving 
point to the load is nearly 1} meters long and additional 
support is needed. It was established experimentally 
that the smallest usable bead of polyfoam (that which 
would adequately support and space the line) had the 
effect of an added shunt capacitance across the line. 
The effect could be minimized by placing the beads at 
positions of minimum impedance along the line where 
the shunt effect would be small. However, for measure- 
ments involving continual end-load variation, it was 
not feasible to reposition the supporting beads for each 
new measurement. Therefore, the use of a single poly- 
foam strip supporting and spacing the entire length of 
line was decided upon (see Fig. 2). With such a strip, 
uniform and consistent measurements were made along 
the line and the results were interpreted in terms of the 
measured transmission-line wavelength, \”. The wave- 
length in a medium of polyfoam, \,, was calculated as 
follows: 
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where €,/e=1.05=the relative dielectric constant of 
polyfoam 


40 
Mp = —— = 39.02 cm 
1.025 
\=free space wavelength. 

The measured \” must be between the values for free 
space and polyfoam; that is 39.02 cm<\”<40 cm as 
only a small portion of the space between line and image 
plane is filled with polyfoam. The measured value of \” 
was 39.62. The values of phase constant, 8, and char- 
acteristic resistance, K,, were adjusted to the following 
new values for use in computing impedance. 


120 


b 
cosh} (~) = 
WV € ett 2a 


eft dr” peg 
= (=) = relative effective dielectric constant. 
€ 


154.1 


l| 


2a 
(ie 0 ==) (Oboe IK 


These new values apply to impedance measurements for 
the circular loop. The square loops were not supported 
in this way and the original values of 6=0.157 and 
R,=155.4 ohms apply to them. The image-plane system 
described above is a modification of that discussed by 
Angelakos® and Conley.’ 


IMAGE PLANE 


TRANSMISSION LINE 


POLYFOAM 
SUPPORT 


Fig. 2—Cross section of supported transmission line. 


Measuring Equipment 


The arrangement of the measuring equipment is 
shown in the block diagram of Fig. 3. A standard fre- 
quency generator (Measurements Corporation, Boon- 
ton, N. J., Model 84) having a maximum output power 
of 200 milliwatts was used to produce a frequency of 
750-mc. A double-stub tuner and a single-tube tuner are 
attached to the output. At the correct 750 mc setting 
the generator was found to produce a substantial 
520-mc signal in addition to the desired 750 mc. The 
single stub is used to detune the 520-mce signal and the 
double stub as a matching device for the 750-mce signal. 

6D. J. Angelakos, “Current and Charge Distributions on Anten- 
nas and Open-Wire Lines,” Cruft Lab. Tech. Rep. No. 98, Harvard 
University, Cambridge, Mass., 1950. 

7 P. Conley, “Impedance Measurements With Open-Wire Lines,” 


Cruft Lab. Tech. Rep. No. 35, Harvard Univ., Cambridge, Mass., 
March, 1948. 
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Fig. 3. 


A waveguide-beyond-cutoff attenuator is used to vary 
the power fed to the line. The rf signal is picked up from 
the transmission line by a charge or current probe and 
mixed with a 1000-cps modulation signal through a 
crystal. A U. S. Navy Radio Receiver Model RDO 
tuned to 750 mc is used as an amplifier and detector. 
The detecting system was calibrated and all meter 
readings corrected against the calibration curve. The 
1000-cps output from the receiver was filtered through a 
noise filter and observed on a Ballantine voltmeter and 
an oscilloscope. 


Measuring Procedure 


Impedance measurements® were made using the 
standing-wave ratio method for antennas which were 
well matched; the width-of-the-minimum method was 
used when twice-the-minimum power point could be 
located accurately. The hyperbolic solution of the 
transmission-line equations shows that 


él 
sinh (ad + p) = sing (5) 


and 
coth p = vswr 


where 


y =a+j6 = propagation constant 

6/=width of the minimum at twice the minimum 
power point 

d=distance from line termination to probe 

p=terminal attenuation function. 


For high standing-wave ratios it is apparent that ad 
and p are of the same order of magnitude, and that 
attenuation must be considered. A half-line attenuation 
constant of 2.5710-5 nepers/cm was calculated and 
used in many of the final computations. All measure- 


8D. D. King, “Impedance measurements on transmission lines,” 
Proc. IRE, vol. 35, pp. 509-513; May, 1947, 
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ments were made on a half-loop configuration which, 
when referred to the characteristic resistance of the full 
transmission line, resulted in the input impedance of the 
complete loop. 

Current distribution measurements were not made on 
square loops, since it was not feasible to maneuver a 
probe around a series of right-angle bends. Distributions 
on circular loops involved slotting these loops parallel 
to the plane of the loop and moving a shielded-loop 
probe®? around the loop contour. Measurements near 
the junction were avoided since errors due to coupling 
to the transmission line currents would be involved. 

The slotted line shown in the block diagram of Fig. 3 
was used for phase measurements. If this line is well 
matched, the error due to vswr >1 is kept at a minimum. 
A signal is picked up by the search probe at an arbitrary 
point on the loop and mixed with a signal similarly 
coupled out of the matched line. If the probe position 
on the matched line is varied to give a minimum signal 
on the Ballantine voltmeter and oscilloscope, the two 
signals are 180° out of phase. For a new position of the 
antenna probe another position of the matched-line 
probe will give a similar minimum signal. The difference 
in position of probe along the matched line from some 
arbitrary reference position may be expressed in degrees 
and is a measure of the difference in phase between sig- 
nals at the two positions on the image line. The position 
of the minimum can only be obtained accurately if the 
minimum is sharp. Therefore the amplitude of the signal 
from the matched line must be very nearly that of the 
signal from the search probe on the antenna itself. This 
was most easily done by varying the probe depth in the 
phase line and calibrating this change directly in de- 
grees of phase shift. A small correction was then applied 
to the phase angle. 


Antenna Construction 


Loop antennas of {-inch brass tubing (radius =0.318 
cm) were designed to cover a range of circumference or 
perimeter from \/2 to 3. The conventional parameter 
for the thickness of an antenna is 


Pie 
Q = 2\In—, 
a 


where 


2h =circumference or perimeter of loop 
a=radius of conductor. 
With constant radius the longer antennas are electri- 
cally thinner. One set of data was taken on square 
loops of 0.106 inch conductor diameter driven from one 
corner. The physical thickness of these loops was such 
that the conductor could be readily shortened and re- 
bent from the 2) to the \/2 loop size, and from the re- 


§R, King, “Theory of Linear Antennas,” Harvard University 
Press, Cambridge, Mass., pp. 129-131, 1956. 
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b= 1246 cm 
a= d/2=.318 cm 
perimeter 


he 2 
p'= actual conductor length 


2 
= h-.88Il cm 
b= 1.246 cm 
a= d/2=.134 cm 
perimeter 


2 
+ actual conductor length 


2 
= h -—.881 cm 


b= 1.246 cm 
a= d/2= 318 cm 
_ perimeter 


2 
= actual conductor length 
2 


Sey De 
C. SQUARE LOOP DRIVEN FROM CENTER OF ONE SIDE ~~ 2 °™ 


1.246 cm 
d/2 = 318 cm 
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2 
;_ actual Conductor length 
Bee a Dae fa 
b 
=h- 2 cm 


D. CIRCULAR LOOP 


Fig. 4—Loop configuration and notation (image-plane method). 


sults a fairly accurate prediction could be made about 
the input impedances of thin loops. 

The critical dimensions and mounting procedures are 
given in Fig. 4. Although all loops were constructed on 
the basis of loop size representing the complete square 
or circular configuration, it was decided to let the an- 
tenna size be determined by the actual length of con- 
ductor. In the loop case this implies an effective 
perimeter or circumference slightly smaller than the 
actual one. Results obtained will be referred to these 
effective dimensions. 

Both types of square loops were constructed using a 
set of angle bends to which various lengths of tubing 
could be fitted. Contact to the image plane was made 
by tapping holes in the plane and screwing a slotted 
brass plate directly to it. A few sets of holes permitted 
the mounting of the complete set of antennas and 
neither the holes in the ground plane nor the screw 
heads extending slightly above the surface of the plate 
were found to have any effect on the measurements. To 
avoid further defacing of the ground plane the circular 
loops were mounted at a single contact point by a 
screw from the back of the screen, and connection to 
the transmission line was made by a series of appropri- 
ate lengths of tubing supported by strips of polyfoam. 
Had the extent of the undertaking been realized at the 
beginning of experimentation, more consistent mount- 
ing methods could have been used. 
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PRESENTATION OF DATA 
Input Impedance 


Impedance data have been obtained for several loop 
sizes between A/2 and 3d. These experimental results are 
available in tabulated form but appear here in Figs. 5 
through 11. The procedure for procuring the final 
impedance values is partially mathematical and par- 
tially graphical. Standard Smith chart representations 
of the normalized impedances were obtained from 
measurements of standing-wave ratio and position of 
voltage minimum. Fig. 5 (next page) compares im- 
pedance spirals for loops of two different conductor sizes 
which, for purposes of comparison, may be described 
as thick and thin. Another Smith chart, Fig. 6, compares 
loop measurements for three loop configurations, keep- 
ing the conductor size constant. 

The actual plotting of impedances in Smith chart 
representation is facilitated by the spiral form which 
individual curves take. The need for a vast number of 
experimental points in the vicinity of peaks and dips is 
eliminated. The normalized values of resistance and 
reactance can be read directly from the Smith chart, and 
when multiplied by the characteristic resistance of the 
transmission line along which the measurements were 
made, produce the impedance of each antenna meas- 
ured. In addition, values are easily extrapolated to 
yield additional Boh’ =27h'/X or Bo" h’ =2rh’/d” values. 
The linear or semilog plots of the final impedances were 
obtained by this partially graphical method, the irregu- 
larities of the curves requiring some extrapolated values. 
Reading directly from a Smith chart in regions of high 
standing-wave ratio can produce errors caused by the 
limited accuracy of that part of the circle diagram. For 
such values it is better to plot the normalized admit- 
tance and compute the impedance from it. 

As a supplement to the Smith charts curves have 
been plotted for the thin square loop driven from one 
corner in Fig. 7 (p. 616); the thicker square loop driven 
from one corner in Fig. 8; the square loop driven from 
the center of one side in Fig. 9; and the circular loop in 
Figs. 10 and 11. Figs. 10 and 11 also contain the theo- 
retical impedances and admittances for antennas of 
comparable dimensions obtained by Storer.* 


Current Distribution 


Current distributions on circular loops have been 
plotted as a function of probe position around the loop. 
The distance around the half-circumference of each 
loop is calibrated in degrees (the angle which the probe 
position forms with the image plane) with the driving 
point being considered as ®=0° and the short circuit 
@=180°. The distance around the half-loop is also 
calibrated in fractions of a wavelength. 

The search probe used to determine the amplitude 
and phase of the current was mounted in a slot in the 
plane of the loop. Therefore measurements very near 
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Fig. 5—Impedance of square loop driven from one corner. 
\ 


= 180° were not considered accurate due to the small 
circular short-circuiting plate and the coupling of the 
probe to its own image. Measurements near the driving 
point were even less reliable since the orientation of the 
probe was such that a maximum of coupling to the 
transmission line occurred (see Fig. 12, p. 617). The ap- 
parent distribution then was a result of a combination 
of currents on both antenna and transmission line. Quite 
arbitrarily, results obtained by the probe within 2 cm of 
the junction were ignored, and where feasible the curves 
were extrapolated to ®=0°. 

The most reasonable way of presenting the amplitude 


of the current would be to assume a driving voltage of 
1 volt and normalize | Z| in amperes/volt to the input 
admittance of each antenna. Since the reliability of the 
input current data was known to be quite poor, the 
current amplitude curves were presented as relative to 
the current maximum obtained at 6=180°. In Fig. 13 
(p. 618) this maximum was arbitrarily given the value of 
one. In Fig. 14 where comparison with theory is possible, 
a more reliable current maximum was available, and 
two sets of experimental data were normalized to the 
theory at ®=90°. 

The position ®=90° is equally distant from both the 
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Fig. 6—Impedance of three-loop configurations. 


driving point and the short circuit and is likely to pro- statement about the effect of conductor size on the 
vide the most reliable data. This point was taken as a impedance of loop antennas. It is apparent that the 
reference for the phase measurements and, for purposes’ thinner antenna has a higher input impedance at anti- 
of presentation, the phase angle, 0, is set equal to 0° at’ resonance than the thicker antenna, and its resonance 
&=90°. occurs for a smaller loop perimeter. A superposition of 
the semilog plots of Figs. 7 and 8 would illustrate this 

CONCLUSION more clearly. 
The input impedances of the square loops of different It is interesting to compare the input impedance of a 
thickness shown in Fig. 5 provide a basis for a general loop antenna with a dipole of corresponding length and 
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conductor size. Kennedy and King!® include circle dia- 
grams and linear plots which should be compared with 
the experimental results from this paper and the circu- 
lar-loop theory of Storer.4 The resonant and antireso- 


10 Pp, A. Kennedy, and R. King, “Experimental and Theoretica 
Impedances and Admittances of Center-Driven Antennas,” Cruft 
Lab. Tech. Rep. No. 155, Harvard Univ., Cambridge, Mass., April, 
1953. 


Fig. 11—Admittance of circular loops. 


1956 


ANTENNA 


IMAGE PLANE 


Fig. 12—Probe orientation. 


nant antenna sizes compare in the following way: the 
first antiresonant loop has a resistance and reactance of 
very large magnitude, and a comparison of the input 
impedance spiral of the loop with that of the dipole 
shows the sequence of resonance and antiresonance to 
be interchanged, but to be quite comparable in general 
outline. Note, however, that the loop antennas are more 
capacitive than dipole antennas of corresponding size. 

A comparison of loop antennas of all three configura- 
tions is found in the Smith chart of Fig. 6. It is apparent 
that the curves lie quite close together when 2h=)/2, 
d, and 3/2. When 2h = 2), there is a wide spread in input 
impedance for the three configurations—the square 
loop driven from the center of one side having the 
largest impedance, the square loop driven from one 
corner having the smallest impedance, and the circular 
loop falling in between. The curves approach each other 
again at 5/2 and are close at 3X. 

A plausible explanation is found for this curve separa- 
tion if a cosinusoidal distribution of the current about 
each of these contours is considered as in Fig. 15. For 
the loop size 2h=2), configuration (a) has a current 
minimum at all four corners. Therefore, a maximum 
concentration of charge exists at these points. Con- 
figuration (b) represents the opposite condition and 
(c) can be said to have a maximum concentration of 
charge in several bent regions. The condition in (a) and 
to a lesser extent in (c) can be interpreted as an effective 
lengthening of the antenna or compression in the dis- 
tribution pattern. 

The amount of compression is then a measure of the 
capacitive loading at the corners or bends. The input 
current will be affected by the compression in distribu- 
tion pattern, and the input impedance will vary inverse- 
ly as the current. 

A plot of current distribution (Fig. 16) obtained ex- 
perimentally by Kaliszewski is of value. The com- 
pression in distribution pattern is discernible from the 
shifted minima—the “center-driven” square loop being 
shifted the most and the “corner-driven” square loop 


1 T, Kaliszewski, “Complementary Loop Antennas,” Cruft Lab. 
Tech. Rep. No. 241; February, 1956. 
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being shifted the least. The currents at the driving point 
then differ by an amount proportional to the shift. The 
shift of the circular loop is exaggerated since the only 
available data are for a loop of circumference 2.2\. The 
relative values of the input current from Fig. 16 mean 
little since there is some error in the peak voltages. How- 
ever, a qualitative estimate of | Zo| can be made from 
the shifted peaks. 

More could be said about this corner effect if each 
antenna size measured could be examined for the re- 
lationship between current minima and bends in the 
antenna. In some cases such a combination of current 
minima and current maxima at the corners will occur 
that all three loop configurations will have distributions 
shifted essentially by the same amount, and no ap- 
preciable separation in the impedance spirals will be 
noticeable. It would be hard to say anything quantita- 
tive about small effects which are of the same order of 
magnitude as the experimental error, so only the an- 
tenna size for which the corner effect is greatest has 
been considered. 


Comparison with Theory 


In Fig. 10, the experimental input impedances of a 
set of circular loops are compared with Storer’s theo- 
retical results. No attempt has been made in the theory 
to consider the actual driving conditions, but, never- 
theless, fairly good agreement exists. A plot of the ad- 
mittance in Fig. 11 is perhaps more revealing. The con- 
ductances appear to be in good agreement, but some 
discrepancy exists in the susceptance terms. By adding 
a constant negative capacitance to the theoretical 
susceptance, the effect of coupling between line and 
termination can be taken into account. An exact de- 
termination of this quantity has not been made, but a 
qualitative estimate using data available from the cor- 
responding problem in the dipole suggests that the 
correction involved would promote better agreement 
between theoretical and experimental susceptance 
curves. 

A comparison with theory of amplitude and phase of 
the current on a circular loop of particular circumfer- 
ence is provided in Fig. 14. There is some scatter of 
experimental points, but in general good agreement with 
theory exists for both amplitude and phase. The curve 
of a cosine distribution is superimposed and it appears 
to be a reasonable approximation of the actual distri- 
bution. 


General Observations 


If the curves of Figs. 7, 8, 9, and 10 are observed near 
resonance, an estimated bandwidth of 150 mc with 
standing-wave ratio <1.5 can be calculated. (For this 
estimate a constant antenna length was chosen and re- 
sistance was plotted as a function of frequency.) 

The current amplitudes for the large loops of Fig. 13 
show a similarity to a traveling wave form for <90°. 
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Fig. 15—Assumed consinusoidal current distribution (2h =2)). 
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Fig. 16—Current distribution—comparison of three configura- 
tions (complementary slot methods). 


The current at points far from the short circuit re- 
sembles that of an infinitely long transmission line. A 
similar observation has been made by Storer‘ from his 
theoretical solution of the current distribution on large 
loops. 
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Systematic Errors Caused by the Scanning of Antenna Arrays: Phase 
Shifters in the Branch Lines* 
L. A. KURTZ} anp R. S. ELLIOTT} 


Summary—By choosing a suitable equivalent circuit representa- 
tion for an array-type scanning antenna with the phase shifters in 
its branch lines, a general expression is derived for the radiation 
pattern in terms of the active impedances of the radiating elements, 
the incremental phase shift between elements, and the desired 
aperture distribution. If the active impedances of the radiating 
elements vary with beam position as a result of mutual coupling, or 
if the active impedances are constant, but different from the character- 
istic impedances of the branch lines, then an infinite series is required 
to represent the radiation pattern. The first term of the series is the 
desired pattern and the remainder can be defined as the systematic 
error. Individual terms of the series represent beams with relative 
angular positions which correspond to multiples of the interelement 
phase shift and with relative magnitudes which are dependent upon 
the deviation from a match of the active impedance of the radiators. 
The systematic error, if uncontrolled, can prevent the achievement 
of low sidelobe level. Experimental evidence is presented in support 
of the analysis. 


INTRODUCTION 


HE RESEARCH reported in this memorandum 
ose initiated by experimental findings with a 

two-dimensional scanning array. This antenna 
was fed by a standing wave through shunt coupling to 
a main line and the beam was scanned away from 
broadside by inserting phase shifters in the branch lines. 
It was noticed that the pattern steadily and seriously 
deteriorated with increasing scan angle, and developed 
high, unwanted sidelobes. The hypothesis was made 
that this deterioration was caused by the changing 
mutual coupling between elements in the array. A 
theory was evolved, based on this hypothesis, which 
satisfactorily explained the principal perturbations in 
the pattern and linked pattern degradation to variations 
in the branch line impedances. 

To test the theory, a versatile array was constructed 
for which mutual coupling in the plane of scan was 
virtually eliminated. Mutual coupling was then simu- 
lated by a range of controlled impedances achieved by 
varying the number of radiating elements in each stack 
of the array. Adjustable phase shifters were inserted in 
the branch lines, and radiation patterns were measured 
as functions of active impedance and scan angle. Sub- 
sequent sections of this report will be devoted to a de- 
velopment of the theory, discussion of the experiment, 
and a comparative analysis. One significant outgrowth 
of the analysis will be seen to be the specification of a 
tolerable mismatch for a prescribed sidelobe level. 


* Manuscript received by the PGAP, November 25, 1955; revised 
manuscript received, April 16, 1956. This research was performed 
while the authors were members of the Hughes Research Labs., 
Culver City, Calif. 

+ Rantec Corp., Calabasas, Calif. 


THEORY 


The theory to be presented will be limited to arrays 
which are fed by a standing wave in the main line, 
which contain phase shifters only in the branch lines, 
and which are scanned in a plane parallel to the main 
line. The radiating elements can be either slots or di- 
poles; any polarization is permissible, and the coupling 
between the main line and a branch line can be either 
series or shunt. For shunt coupling, an adequate 
equivalent circuit is suggested by Fig. 1. In Fig. 1, the 
junctions are spaced 180 electrical degrees apart and 
the degree of coupling of the mth branch is controlled 
by Am, the turns ratio of the coupling transformer. The 
phase shifters are represented as variable lengths of 
transmission line of characteristic impedance Zo. The 
fixed electrical length of the first branch line is 0; and 
(0,+m@) is the variable electrical length of the mth 
branch line. Thus 6=0 corresponds to broadside opera- 
tion. 


O,+me 


am 
Am — COUPLING TRANSFORMER TURNS RATIO 


| | 
| @ —PHASE SHIFT BETWEEN SUCCESSIVE BRANCH LINES, DEGREES | 
| Z{) IMPEDANCE TERMINATING m'” BRANCH LINE | 
| | 


Fig. 1—Simplified equivalent circuit of scanning array. 


The active impedance terminating each branch line, 
Zr”, may consist of mutual-impedance as well as self- 
impedance. The principal mutual contributions come 
from the two adjacent neighbors. For long arrays with 
conventional aperture distributions, the magnitude of a 
given radiating current, Zr” differs only slightly from 
the neighboring values Ze”! and Ir”*1. If a uniform 
progressive phase shift is employed to scan the beam, 
the phase of Jr” is approximately halfway between that 
of Ipr™! and Ip”t!. Therefore, Zr” does not differ 
greatly from one branch line to another, if the self- 
impedances are alike, which is usually the case. This 
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argument becomes weak for the branch line impedances 
near either end of the array, but for long arrays the 
influence of the end elements is relatively slight. Hence, 
the assumption is made that Zr” is the same for all m 
and the superscript will be dropped in the ensuing 
analysis. 

The additional assumption is made that all system 
components are perfect, with the one exception that 
Zr may not necessarily be equal to Zo. Thus all junc- 
tions are assumed precisely placed, all transformers 
exact, all phase shifters lossless and perfectly matched. 

Let E, be the main line voltage common to all junc- 
tions, and let Jp” be the current flowing through the 
active impedance, Zr, of the mth branch line. Using 
conventional transmission line theory, the relation be- 
tween current and voltage is found to be 


AyEs 
Zr cos (8; + mé) + 7Z, sin (6; + mé) 


ili = 


(1) 


If d is the spacing between radiators and ¢ is the 
azimuthal angle measured from broadside, the array 
factor in the plane of scan is proportional to F(¢), 
where 


M 
F(¢) = yD Tpmeimka sin ¢ (2) 


m=1 


and M is the number of branch lines in the array. 
If the reflection coefficient, I’, is defined by the rela- 
tion 


Zr Ta Zo 
is , 2 (3) 
Zr + Zo 
(1) may be rewritten 
ZEA “fp Ee 3(91+ m8) 
Howe = ; (4) 


ZR ata Zo 1 + Te7 27 (81+ mo) 


By expanding the factor [1+Te-2/@+~®]—-1 in a bi- 
nomial series, the expression for Ip” becomes 


2A pide i 
ys, (—1) "DP e-F n+) (1+ m8) (5) 


Ip™ = 
ZR a Zo n=0 


and, therefore, the field pattern is given by 


M 
F(¢) = B yy A meimlka sin ¢—8] 


m=1 


M 
a (Te 724) B Dy A mei mika sin 6—36] 


m=1 


M 
+ (Te??*)2B Sy A meimtkd sin 58] 


m=1 
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or 


M 


F(¢) 2 B ys (—1)"D%e- 720 yy A meimtka sin @(2n+1) 0] (7) 


n=(0 m=1 
in which 
2E,e~ 7% 
B a ae eS 
Zre+t+Zo 

Eq. (7) is quite general even though it was derived 
for the specific case of shunt loading of the feed line 
with the amount of loading being controlled by ideal 
transformers. A similar expression can be obtained for 
shunt loading of the feed line wherein the coupling is 
controlled by quarter wave transformers of character- 
istic impedance, Zo”. In this case, the term (—j)?"*! 
would replace the term (—1)” and Z»)” would replace 
Am. If series loading of the feed line were employed, 
then J,, the common feed line current, would replace 
E,, all else remaining unchanged. Essential features of 
(7) are independent of the type of loading of the feed line. 

Eq. (6) can be interpreted generally in terms of 
multiple reflections within the branch lines. Successive 
terms of (6) are reduced in magnitude by ie the re- 
flection coefficient of the impedance terminating the 
branch lines. Eq. (6) implies that the effective sidelobe 
level cannot be less than |r|. In addition, successive 
terms differ in interelement phase shift by 26, twice the 
interelement phase shift introduced into the branch 
lines. Therefore, the second term in (6) represents that 
energy which travels out along the branch lines, is re- 
flected by the mismatch of Zr, travels back along the 
branch lines, to be reflected again because of the small 
coupling between branch line and main line, travels out 
along the branch lines again to Zr, where it is radiated. 
Similarly, the third term in (6) represents that energy 
which requires two round trips in the branch lines be- 
fore being radiated, etc. Hence, the array factor given 
by (6) represents many individual beams. Successive 
beams arise from the proper aperture distribution but 
have successively smaller magnitudes and spatial po- 
sitions which are successively shifted relative to the 
main or zero order beam. All beams except the n=0 
beam constitute the systematic error. 

Certain special cases are of interest. For instance, 
if [=0(Zr=Z)), then there will be only one beam for 
all values of interelement phase shift. Likewise, if 6=0, 
there will be only one beam for all values of I. In the 
region between these two extremes, (6) is of value in 
interpreting measured radiation patterns in terms of 
the design aperture distribution, the reflection coef- 
ficient, and the interelement phase shift. 


_ | When the phase shift is so small that the principal lobe of the 
first error pattern overlaps the main beam of the desired pattern, 
the definition of sidelobe level becomes tenuous. What is observed 
is a shoulder on the main beam. At broadside, the principal lobes of 
all error patterns coincide with the desired main beam. Although 
the systematic error is still present, it is undetectable and the side- 
lobe level reverts to its design value. 
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EXPERIMENTAL ARRAY 


Experimental verification of (6) requires an array in 
which the interelement phase shift, 6, and the active 
impedances, Zr, can be independently and accurately 
controlled. Such control has been obtained with the two- 
dimensional array of edge slots shown in Fig. 2. The 
main line waveguide contains a linear array of edge 
slots cut to yield a 32-db Dolph-Tschebyscheff pattern. 
These slots are coupled to the branch line waveguides 
as shown. Each branch line contains a standing-wave 
fed array of edge slots, also cut to give a 32-db Dolph- 
Tschebyscheff pattern. 


Fig. 2—Experimental array with dielectric plugs inserted in the 
branch lines to provide the required phase shift. 


. The required phase shift, m@, was introduced into 
each branch line by inserting a block of polystyrene in 
the waveguide. Quarter wave transformers in the form 
of notches were cut in the ends of the polystyrene blocks 
to reduce the discontinuity at the dielectric interface. 
The phase shifting blocks introduced a vswr of less than 
1.10 and the error in phase shift was less than three 
degrees. With all blocks positioned for 6=0°, the E- 
plane and H-plane patterns of the array are shown in 
Fig. 3. The deviations of these patterns from those 
having uniform 32-db sidelobes characteristic of the 
Dolph-Tschebyscheff aperture distribution result from 
the random errors in the construction of the array and 
from the element factor. 

The impedance, Zr, terminating each of the branch 
lines is that of a linear array of edge slots, and since 
these slots are standing-wave fed, Zr can be varied by 
covering groups of slots with metal-foil tape. The range 
of impedances which were used is shown in Fig. 4, 
(next page) where the numbers indicate which rows of 
slots are radiating. The impedance values of Fig. 4 are 
the averages of the sixteen individual branch line imped- 
ances. These values do not include the effects of mutual 
coupling between arrays. This mutual coupling, how- 
ever, has been reduced to a negligible value by chokes 
located between array waveguides. These chokes re- 
duce the cross polarization in the radiation patterns by 
more than 15 db, and the same surface currents which 
produce cross polarization also produce mutual coupling 
between waveguides. 
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Fig. 3—Measured radiation patterns f=9375 mc. 


A measured radiation pattern is represented by the 
product of the array factor given by (7) and an element 
factor. Fig. 5 shows the element pattern (#7 plane) of 
one of the center linear arrays of the two-dimensional 
array of Fig. 2. Out to an angle of +70° this pattern 
represents, within +1 db, all the element patterns 
except the outside halves of the end element patterns. 


MEASURED RADIATION PATTERNS 


Measured H-plane radiation patterns for the range of 
impedances shown in Fig. 3 and for increments of 20° 
in 6 have been recorded, and representative? patterns 
are given in Figs. 6 through 9 (pp. 623-626). Relative 
gains are not shown because all patterns are normalized. 
The beams are numbered to correspond to the terms in 
(7); 2=0 corresponds to the main beam. The primed 
numbers indicate the second occurrence of a given beam 
due to the element spacing of 0.70) in the H plane. The 
best example of this recurrence of a beam is shown in 
Fig. 9 where 0 = 140°. 


2 For a complete set of patterns the reader is referred to Tech. 
Memo. No. 359, Hughes Res. Lab.; May, 1955. 


IRE TRANSACTIONS ON ANTENNAS 


J 


HEH 


AND PROPAGATION October 


LT 
TITY 
wanes 


RADIATIN 


A 


LOTS 


Fig. 4—Average impedance (Zp) of branch line arrays. 
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Fig. 5—Measured element pattern of individual linear array 
of Fig. 2. 


COMPARISONS OF THEORY AND EXPERIMENT 


The positions and relative magnitudes of the various 
beams shown in Figs. 6 through 9 can be predicted 
from (7) through the knowledge of I and @. Both I and 
6 are determinable quantities for the array described 
above. The average values of I’ can be computed from 
the data in Fig. 3 and 6 can be deduced from the physi- 
cal dimensions of the dielectric which is inserted in the 
branch line waveguides. 

Using (7), the position of a beam is given by the fol- 
lowing equation: 


kd sin 6 — (2n + 1)0 = 2p (8) 


with p an integer. A curve of @ vs (2u+1)@+2zxp for 
d=0.70Xo is shown in Fig. 10. A precision measurement 
of beam positions seemed unwarranted because of the 
shifting effect which is a consequence of the partial 
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Fig. 6—Measured radiation patterns vs Zp; 6=20°; f=9375 me. 


superposition of lobes of patterns of different orders. 
However, when beam positions are read from the pat- 
terns of which Figs. 6 through 9 are representative, and 
compared with the values predicted, a remarkably close 
agreement is noticed. This is evidenced by Table I (p. 
627). 

Comparison of the measured magnitudes of the beams 
with the magnitudes predicted from (7) has been made 


by the following procedure. First, the beam magnitudes 
in the measured patterns such as those in Figs. 6 
through 9 were adjusted to remove the effect of the ele- 
ment factor given in Fig. 5. Second, the average differ- 
ence between successive beams in a given pattern was 
determined and entered in Table II (p. 627). As an ex- 
ample, in Fig. 8(D), corrections of 3 db, ¢ db, and 33 db 
are necessary for beams 0, 1’, and 2’, respectively. With 
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Fig. 7—Measured radiation patterns vs Zp; 6=60°; f=9375 me. 


the element factor removed, beam 1’ is 84 db below beam 
0 and beam 2’ is 2(9$) = 193 db below beam 0. Therefore, 
the average value 9 db is entered in Table II for this 
case (9@=100° and Zr=Z,). 

The last step in the procedure was the conversion of 
the average values of Table II from decibel to voltage 
ratios and the comparison of these ratios with the re- 


flection coefficients obtained from Fig. 4. This com- 
parison is made in Table III (p. 627). The agreement be- 
tween the measured reflection coefficient and the reflec- 
tion coefficient inferred from the radiation patterns is 
quite satisfactory. This agreement, together with the ac- 
curate prediction of the beam positions, provides strong 
support for the theory embodied in (7). 
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DISCUSSION 

This study of resonantly spaced arrays with phase 
shifters in the branch lines has revealed several new 
ideas. The most important of these has already been 
mentioned, namely, that the effective sidelobe level 
cannot be-less than |I'|. In many two-dimensional 
arrays, the mutual coupling between branch line radi- 
ators is not so small as it is in the array of Fig. 2, 
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Fig. 8—Measured radiation patterns vs Zz; 0=100°; f= 
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especially if the other polarization is being employed. 
With increased mutual coupling, the active impedances 
of the branch lines are strongly affected by the inter- 
element phase shift. Thus, in many two-dimensional 
arrays, as the phase shift increases so does the reflection 
coefficient. The influence this has on the systematic 
error is equivalent to moving diagonally across Table II 
and indicates a progressive deterioration of the sidelobe 
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Fig. 10—Beam position vs phase shift. 


level as the beam is scanned further away from its 
central position. Several ways have been found to over- 
come this defect. One method places the phase shifters 
in the main line whereas another maintains the phase 
shifters in the branch lines but uses a traveling wave to 
excite the branch line radiators instead of a standing 
wave. Both of these methods have been tried experi- 
mentally and gave excellent results. The sidelobe level 
remained acceptably low over a wide angle of scan. The 
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TABLE I 
Progressive Beam Beam race Progressive Beam Beam Position 
Phase Shift No. Phase Shift No. @ 
e u Calculated Measured 0 Z Calculated Measured 
20° 0 4.50 Sis) 100° 0 23.40 24.1 
1 VST) 14.1 {3 —13.75 —12.9 
2 23.40 DSM 22 Bons 37.8 
3 33015 32.0 3 — 4.55 — 3.0 
4 45.60 48.0 1205 Opes 28.45 Daf i 
40° 0 9.15 10.4 Qi —72.25 | —69.5 
1 28.45 29.0 1} 0.00 — 1.0 
2, 52.50 53} 65) 2ii —28.45 —29.0 
ile —72.25 —67.3 140° 0 Bin (15) 34.0 
é 2} —40.05 —38.2 0: — 60.80 —57.4 
3} —18.51 —19.0 13 13):275 14.8 
60° 0 ikea Zhs) 15e2 Qii — 4.55 — 4.5 
1 45.60 46.4 3iil —23.40 —22.4 
1g —45.60 44.1 160° 0 40.05 St Tf 
2: —13.75 —13..3 Qi —52.50 —52.7 
80° 0 18.50 19a 1: 28.45 26.3 
] Ds, 69.7 2ii 18.50 18.0 
3 —40.05 — 38.0 Zilli 9.15 8.5 
4 0.00 1.0 4iv 0.00 — 1.5 
1i —28.45 —27.5 SY — 9.15 — 9.0 
Dp 9.15 9.8 6vi —18.50 —19.0 
oF 52.50 5200 
TABLE II 
Average Difference of Beam Level—DB 
Zr 
G—207 C= 4Oe 6=60° @=80° HOO 6=120° 96=140° 6=160° Average 
Zi 25 PIN SS) 24.5 31 28 DAR SS 20.5 19 24 
Z2 LieS 16 16 19 Dili) 1525 16 15 17 
Z3 13 12 iP 18) 14 1S 12 Jol 1235) 
LZ4 9 8 8 8.5 9 os) 8 8 8 
Zs 4 Soe) 3 3 Dee) 359 Sho) Shas) 35) 
TABLE III ent from that of the main beam. Eq. (7) is useful in pre- 
Roiearioa Ceca E dicting this direction. It is also possible that the reflec- 
tion coefficient may be inferred from the radiation 
Impedance Pattern pattern in cases where it is otherwise unobtainable. 
Zr Measurements Measurements 
Pig. 4) EO! CONCLUSION 
os en ane For resonantly fed scanning arrays with phase shifters 
Za 0.26 0.24 in the branch lines, the sidelobe level has been found to 
a on feee depend on the interelement phase shift and the imped- 


traveling-wave feeds gave a slightly lower efficiency but 
provided the added improvement of an increased band- 
width. These results will be reported in detail at a later 
date. . 

Several additional features of interest were uncovered. 
For example, contrary to anticipation, the systematic 
error is commonly observed in a direction quite differ- 


ance mismatch in the branch lines. Spurious lobes 
created by multiple reflections in the branch lines are the 
controlling factor in the sidelobe level. A theory has 
been evolved to explain this phenomenon and has been 
corroborated by extensive experiments. These findings 
have a practical application to arrays in which mutual 
coupling causes the active impedances of the branch 
lines to depend on scan angle. For such arrays a serious 
deterioration of the sidelobe level can result with scan 
and other feeding methods are preferred. 


See 


628 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


October 


A High-Performance Conically-Scanning X-Band 
Antenna of Novel Design”* 
J. G. McCANN{ and R. J. STEGENT 


Summary—A conically-scanning antenna is described which con- 
sists of a scanner mechanism, rear feed, and 8-foot diameter parab- 
ololoid producing a pencil beam having equal E- and H-plane 
beamwidths. The polarization may be quickly changed from linear 
to circular to vertical polarization. The sense of the polarization is 
readily changed from left-hand to right-hand. The replacement of one 
section of transmission line by a half-wave plate changes the unit 
from a nutating polarization scanner to a twice-scan-speed rotating 
polarization scanner. The scanner incorporates some novel features 
such as a symmetrical rear feed having equal E- and H-plane pat- 
terns, high-power quarter-wave and half-wave plates, and an orthog- 
onal mode absorber. 

This paper describes the techniques which were employed to 
obtain the above mentioned performance. The problems encountered 
and their solutions are described. 


INTRODUCTION 
HE DESIRABILITY of maintaining accurate 


tracking information on missiles and aircraft has 

resulted in a demand for better radar sets. These 
improved radar sets must function as both beacon and 
skin trackers with improved sensitivities and accuracies. 
Modified or newly designed antennas which are a part 
of the new radars must handle higher powers and be 
more versatile than formerly available types. The im- 
proved performance of nutating polarization antennas 
over the rotating polarization types in general dictates 
their use. Circular polarization and also other types are 
quite desirable for certain applications. Vibration free 
scanning is most desirable so that units which are 
dynamically balanced must be used. 

The antenna to be described here operates from 9100 
to 9600 mc and includes two interchangeable circularly 
symmetric feeds which radiate linear or circular polar- 
ization equally well. The transmission line contains 
components which may be readily adjustable to provide 
for either horizontal, vertical, right-hand circular, left- 
hand circular, or twice can speed rotating linear polar- 
ization. The complete antenna system is designed to 
handle 250 kw peak power unpressurized. Round wave- 
guide is used throughout most of the transmission line 
so that a rotating feed may be used. Because the moving 
parts rotate rather than nutate, the mechanical com- 
ponents are simplified and near perfect dynamic balanc- 
ing is achieved, thereby minimizing vibration and re- 
ducing the attendant “noise.” 


DESIGN 


General 


The conically scanning antenna consists of a wave- 
guide transmission line, torque tube, feed, scanner as- 


* Manuscript received by the PGAP, October 3, 1955; revised 
manuscript received February 23, 1956. 
7+ Canoga Corp., Van Nuys, Calif. 


sembly with drive motor, and reference generator. The 
complete unit is shown in Figs. 1 and 2 and is to be 
mounted in an 8-foot diameter paraboloid whose focal 
length is 27.5 inches. The transmission line and torque 
tube include taper sections, orthogonal mode absorber, 
quarter- and half-wave plates, and rotary joints. On 
transmitting, the rf energy enters the rectangular 
waveguide at the rear of the unit; passes, in order, 
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Fig. 1—Cutaway drawing of the complete unit and paraboloid. 


Fig. 2—Conical scanning unit including two feeds, half-wave 
plate and pressurizing hood. 


through a taper section from rectangular to round 
waveguide, an orthogonal mode absorber, a rotary 
joint, a quarter-wave plate section, and a second rotary 
joint. At this point the rf either passes through a half- 
wave plate section or a hollow waveguide, and then 
through a taper section to dielectric loaded round 
waveguide. A splash plate directs the energy through a 
pressurizing cover and onto the paraboloid which colli- 
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mates it into a beam. The circularly symmetric splash 
plate of the feed is offset from the axis of the paraboloid 
so that the center of radiation is displaced from the 
axis of rotation. This displacement causes the center of 
radiation to move in a circle about the focal point of the 
reflector, the radiated beam of the antenna describing a 
corresponding circular motion in space. By using either 
of two feed units one may obtain a 3-db or 1-db cross- 
over. The complete unit including the radome may be 
pressurized to 30 psi absolute. This technique eliminates 
the need for a rotating pressure seal. The various types 
of polarizations are obtained by simple adjustments of 
parts of the unit. 


Polarizations Available 


The rf input waveguide, taper section, and orthogonal 
mode absorber may be rotated to any one of four po- 
sitions. Vertical polarization is obtained by positioning 
the rectangular waveguide so that the electric vector is 
vertical. The quarter-wave plate section is horizontal 
and, therefore, is normal to the electric field and conse- 
quently has no effect on the linearity of polarization. 
Rotating the input waveguide 90° results in horizontal 
polarization. The quarter-wave plate section being 
parallel to the electric field again does not effect the 
linearity of the polarization. 

If the rectangular waveguide is rotated to a position 
45° from the vertical, the polarization will be circular. 
Right-hand polarization! results when the electric 
vector, as viewed from the rear of the unit, is up to the 
left; left-hand, when the £ vector is up to the right. 

Linear polarization rotating at twice the scan speed 
may be achieved by use of a half-wave plate rotating 
with the feed. The rf input waveguide is positioned so as 
to result in linear polarization leaving the quarter-wave 
plate. This linear polarization, on propagating through 
the rotating half-wave plate, will be converted to linear 
polarization rotating at twice the scan rate. It appears 
that twice scan rate rotating polarization may provide 
better tracking of a linearly polarized beacon than 
ordinary rotating polarization. 


MICROWAVE COMPONENTS 


Input Section 


The input section consists of an inductive iris, a taper 
section from rectangular to round waveguide, an orthog- 
onal mode absorber, and one-half of a rotary joint. 
The function of the inductive iris is to match out the 
inherent discontinuity of the taper section and orthog- 
onal mode absorber. The maximum vswr of the input 
section over the band 9.0 to 9.7 kmc without the iris 
was 1.20 and with the iris was 1.08. 

If it is assumed that a maximum value for the axial 


1 The sense of polarization is defined by IRE Standards, for an 
observer watching the departing wave: The polarization is “right- 
handed” or “left-handed” according to whether the £-vector traces 
out the ellipse in the clockwise or counterclockwise sense, respec- 


tively. 
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ratio is 0.6 db and the antenna vswr is under 1.4 then 
it is easily shown that the vswr of the orthogonal mode 
absorber must be less than 1.5. The absorber must 
dissipate all of the reflected power from the feed which, 
for 250 watts average power and a vswr of 1.4, would 
be 10.5 watts. After considering the relative advantages 
and disadvantages of a number of different types of 
absorbers, it was decided that a polyiron loaded slot 
should be used. Measurements indicated that a long 
axial polyiron loaded slot in round waveguide would 
have very low loss (less than 0.1 db for 12 inches of slot) 
and reflection when oriented so as not to couple to the 
incident wave. The polyiron is placed on the outside of 
the round waveguide, thereby minimizing attenuation 
of the incident wave. When oriented for maximum 
coupling the reflection from the device was still very 
small. The attenuation was about 1.1 db per inch. Two 
slots each 3.5 inches long and diametrically opposite 
each other are used. The total one-way attenuation is 
7.7 db which represents a vswr of 1.46 to the reflected 
wave. The polyiron is in contact with the brass wave- 
guide and can therefore adequately handle the power 
to be dissipated. 

The rotary joint is a half wavelength folded coaxial 
type in which the break in the metal occurs at a current 
minimum. The characteristic impedances of the two 
quarter-wavelength sections are designed to give broad- 
band performance. 


Quarter-Wave Plate Section 

The quarter-wave plate consists of a length of guide 
having circular cross section at either end and departing 
from circular symmetry in the intermediate section. 
This part was broached from a thick walled tube using a 
special broach having a circular cross section with a flat 
side. After broaching, the ends of the tube were turned 
to a circular cross section. The ends of the flat portion 
were turned to a diameter less than the tube diameter 
to provide a matching section. The length of the match- 
ing section was determined experimentally, being about 
a quarter guide wavelength. The length of the center 
section was adjusted to obtain 90° of relative phase 
shift between the two orthogonal modes normal to and 
parallel to the flat portion. The axial ratios as a function 
of frequency for three different length quarter-wave 
plates are shown in Fig. 3. 


Half-Wave Plate Section 


The design and fabrication of the half-wave plate is 
similar to that of the quarter-wave plate. The center 
section is made longer so that there is a differential phase 
shift of 180° between the principal orthogonal modes. 
Fig. 4 shows the axial ratio as a function of frequency. 


Feed Unit 

The feed unit shown in Fig. 5 consists of a section of 
air-filled round waveguide, a taper section to teflon- 
loaded waveguide, and a teflon supported metal splash 
plate. 
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Fig. 3—Axial ratios for three different length plates. 
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Fig. 5—Section view of feed. 


SERRATED SPLASH PLATE 


The first section of air-filled guide is cut at a small 
angle with the normal to the axis of the waveguide. 
Rotation of the torque tube causes the center of radia- 
tion to describe a circle about the axis of the paraboloid. 
The collimated beam from the paraboloid traces out a 
similar path in space. 

The teflon taper begins at the outer wall of the larger 
diameter round waveguide for several reasons. This 
construction, in conjunction with the splash plate sup- 
port, is used to hold the complete feed unit together. 
The outside taper also reduces the fields at the air- 
dielectric interfacial space, thereby reducing the possi- 
bility of corona or arc-over. The length of the taper is 
adjusted so that the mismatch due to incremental 
changes in cross section are minimized. 

The teflon loaded round waveguide was designed to 
have the same cutoff wavelength as the air-filled wave- 
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guide. Teflon is used to reduce the diameter of the ring 
source as well as for strength. Since a ring source has 
different centers of radiation in the & and /Z planes it 
cannot be focused properly in a paraboloid. This de- 
focusing results in poorer beam formation than with a 
point source feed. The diameter of the teflon loaded 
waveguide is made as small as possible to minimize this 
defocusing. 

The radiating portion of this feed is a modification of 
a Naval Research Laboratory design.? The teflon in 
the waveguide is extended beyond the énd of the tube 
and supports the serrated splash plate. The depth of the 
serrations may be adjusted to control the H-plane 
primary patterns. These serrations have little effect on 
the #-plane pattern, the outer diameter of the splash 
plate determining its shape. The match of the feed to its 
transmission line has been optimized over the band by 
properly shaping and positioning the various parts. The 
matching of the feed was obtained chiefly by adjusting 
the distance from the splash plate to the end of the 
metal tube and by introducing additional reflections 
from the flat at the center of the splash plate. These re- 
flections are at a position which tends to cancel other 
reflections over the complete band. The splash plate 
support and the waveguide dielectric are threaded and 
screwed together. Teflon pins are inserted axially to lock 
the parts together. The splash plate retains these pins, 
the splash plate being held in by a teflon nut. 


Radome 


The pressurizing hood is a thin wall type which has 
been designed to allow the wave to pass through at or 
near normal incidence. The hood is fabricated of a fiber- 
glas cloth impregnated with a low-loss polyester resin 
and will withstand 30 psi pressure. The radiation pat- 
tern is only very slightly affected by the presence of the 
hood. The signal loss at the peak of the main beam due 
to the hood is 0.5 db. 


Polyiron Absorbing Disk 


The polyiron disk shown in Fig. 1 absorbs any power 
that may tend to leak into the region of the front bear- 
ing, thereby minimizing the possibility of bearing burn- 
out due to rf arcing. 


Secondary Radiation Patterns and VSWR and Power 


The maximum vswr of the complete unit for twice 
scan speed rotating linear polarization is 1.32. The 
maximum vswr’s when radiating nutating linear and 
circular polarizations are 1.32 and 1.17 respectively. The 
vswr of the circularly polarized feeds are obviously 
better than the linearly polarized units because reflected 
energy from beyond the quarter-wave plate is returned 
as an orthogonal mode and dissipated in the absorber. 

The measured beam crossover levels for the two dif- 
ferent offsets are 3.2 +0.4 db and 1.0+0.2 db. The offset 


2 F. L. Hennessey, “A Broad Band Axially Symmetric Vertex 
Feed,” Naval Res. Lab. Rep. No. 4071; October 22, 1952. 
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Fig. 6—Offset plane patterns showing beam crossover (frequency 
9100 mc). 


plane secondary radiation patterns are shown in Fig. 6. 
All side lobes were below 25 db. The patterns were 
measured with a transmitter-receiver distance of 1600 
feet. The complete unit was high power tested without 
pressurization. There was no sign of arcing or corona at 
the maximum available power of 250 kw. 


Drive System 


The conical scanner is driven by a synchronous motor 
rated at 4 hp, 115 v, 60 cycles, 3 phase, 1800 rpm. This 
drive motor is connected by means of a Gilmer “Tim- 
ing” belt to the choke and sprocket* assembly of the 
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torque tube section as shown in Fig. 1. The use of this 
belt drive system insures quiet, vibrationless operation. 
The need for periodic lubrication, oil and grease seals, 
and other difficulties associated with gear drives is elimi- 
nated. The belt and sprockets have teeth which elimi- 
nate slippage and insure that conical scanning occurs at 
the synchronous motor speed. An idler pulley is pro- 
vided for adjusting belt tension. 


Reference Generator 


The reference generator is a two-pole generator de- 
livering 67 v into an 1800 ohm load. It may be elec- 
trically phased by loosening base clamping screws. 


CONCLUSION 


A high-performance conically scanning antenna has 
been designed for use in conjunction with a high-pre- 
cision, high-power, X-band tracking radar. The antenna 
provides a variety of different types of polarization, 
with two different beam crossover levels. The feed ele- 
ments rotate and therefore may be completely dynam- 
ically balanced, thereby providing vibration free opera- 
tion. Highly precise tracking is achieved because of the 
reduction in error due to vibration of the antenna. 
High power components are used throughout and it is 
believed that this unit will handle one megawatt when 
pressurized to 30 psi absolute. 
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Slot Admittance Data at K, Band* 
MAURICE G. CHERNINt 


Summary—Admittance data on transverse edge slots in RG-96/U 
waveguide can be obtained by a technique called the moving lossy 
short technique. By this technique the radiation attenuation of a test 
section of identical slots can be determined. It is then possible to 
specify both the slot inclination angle and the slot depth of cut re- 
quired to yield any conductance within the range of measurements. 
These data were used to design an experimental 30-slot array with 
a 25-db Taylor aperture distribution. This array was successful, 
yielding sidelobes near —23 db. Subsequently, the same data were 
used in designing an 8-foot array with 432 edge slots having the same 
aperture distribution. This array had a half-power beamwidth of 
14 minutes and sidelobes of the order of —24 db. These results com- 
pared favorably with design objectives. 


INTRODUCTION 


ECENT developments in system components at 
R K, band have made the need for antennas in this 
frequency band more pronounced. The long, 
linear array with transverse edge slots is an antenna 
type which has had successful use at the lower fre- 
quencies, for example at X band. To construct similar 
linear arrays at K, band, it is necessary to obtain in- 
formation on the electrical characteristics of the slots 
which are to be used. Such information can be obtained 
by scaling data obtained at X band, although this 
method is not entirely satisfactory due to the dimen- 
sional differences of the waveguide used at the two fre- 
quencies. A second method is to measure the slots at 
K, band as they are to be used. Since this method does 
not require the extrapolation of data from other fre- 
quencies, it offers a more direct approach to the 
problem. 

An adaptation of a measuring technique suggested by 
Oliner! was used to measure the radiation and dissi- 
pative attenuation of sections of K,-band waveguide 
with identical transverse edge slots in one of the narrow 
walls. This technique, called the moving lossy short 
technique, has as its principal advantage over the 
standard slot measurement procedure the fact that the 
vswr which is being measured can be adjusted over a 
wide range by varying the loss in the movable short. 
The characteristics of the unknown to be measured de- 
termine the optimum amount of loss to be incorporated 
in the short. Another advantage over the standard slot 
measurement procedure is that the moving lossy short 
technique “automatically” averages 8 to 10 (or more) 
data points, compared with the one datum point per 
frequency which is usually taken with the standard 
technique. 


* Manuscript received by the PGAP, October 3, 1955; revised 
manuscript received, April 18, 1956. 

+ Hughes Aircraft Co., Inc., Culver City, Calif. 

1H. M. Altschuler and A. A. Oliner, “Microwave measurements 
with a lossy variable short circuit,” Res. Rep. R-399-54, Polytech. 
Inst. of Brooklyn, Bklyn., N. Y., September 16, 1954. 


It should be noted that the determination of the at- 
tenuation of a slot or group of slots requires no knowl- 
edge of the admittance of the slot. Hence, it is not 
necessary to know if a slot is resonant to determine the 
frequency at which radiated power is maximized. How- 
ever, since the conductance of a slot is a more common 
reference term, the data are presented in this form. 


THEORY 


The philosophy behind the moving lossy short tech- 
nique is that the attenuation produced by a test section 
is related to certain of its scattering matrix elements and 
can be determined more precisely from the scattering 
matrix elements than directly. The total attenuation 
(A,) of a test section is equal to the dissipative attenua- 
tion (Aq) plus the reflective attenuation (A,). The dis- 
sipative attenuation includes both radiation and wall 
loss. If the assumption is made that the test sections are 
symmetrical, then these attenuations can be expressed 
in decibels as follows: 


A, = — 10 logio bea (1 — | Sul) | 
Su? 
A, = — 10 log io te |r. |? 
Leaks 
Ag = —=10 logio fe | a : (1) 
; ; O11 


These equations can be developed in the following 
manner. Network theory defines the matrix element Sj. 
as the transmission coefficient from terminal 1 to ter- 
minal 2. 


Probe 


CIN EZ 
Slotted Line |Test Section Z 
1 a 


ss] Movable Lossy Short 


The insertion loss is then equal to —20 logio | Sie in 
decibels. However, when there is loss present in the 
short circuit, one half of this loss appears in the insertion 
loss as determined from | Sia . Hence, to determine the 
insertion loss of the given test section only, the effect of 
the lossy short must be accounted for; that is, 


A, = total insertion loss — 1/2 loss in termination 


= — 20 logio | Sis | + 10 logio | its 


is the reflection coefficient of the short. 


where | If 
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| Sul? 
ON 


| s 


Aj = 10 log io 


From Deschamps’ construction technique? 


Sees 
Heer Sarl 
From physical symmetry 
[Su] = | Sr]. 


Substituting for [Sis] in the expression for A, 


A; — 10 logio bea (1 cae [sul |. (2) 


: | 

This equation is the expression for the total attenua- 
tion of a test section and represents attenuation due to 
radiation, dissipation, and reflection. The reflected 
power from the discontinuities of the slots is propor- 
tional to the square of the reflection coefficients. Since 
the reflection coefficient of a discontinuity is the first 
element in its scattering matrix, the reflected power 
from the test section is proportional to |Su|2. The 
power that is not reflected is 1—| S|”. The attenuation 
of the section due to reflection is 


— 10 logy [1 — | Su 


A, = 10 logis — APEC) 


RRSP, 


Since the reflected power is the entire reflected power 
(that is, includes power reflected from the lossy short), 
(3) for A, must be corrected. 

If the reflection coefficient of the test section is desig- 
nated as kale then? 

| Sul? 


|T 


9 


/ 
ees 


) 


and 


j Sur]? 
(4) 


A, = — 10lo = 
«| |r, [2 


The equation for the dissipative attenuation may be 
obtained by algebraic manipulation of (2) and (4) in 
A,=A,+Aa. 


EXPERIMENTAL PROCEDURE 


To apply the moving lossy short technique, numerous 
test sections were manufactured, each with twenty 
identical slots. The sections differed from each other in 
slot inclination angle or in depth of cut. The slot spacing 
was chosen in such a manner that the same admittance 
data might be used in the design of a long array which 
would meet small “squint angle” requirements. The 
waveguide wall in which the slots were cut was milled 
to a 0.015-inch thickness before the slots were cut. The 


2G. A. Deschamps, “Determination of reflection coefficients and 
insertion loss of a waveguide junction,” J. Appl. Phys., vol. 24, pp. 
1046-1050; August, 1953. 

3H. M. Altschuler and A. A. Oliner, op. cit., 
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milling made it possible to control the wavelength at 


’ which a test section displayed its maximum attenuation 


by varying the slot depth of cut. 

Two sets of data were taken at each frequency for 
each test section measured. The first measurements 
were taken with the movable lossy short connected 
directly to a slotted line. The test section to be measured 
was then inserted between the slotted line and the lossy 
short. For both measurement runs, the movable lossy 
short was moved to eight equally-spaced positions so 
that the term S, (an element of the scattering matrix 
of the test section) could be obtained. The position of 
the voltage minimum and the vswr in the slotted line 
were recorded for each position of the short. The eight 
positions of the short covered a movement of exactly 
one-half guide wavelength. For optimum accuracy, the 
guide wavelength was measured by the short driving 
mechanism with the slotted line probe fixed. 

From these data, the amplitude and phase of the re- 
flection coefficient, [, were calculated for each position 
of the lossy short and plotted on standard polar co- 
ordinate paper. When the calculations for the data ob- 
tained with the test section absent were plotted, a circle 
resulted which has its center at or very near the origin 
of the plane. If the circle is centric, the radius is equal to 
ies i This circle is, in effect, a calibration of the lossy 
short at that particular frequency. The corresponding 
calculations when the test section was present in the 
setup yielded a circle in the reflection coefficient plane 
which has a radius and displacement from the center of 
the plane related to both the power being radiated by 
the slots and the power being dissipated in wall losses. 
The radius of this circle is p. 

The term |.S;|, an element of the scattering matrix 
of the test section, is then obtained by geometric con- 
structions on the [’-plane circles resulting from the data 
recorded with the test section present. The method of 
constructions is illustrated in Fig. 1. The underlying 
theory, as well as the procedures for obtaining the other 
elements of the scattering matrix (amplitude and phase), 
is presented by Deschamps. 

Through the use of (1) through (4), the measure- 
ments made on each test section were translated into 
attenuation per section. The attenuation in turn was 
converted into conductance per slot in this manner. 

Each test section is considered as 


Pr 


> i 


ay 


4 Equal spacing of the measurement positions is not necessary 
when the movable lossy short is connected directly to the slotted line 
(test section absent). However, in the experimental work the posi- 
tions used were equally spaced for both sets of measurements because 
the data taking was facilitated by maintaining the same test setup 
for all runs. 


634 


Fig. 1—Method of construction of T-plane circles. 

(a) 1) Plot reflection coefficients on polar coordinate paper 
(always an even number of points) and draw best circle 
through all points. 

2) Number points to give two equal groups of points. 
3) Label center of circle C and intersection of polar coordi- 
nates O. 
4) Connect points 1 with 1’, 2 with 2’, etc. 
5) Label intersection of chords O. 
(b) 1) Connect C and O. t 
2) Raise perpendiculars at C and O to intersect circle periph- 
ery at opposite sides. Label one point of intersection A, 
and the other B. 

3) Draw AB. 5 

4) O’ is the intersection of COand AB. 

5) Draw OO’. Line OO’ is Su. 


where 
p;=power input to section, 
Po=power output from section, 
pb, =power radiated from section, 
pi=power dissipated in wall losses. 


From an energy conservation outlook, it may be said 
that 


P= ipe apie Pr 


or 
(oe 
Pi Pi pr 
Ley eee (5) 
Pi Pi Pi 
where 
Pp Bit 
ered (6) 
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and B is the wall loss attenuation factor of the waveguide 
without slots (expressed as a power ratio). Eq. (5) may 
be applied to elements of the test section, each element 
extending from midway between two slots to the next 
adjacent midway position as shown below. 


<> 
| 


Since g, the conductance per slot, is equivalent to the 
ratio of the power radiated by a slot to the power in the 
guide at that point, the value of g can be calculated from 
(5) if the ratios p/p; and pi/p; for each slot and its 
associated spacing are known. The test section factor 
p./ pi, or rather its inverse, is the dissipative attenuation 
and is determined from the measurements. It is ex- 
pressed in decibels and represents the dissipative atten- 
uation of twenty slots and their associated slot-to-slot 
spacings. Dividing this attenuation figure in decibels by 
20 and converting it to a power ratio make it suitable 
for insertion into (5) for the calculation of conductance 
per slot. The wall loss factor, p:/p;, the remaining term 
which must be determined, can be calculated from (6) 
if 8 is known. This term can be read from a handbook 
and converted into the expression for one-slot spacing. 
The conductance per slot, g, can then be calculated 
from (5). 


RESULTS 


For each test section a set of circles in the reflection 
coefficient plane was plotted over a frequency range 
from 8.2 to 8.9 millimeters. Since each test section was 
unique in slot inclination angle or in slot depth of cut, 
complete information was acquired on the effect of 
these parameters on the admittance characteristics of 
slots at K, band. From these [’-plane circles “resonance” 
graphs similar to that shown in Fig. 2 were constructed 
with conductance per slot plotted vs wavelength for 
both fixed slot inclination angle and different slot depths 
of cut. The information from these “resonance” graphs 
were compiled in a “summary” graph to show the 
manner in which the conductance per slot vs wavelength 
varied as both the slot inclination angle and depth of cut 
were changed. Finally, it was possible to obtain two 
curves from which the slot inclination angle and depth 
of cut could be specified to yield a certain conductance 
at a particular frequency. Once an array has been de- 
signed and the individual slot conductances determined 
which are required to produce the desired aperture dis- 
tribution, then the necessary slot angles and depths of 
cut may be read off the pertinent graphs. 

A 30-slot, nonresonant linear array which would yield 
a —25 db sidelobe level was constructed to test the data. 
Transverse slots were cut in the edge of RG-96/U 
waveguide and spaced 0.2424 inch apart. The slot wall 
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Fig. 2—Sample “resonance” graph plotted from I-plane circles with 
conductance per slot vs wavelength for slot inclination angle 


of 10 degrees and slot depths of 0.025, 0.028, and 0.030 inch. 


Fig. 3—432-slot array. 


was milled to a thickness of 0.015 inch. The array was 
designed to deliver 30 per cent of the input power to the 
load. Radiation patterns of the array, taken with a horn, 
show 23-db sidelobes, close enough to the design objec- 
tive to satisfy the experimental requirements. The vswr 
of this array was under 1.12 over the entire frequency 
band of 8.3 to 8.8 millimeters. 

An 8-foot, nonresonant array with 432 transverse 
edge slots was subsequently designed as a line source. 
It had been decided that this array was to have a 25-db 
Taylor’ aperture distribution, and 5 per cent of the array 


5 T. T. Taylor, “Design of line sources for narrow beamwidth and 
low sidelobes,” Tech. Memo. No. 316, Res. and Dey. Labs., Hughes 
Aircraft Co., July 31, 1953. 
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Fig. 4—Electrical characteristics of 432-slot array. (a) VSWR. 
(b) £-plane radiation pattern. 


input power was to be delivered to a matched load. As 
with the experimental 30-slot array, the slot wall was 
milled to a thickness of 0.015 inch. This array is shown 
in Fig. 3 in a support on the pattern range. The vswr 
vs wavelength and the #-plane radiation pattern are 
shown in Fig. 4. The half-power beamwidth is approxi- 


636 


mately 14 minutes, as measured with a precision rotary 
table. Both the sidelobe level of the array and the half- 
power beamwidth compare favorably with the design 
objectives. 

In the interests of improved accuracy, modifications 
had to be made to existing laboratory and machining 
equipment. In the measuring setup, the slotted line was 
modified to enable the position of the carriage to be 
measurable as close to 0.0001 inch as possible. This 
accuracy was accomplished by providing a 40-thread- 
per-inch drive with a 0.0001-inch least count Ames 
gauge attached to the carriage. The moving lossy short 
was also modified for a more controlled positioning of 
the short by replacing the standard micrometer drive 
with a unique driving mechanism plus an Ames gauge 
with a least count of 0.0001 inch. 

Machining problems arose as a result of the tight 
tolerances necessary at such high frequencies as those 
in the A, band, where, for example, a 0.0012-inch error 
in slot spacing can cause a phase error of 1 degree. 
Precision waveguide milled from coin silver (90 per cent 
silver and 10 per cent copper) was used in preference to 
brass waveguide with an internal silver lining, the only 
other waveguide available for the specific frequency. 
The silver lining tended to delaminate during milling 
operations. The actual cutting of the slots to tolerance 
with the exact design spacing between them was ac- 
complished through the use of a specially designed 
holding fixture. The device forced a precision-milled 
gib against one of the wide dimensions of the guide and 
enabled it to be stressed equally during the cutting 
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process (Fig. 5). The slots were undercut several 
thousandths of an inch. To finish the slots, the guide 
section was mounted on the table of a precision travel- 
ing-microscope and the slots were brought into tolerance 
and into shape (sharp, square corners) by hand with the 
aid of small, specially designed files. 


Fig. 5—Slot milling operation. 
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Radiation by Disks and Conical Structures* 
A. LEITNER} anv C. P. WELLSt+ 


Summary—The Lebedev integral transform is applied to a class 
of boundary value problems in the theory of diffraction and anten- 
nas, including circular disks, apertures and hollow conical structures. 
It is found that the conventional Wiener-Hopf technique, together 
with this transform, does not explicitly solve these problems. Instead, 
one is led to an infinite system of linear equations for the representa- 
tion of the unknown transform function. 


INTRODUCTION 


N THIS PAPER, we study the use of the method 

of integral transforms in a class of radiation and 

antenna problems involving mixed boundary con- 
ditions on coordinate surfaces of the spherical system. 
The radiating objects of principal interest in this class 
are the hollow circular cone and the biconical antenna, 
both of finite slant height and arbitrary apex angle. If 
these electromagnetic radiation problems possess azi- 
muthal symmetry, they are reducible to scalar form. 
Let 4 be the half-apex angle of the cone. The problems 
can be formulated as two-part boundary value problems 
on the surface 6 =@. It was hoped that the use of the 
appropriate integral transform theorem combined with 
the Wiener-Hopf technique could produce an explicit 
solution as was demonstrated in the electrostatic prob- 
lem of the conducting cone.! Here we use the word 
“explicit” in the sense that the coefficients in the series 
of wave functions representing the solution can be de- 
termined without recourse to an infinite system of linear 
equations. 

The transition from the electrostatic to the wave 
problem involving conical structures is tantamount to 
transition from the Mellin to the Lebedev transform. It 
appears that this extension is of considerable interest 
among applied mathematicians. 

We find that the method “fails’’ in the sense that the 
conventional Wiener-Hopf technique also leads to an 
infinite system of linear equations, although we find it 
to be of particularly simple form involving factorials 
only. In the following sections we present a test of these 
methods on a simpler acoustic problem involving the 
circular disk. The disk is a cone of half-apex angle 7/2, 
and Karp! has shown that problems involving cones of 
arbitrary angle can be solved by the Wiener-Hopf tech- 
nique if this can be done for the disk; also, the wave 
functions of the azimuthally symmetric electromagnetic 


* Manuscript received by the PGAP, November 30, 1955; re- 
vised manuscript received June 18, 1956. This work was performed 
under sponsorship of the Office of Ordnance Research, U. S. Army. 

Mich. State Univ., Lansing, Mich. : 

1S, N. Karp, “Natural Charge Distribution and Capacitance of 
a Finite Conical Shell,” Res. Rep. EM-35, Math. Res. Group, N. Y. 
Univ., New York, N. Y., 1951. 


field are essentially the same as those for acoustic fields 
of like symmetry. One may conclude that the solution 
for conical structures may produce an infinite linear 
system of a form much simpler than those obtained 
with other techniques. We hope to present these results 
in a later paper. 


STATEMENT OF THE PROBLEM 


We seek solutions of V?-u+k?u =0, R=27/) where d is 
the wavelength and u the velocity potential. We impose 
the radiation condition, r(du/dr+iku)—0, ru bounded, 
as re, 

The problem is that of the radiation by the freely 
vibrating rigid circular disk 


(0u/02Z)o—2/2 = 


(beara = 0, yr 


v, rs a, (1) 
(2) 


The function uw also represents the wave scattered by a 
stationary rigid circular disk from a normally incident 
plane wave. 

This problem has been explicitly solved by use of the 
oblate spheroidal wave functions, (cf. Bouwkamp and 
other authors”). It has also been treated in cylindrical 
coordinates, using the Hankel transform and the exact 
inversion of a certain infinite linear system, by Sommer- 
feld.* 

In spherical coordinates the above function u can be 
represented as a series of the type 


IV 
a 


Ye gurl ’Zngrya( Er) Px(Co8 8) 


n=0 


where Z,41/2 are cylinder functions, P, Legendre func- 
tions and g, unknown coefficients. These series may be 
converted to contour integrals surrounding the positive 
real axis of the complex order plane. The next step is to 
deform this contour into one lying inside an infinite 
strip of finite width surrounding the imaginary axis, but 
it can be shown that this leads to divergent integrals in 
the problem stated above. However, if we use the device 
of transforming kR——v7y, with y real and positive, this 
deformation is permissible without incurring change in 
value of the integrals involved. The transformation 
leads to modified Bessel functions of real, positive argu- 
ment. The Wiener-Hopf technique can now be applied. 


2C. F. Bouwkamp, “Theoretische en numerieke behandeling 
van de buigung door een ronde opening,” Dissertation, Groningen; 
1941. See also: “Diffraction Theory,” Reps. Progr. Phys., Phys. Soc., 
Gt. Brit., vol. 17, pp. 35-100; 1954. 

3 A. Sommerfeld, “Die frei schwingende kolbenmembran,” Ann. 
Phys., vol. 42, pp. 389-420; 1942. 
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The idea is to solve the intermediate problem with pure 
imaginary k, then to reconvert the contour around the 
real axis and, finally, to transform to real positive k. 
This device was first used by Oberhettinger‘ in a similar 
approach to the problem of diffraction by a wedge. 

The transformation to modified Bessel functions 
alters the radiation condition, stated above, to exponen- 


tial decay of the function like e~7"/r as r>~. 


SOLUTION 


We represent (yr, @) as follows: 
r'u(yr, 0) = { ug(u)Payrsu( | cos8| )Za(yr)du, (3) 
ii 


where J_, is the modified Bessel function and g(u) the 
unknown function to be determined by boundary con- 
ditions (1) and (2); ZL is a contour inside an infinitely 
long strip of finite width about the imaginary p axis 
from Im=— © to Imup=+~; P_s+, (cos @) is the 
Legendre function, the appropriate angular function 
when @S7/2, but, when 7/2S0S7, P_124,(—cos @) 
must be used. 
Imposing the boundary conditions (1) and (2): 


~ 
IIA 


eer bie fi Datei ser Olean dual) 
L 


IV 


a: 0 = f ugla)PvalOFsOridu. (5) 
qa 


This is a set of dual integral equations for g(u). Their 
form can be simplified by use of the Lebedev transform 
theorem,? 


Hep { rane Grat (6a) 
16 


ity = if BER aere, (6b) 


where K, is “the modified Hankel” or MacDonald func- 
tion. The theorem is valid provided both integrals con- 
verge and provided A(w) is an even function of 4p, 
analytic in a strip of finite width containing Z and the 
imaginary axis, with decay at least as rapid as 
aA nal 7 | 3/2 as 7, the imaginary part of pw, goes to 
+ « in this strip. 

Using this transform theorem we represent the in- 
homogeneous term in (4) as follows 


ae ent f 1st ord (7a) 
= J wA(u)l_,(yr) du, a 
(Ez we ib : 
4 F. Oberhettinger, “Diffraction of Waves by a Wedge,” Comm. 


Math. vol. 7, pp, 551-563; August, 1954. 

. Lebedev, “On representation of arbitrary functions as 
eee: over MacDonald Functions of Complex Order,” Dokl. Ak. 
Nauk. SSSR (N.S.) vol. 58, pp. 1007-1010; 1947. (In Russian.) 


Lure Reais 


ANTENNAS 


AND PROPAGATION October 


with the result that 


NG) = Grote f “AK (yr)dr. (7b) 


The integral in (7b), which defines A(u), converges in 
the strip —3/2<Reyw<3/2 and has the appropriate de- 
cay such that the conditions of Lebedev’s theorem are 
satisfied; it can be continued outside this strip by the 
use of the Lommel functions. It is found that A has the 
growth of K,(ya)/u as |u| 2 and possesses simple 
poles at w= +(2n+3/2) with the residues 


(21/2iv/ary8!?)(—)"P(n + 3)/n!, a = 0,1,2,---. (8) 


The dual integral equations can now be written 


POs fale) Pavers) — A(u) |Z_u(yr)du, (9) 


x 
IIV 


a: 0 = f gla)P saia()T ord, (10) 


and in this homogeneous form the first step of the 
Wiener-Hopf technique may be applied: 7.e., (9) and 
(10) will be formally satisfied, if it can be demonstrated 
that the integrands are analytic functions of at least 
algebraic decay in halfplanes to the left and right, re- 
spectively, of the strip containing ZL, with these half- 
planes of analyticity overlapping in the strip. 

It is precisely at this point that the method requires 
modifications of the integral equation, because for the 
second of them no such halfplane of analyticity can be 
found in its present form. Since this circumstance is un- 
conventional we present here certain relations on which 
the proof of the assertion is based: The factors of the 
kernel J_,(yr) are, as is well known, the integral trans- 
forms of the boundary values u(yr, 7/2) when ¢ Sa, and 
du(yr, w/2)/dn when r=a. Thus 


wi|g(u) P_1j24y'(0) — A(u) | 


I if r'![du(yr, r/2)/dn|K,(yr)dr, 


a 


wien Pag ewes aii Muy, 2/2) Kylndr 


0 


Although uw(yr, 7/2) and its derivative are unknown, 
their behavior as functions of 7 is well known from 
general diffraction theory, including that of singularities 
near edges, such as at r=a. Thus the functions defined 
by the integrals and their analytic continuations can be 
completely analyzed (the first integral above is an en- 
tire function of w). One can show that g(u) has simple 
poles at the points w= +(2n+4), +(2n+3/2), in fact, 
the residues at the latter set of points are known in 
terms of the residues of A(w) given in (8). Furthermore, 
these integrals determine the asymptotic growth of 
these functions completely as K,(ya)/u* when |u| > 0, 
where a>}. 
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From these results one may derive that the integrand 
in (9) is analytic on a /eft-halfplane Reuw<+} with 
asymptotic behavior p!/?-*(r/a)—, a decay in the strip 
and on its left side, since rSa. However for the inte- 
grand of (10), in which r=a, no such halfplane of analy- 
ticity can’be found unless its form is changed, in such a 
way that the kernel in (10) becomes the modified 
Hankel function. This can be done by the decomposition 


8(u)Prje4u(0) = (7/2) [A(—p) — A(u)|/sin mu. (11) 
Substituting into (10), one obtains 
at oa a fu h (u) K,(yr)du.® (12) 
L 


where 
Whit) = (ry f r1Pu(yr, w/2)I,(yr)dr. (13) 
0 


From this definition of h(w) it can be determined that 
h(w) is free of poles on a right halfplane Reu> —4 with 
growth J,(ya)/p/? as | | —o in this halfplane; that 
h(w) has simple poles at w= —(2”+3), with unknown 
residues; and that the integrand in (12) is analytic on 
this left halfplane with decay yu /?(r/a)-* as | | —'0 
there (since r2a). Overlapping halfplanes are now de- 
termined for both integral equations, and (12) is also 
formally satisfied. 

It now remains to determine the unknown function 
g(u) or the related function h(u) by function-theoretical 
means. With this in mind it is convenient to write the 
integral equations in the form 


A 


a0 { Gur Gel oedds, 14) 
L 


~ 
IV 


a: 0 = f HHw)[(0/2Ka(7)/T) Ida, (15) 


where 


G(u) 
At (u) 


I 


w[g(u) P1/24u(0) — A(u) |(ya/2)-*/T(A — »), (16) 
T+ p)(ya/2)*A(u). (17) 


The purpose of introducing these factors is twofold: first, 
the expressions inside the brackets in (14) and (15) now 
have the asymptotic behavior (r/a)~* in the left and 
right halfplanes, respectively, and therefore approach 
at least a constant value (in the strip) or decay expo- 
nentially since r2a, resp. At the same time, these 
brackets are analytic in the respective halfplanes. 
Secondly, the functions G~(u) and H*(y) are not only 
analytic in left and right halfplane, respectively, but 
also possess algebraic decay in these halfplanes of analy- 
ticity. We shall call such functions + and — functions 


respectively. 


6 Note that this modification has introduced the Hankel function 
for the region 72a, and that this function explicitly satisfies the radi- 
ation condition. 


Leitner and Wells: Radiation by Disks and Conical Structures 


639 


Now, although the functions G~ and H+ are analytic 
in two different halfplanes, these two regions overlap in 
the strip containing L . Furthermore, by virtue of (11), 
(16), and (17), there is a relation between these func- 
tions which is valid at least in the common strip. We 
now apply the crucial step of the Wiener-Hopf tech- 
nique, namely, to exploit these relations for the con- 
struction of a function which is analytic on the entire 
“plane, while at the same time possessing algebraic de- 
cay—that is, a function which by Liouville’s theorem is 
identically zero. One finds 


2 


a) De) ite 
= Jaw - 5 (2) aah 
(= b f 


(18) 


The superscripts +, — are omitted from H(u) to avoid 
confusion. [H(—y) is a minus function.] 

The term on the left-hand side of (18) is a minus 
function. The term containing H(u) on the right-hand 
side is a + function. But the remaining two terms are 
mixed functions, 7.e., they possess simple poles on both 
halfplanes. Inspection shows that, while these two 
terms do have poles on the negative real axis, their 
asymptotic behavior as | | —o in other directions in 
the left halfplane is one of decay. Thus we may extract 
these poles by Mittag-Leffler series with the same resi- 
dues. The series themselves are plus functions. The dif- 
ference between the original functions and the appro- 
priate Mittag-Leffler series, however, are minus func- 
tions. Therefore these difference terms may be trans- 
posed to the left of (18), which side of (18) will be a 
minus function; and since this will equal the functions 
which remain on the right hand side, either side is 
identically zero by Liouville theorem. Thus one obtains 


cee ace) 


ioe} 


3 ) 
SS Ry [ (ut n+) _ Le Sn/(u + 1) 


0 


(19) 


ce 3 
SPE (n+ 20+ 5), 


0 


where the 7, are the known residues at u= —(2n+3/2) 
of the term in (18) involving A(w), vzz., 
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8 
Ty = (a3/2y/2ri)(—)"(ya/2)2"/T (2 + -). (20) 


The R, and S,, however, are the residues, at w= —1, —2, 
—3,---, and w=—(2n+3/2) of the term in (18) 
which involves H7(—y), 2.e., they are unknown. They 
have the form 


3 3 5) 
R, = | (20 +- =) c (» ac 5) / n\T? (2n a )| 
2 2 2 


-(ya/2)4"+3/2(2n + 3/2). (21) 
3 nN 1 n 
Sa = | nt ( o ) J ooer Ge ~)| 
4 2 4 2 
‘(ya/2)?"H(n). (22) 


item Loy we set yy equal to 1)°3/2,-2, 35/25 2. 
we obtain an infinite system for the numbers H(n), 
H(2n+3/2) whose coefficients are in the ascending 
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powers of ya and otherwise contain only factorials. If 
determined, then these numbers define ZZ for all other 
values of uw, through (19). 

We have approximated the solution of this system by 
terminated expansions of the unknowns in the ascending 
powers of ya and have obtained series expressions for 
the radiation patterns and scattering cross sections for 
the diffraction problem of the disk; these check com- 
pletely with the existing literature on the problem 
where other techniques of solution were used. We have 
also solved the problem of the vibrating.disk in a rigid 
baffle, and of transmission of sound waves by the circu- 
lar aperture.’ The former can be solved exactly by use 
of the Lebedev theorem. Extension of this method to 
the biconical antenna is in progress. 


7 A. Leitner and C. P. Wells, “Radiation by Disks and Conical 
Structures,” Int. Tech. Rep. No. 1. Contract No. DA-20-018-ORD- 
13354, Michigan State Univ.; 1955. 


Diffraction of Microwaves by Tandem Slits* 
LEROY R. ALLDREDGE+ 


Summary—tThe diffraction of a plane electromagnetic wave by 
two identical slits in tandem is investigated for normal incidence 
with the polarization parallel to the edges of the slits. 

Theory shows that the scattering cross-section coefficient is 
proportional to the imaginary part of the far field forward scattering 
factor. The stationary form of the scattering cross section is de- 
veloped in terms of the incident field and unknown currents on the 
edges of the conductors forming the slits. Calculations using the 
Kirchoff-type approximation in this stationary form for a tandem 
slit separation of 0.157) are in good agreement with the experimental 
values for slit widths greater than 0.5). 

Similar calculations for zero tandem slit separation, correspond- 
ing to a single slit, and for slit widths greater than 0.3, are in good 
agreement with those of the exact theory of Morse and Ruberstein,* 
and as determined experimentally. 

The infinitely long slits are approximated experimentally by use 
of a parallel plate system described earlier by Row.? The experi- 
mental results show an interesting resonance effect as the slit width 
changes. 


* Manuscript received by the PGAP, October 6, 1955; revised 
manuscript received, April 26, 1956. This paper is based on work done 
in partial fulfillment of the requirements for the Ph.D. in physics 
at the Univ. of Maryland, College Pk., Md. The experimental part 
was partially supported by ONR Contract N50ri-76, T.O.1. at Cruft 
Lab., Harvard Univ., Cambridge, Mass. The theoretical work was 
done at the Naval Ordnance Lab., White Oak, Md. 

t Operations Research Office, Johns Hopkins University, Chevy 
Chase, Md. 

1 P.M. Morse and P. J. Ruberstein, “The diffraction of waves by 
ribbons and by slits,” Phys. Rev., vol. 54, pp. 895-898; December, 
1938. 

2R. V. Row, “Microwave diffraction measurements in a parallel 
plate region,” J. Appl. Phys., vol. 24, pp. 1448-1452; December, 
1953. 


INTRODUCTION 


OST diffraction problems have been worked out 
| \ | using the classical Kirchoff theory which ex- 
plains quite well the diffraction of a wave 
through an aperture in a conducting screen when the 
dimensions of the aperture are very large in comparison 
to the wavelength. When microwave techniques are em- 
ployed it is not difficult to show that in the aperture 
near an edge of the conducting screen the electric field 
differs greatly from the incident field and on the shadow 
side of the screen the magnetic field is not zero near the 
edge, so the Kirchoff theory cannot be,correct. 

Most interesting diffraction problems which are 
amenable to exact analytical solution have been solved. 
The principal challenge in diffraction theory is the de- 
velopment of sufficiently accurate approximate methods 
for handling problems involving scatterers of complicated 
geometric shapes which have dimensions comparable to 
a wavelength. 

The variational procedure introduced by Levine and 
Schwinger* has been applied most successfully to planar 
diffraction problems where the scatterer consists of one 
or more apertures in a perfectly conducting plane of in- 
finite extent and where only the total scattering cross 

5H. Levine and J. Schwinger, “On the theory of electromagnetic 
wave diffraction by an aperture in an infinite plane conducting 


screen,” Comm. on Pure & Appl. Math., vol. 3, pp. 355-391; Decem- 
ber, 1950 
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section is found. This paper considers an extension of 
this theory to the case where the scatterer is no longer 
confined to a single plane. The particular problem 
treated here is that of a vertically-polarized plane elec- 
tromagnetic wave falling normally upon a vertical slit 
of infinite extent in a perfectly conducting plane with 
another identical slit in another conducting plane di- 
rectly behind the first and parallel to it. This compli- 
cated diffracting obstacle is simply referred to in this 
paper as tandem slits. The parameters in this investiga- 
tion are the slit width (same for both slits) and the 
tandem separation. 


THEORY 
Scattered Field 


Consider two parallel infinite planes I, and Ty at 
z=0 and zg=a as shown in cross section in Fig. 1. Plane 


gING 


POLARIZED TO PAPER 


an ee 
=3 
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Fig. 1—Cross-sectional view of tandem slit geometry. 


IT. contains the perfectly conducting infinitely thin 
screen Ig, and the slit Ia: and plane IT, contains the 
screen Ty and the slit Ij:. The parallel edges of the slits 
are at x=b/2, —b/2. The resulting problem is one in 
two dimensions, being independent of y, so that the 
planes become lines. The plus and minus superscripts 
indicate respectively the right and left sides of the 
planes hereafter referred to as lines in the corresponding 
two-dimensional problem. 

A plane electromagnetic wave polarized in the y 
direction is incident normally on the slit in the half 
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space z<0; it is desired to investigate the diffracted 
field. The incident wave is described by 


Eixe(7) = § exp (1ké-7) = 49 exp (kz) 


we & 
Hine(¥) = — — exp (ikz) (1) 
Zo 


where k = 27/) and ) is the wavelength, Z) = 376.6 ohms 
indicates the impedance of free space, £, 9, and 4 are 
unit vectors in x, y, and z directions, and a bar over a 
letter indicates a vector. The harmonic time depend- 
ence, exp (—iwt) is omitted throughout. Rationalized 
mks units are used. 

This incident electromagnetic field will establish cer- 
tain tangential components of E and F in the plane 
z=0 on the shadow side of the plane T,. There is a well 
known existence theorem that the field within a source 
free region is uniquely determined by either of these 
tangential components on the bounding surface of the 
region. If the plane T, did not exist and if we were inter- 
ested in the field only for z>0, we could think of our 
region as being bounded by TY, and a semicircle of in- 
finite radius in the space z>0. If this type of approach 
is attempted useful results can be obtained for the case 
of a single slit in terms of the transmission coefficient, 
but this method is no longer of use when the second 
plane with its slit is introduced. 

For this more complicated problem it is profitable to 
start with a consideration of the field scattered by the 
conducting screens forming the slits. Since we have 
postulated perfectly conducting screens, surface current 
densities K(7’) will be induced on the screens and will 
act as the source of the scattered fields. When the gen- 
eral scattering theory’ is reduced to two dimensions we 
may write for the y component of the scattered electric 
field 


Eyse(@) = — i Hom(e |e # |) year" 2) 


where Hy (k| 7 —7'| ), which will usually be indicated by 
Hy without repeating the argument unless considered 
necessary to avoid confusion, indicates the zeroth order 
Jankel function of the first kind, and the integration is 
over both sides of both conductors as indicated by the 
notation C=both sides of conductors I'go+Ty.=b:.s. 
(Tas +Tv2). Write the current as K,=K,S%+K,? and 
the scattered field as E,S°=E,S°°+E,?, where the 
superscripts SC) refer to currents and fields correspond- 
ing to the case where T.2 is complete (no slit). The 
scattered field for an infinite screen is 


WLL0 


= avg MoO B89) a" 


Ty 


{Sees = 


so we may write for the field scattered by the slit system 


4S, Silver, “Microwave Antenna Theory and Design,” McGraw- 
Hill Book Co., Inc., New York, N. Y., sec. 3.10; 1949. 


642 


Sewer Wo apa ns 
Ee) = a if Hy K,8¢o(7#’) aT’ 
Ran 


WiMo Bs 
— =f Hy) K,P(#) dP". (3) 
4 Jo 


The surface current density A,?(7’) is the difference 
between the surface current density which flows when 
the slits are present and that which would flow if there 
were no slits in the conductors (infinite conducting 
planes). For this reason K,?(#’) might be called the per- 
turbation surface current density caused by the slits. 
Looked at as a boundary value problem it can be said 
that the driving source for K,,?(#’) is localized in a finite 
area in the slit region, so that K,?(7’) tends to zero with 
increasing distance from the slits unless special condi- 
tions permit a portion of the energy to be carried away 
{from the slit region with no distance attenuation in the 
energy density. This latter condition may exist if the 
tandem separation is large enough to permit propagat- 
ing modes between the conducting planes. Such modes 
would represent energy scattered into the region be- 
tween the conducting planes. When these conditions 
exist the perturbation surface current density between 
the two planes can be writtten as 


Kaa) ms Git Canu Pate F Bera een 


where the subscript w.g. indicates that component of 
the surface current density associated with waveguide 
modes existing on the inner surfaces of the two conduct- 
ing plates. 

Surface current densities that support propagating 
modes in perfectly conducting waveguides do not give 
rise to electromagnetic fields outside the waveguides. 
Hence the [K,?(7’) ]y.¢. contribution to K,?(7’) in (3) 
can be omitted when calculating £,?(7) for 7 outside the 
region between the conducting plates or when calcu- 
lating the energy scattered by the slit system into the 
space outside of the region between the two conducting 
plates. The remaining components of K,?(#’) tend to 
zero with increasing distance from the slit system so 
that if in (3) it is understood that K,?(#’) does not in- 
clude the waveguide components the asymptotic ex- 
pression 


2 \12 
A)? = ( 5 ) exp (ikr) exp | —ikr’ cos (6 — 6’)} 
in kr 


To 


can be used to get an approximation for E,?(7) in the 
far field outside the region between the conducting 


plates. 
We will use only the special case of forward scattering 
so 6 can be set equal to zero finally resulting in 


1 1/2 
E,?*) = (=) exp (2kr)A(Z, 8), (4) 


18m kr 
where 


A(%, 2) = tp | f Ree ete 
h.s. Vay 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


October 


+ exp (—ika) KZ ar’ 


b.s.lag 


fi if | _Kyseogat’ |. (5) 


These expressions would permit the calculation of the 
far field forward scattered electric field if the perturbed 
currents were known over all the conducting screens 
forming the two slits. 


Scattering Cross-Section Coefficient 


The total power scattered by the conductors forming 
the two identical slits must be derived from currents in- 
duced on the conductors. Experimentally it is the for- 
ward scattered far field as given in (4) that can be 
readily measured in the laboratory. An interesting rela- 
tionship can be developed between this far field and a 
part of the total scattered power. 

As indicated earlier some of the scattered energy may 
propagate down the region between the conducting 
plates in waveguide fashion. This energy is in no way 
associated with the forward scattered far field and is 
therefore omitted from the definition of scattering cross- 
section coefficient discussed below. In what follows it is 
understood that the values of the scattered electric field 
ES°, the scattered magnetic field HS¢ and the perturba- 
tion surface current density K? do not include the com- 
ponents of those factors on [,2* and I'y27 associated with 
waves being propagated between the conducting 
screens. 

With the above limitations the time average of the 
power scattered into the space outside the region be- 
tween the conducting planes can be written as 


1 Cae S) 
Pae®’ = — Re | @-(ES¢ x Ase*yar (6) 
2 c 


where Re indicates the real part of the integral, 7 is the 
outward normal to the conducting surface over which 
the integral is taken and the * indicates the complex 
conjugate. 

In (6), substitute ES° = E8°4 EP and similar expres- 
sions for the magnetic field and for the time average of 
the scattered power, and note that 


on ess, ESCo oes Fine 
JES a ah 
Hsoo = Jai 
on |B ee ESseo = — Fine 
{sco — Fi 
H oe Lee 
and 
on Tyr and lepoee ESCo =— Fine 
and So = pee Fixe 


If we further recognize that, since the tangential com- 
ponent of the total electric field must be zero on the 
perfectly conducting surfaces of Tys and Tyo, the tan- 
gential component of E? must also be zero on these 
surfaces, and that 
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Piet: aes one eh 
2 Tao +Tay— 


+ i é. (Eine x Hine*) dP, 
Pope lane 


g. (Eine 4 Hine) dT 


we can obtain from (6) 


1 = “fi; 
Pie an Aone a Re Bi g. cme x Teva h 
Ly 


a al 


+ g. (Fine SG Hive*) dl 


| 
a 
| 
cts 


g- (Fine x H?*)dT 


g-(Bine x H?*)\dV 


+ 


=e { £- (Eine x H?*)dT 


-+ 


&- (Eine x Tear : (7) 


A few additional mathematical manipulations involv- 
ing the general boundary condition 


nx iH? = Kk. 
the identity is 
n (EX H) = —E-n~xX GE, 
the fact that 
eth Ko 
2 STR eS over T,— 


and the numerical value of Ei’ result in 


1 
Sie | { KP*(#)dP" 
2 b.s.lao 


K?*(r') dT" 


Pee = 


+ exp (tka) 


Das Vo 


Pe f _Rscorya" | (8) 


which in accordance with (5) becomes 


1 j Im A(é, 2) 
Payer? = —— Re [iA*(2, 2)] = ————__» (9) 
2a9 2wiLo 


where Im indicates the imaginary part of A(Z, 4). 

The scattering cross-section coefficient ¢ is defined as 
the average power P,,-” divided by the product of the 
slit width 6 and the incident energy per unit area k/2wpo 
and can be written as 


1 
= — Im A(é, 2). 10 
Saeee (2, 2) (10) 


This definition of ¢ might well be called the scattering 
cross-section coefficient of the slits since Pay-” is the dif- 
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ference between the power scattered by the conductors 
when the slits are present and that which would be 
scattered if there were no slits in the conductors (in- 
finite conducting planes). It must be remembered that 
this definition of o includes only the energy scattered by 
the slit system into the space outside of the region be- 
tween the two conducting plates. 


The Variational Formulation 


The scattering cross-section coefficient cannot be com- 
puted accurately using (5) and (10) because the cur- 
rents are not known, so approximation methods will be 
required to obtain quantitative answers to any specific 
problem. 

The simple Kirchoff approximation assumes that all 
of the perturbation currents K” on the conducting sur- 
faces in (5) are zero. This approximation yields a value 
of o = —2 which is equivalent to assuming that both & 
and H are unaltered in the first slit and that the second 
slit plays no part at all. This value approaches the cor- 
rect value for the case where the slit width is very large 
compared to the wavelength and where the tandem 
separation is small. 

It is quite obvious that such a simple approximation 
as this is incapable of predicting any resonant effects 
with changes in either the slit width or tandem slit 
separation. In what follows it will be seen that an ex- 
pression for o which is stationary with respect to varia- 
tions in K? will yield a Kirchoff-type approximation 
which does predict a resonant effect with changing slit 
width and changing tandem separation. 

Eq. (3) gives the perturbed electric field of the slit 
system in terms of currents. This expression was de- 
veloped for a plane wave traveling in the % direction. 
Consider the equation rewritten with additional sub- 
scripts 2 on the currents to indicate the direction of the 
incident field. The limits on the integral containing 
[ K,,8o(7’) |; can be shown as both sides of (Cai +Tw), if 
it is remembered that this factor is zero except over 
Tui. This form adds to the apparent symmetry. When 
the fact that £,?(7) is zero for 7 on Tae or Tye is used we 
obtain the following integral equation for [K,?(#) |é; 


ff Ho [xyseo7) ha” = fH (K.P) kar’ 
s iS 
(11) 


where the limit S indicates both sides of areas occupied 
by the slits Ta and Ty. 

So far we have always assumed the incident wave as 
traveling in the 2 direction in accordance with (1). Now 
consider the wave as coming from the other direction, 
that is, assume the incident wave to be 


[Eine(7) |_s = 9 exp { —ik(#-2)} 


for 7 on Tee or Tye 


SS Ve. 
Mb, 


0 


Venere) se (12) 


For this excitation, different currents will, of course, be 
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induced on the conducting screen. Multiply (11) by 
these new perturbation currents [K,?(7)]-s and inte- 
grate over both sides of ((ue2+Ts2) with respect to the 
variable (7) obtaining 


n= ffir 


Hy [Ky?(#) |\sd0’ al 


= fa fi [KyP(*) |-2Ho™ [K,S(7’) |edT, (13) 
which can be rewritten as 
= [ fix y(*) |-2Ho © [K,P (7) edt’ a 
= f dt’ hes [Ky?(#) |-sHo™ [K,S@o(7’) |dT 
S+C C 
Fy { if [K,?) 4H [Ky8(7’) Jedd. (14) 
3 


If the excitation is as given in (12) the integral equa- 
tion for [K,?(#) |_s corresponding to (11) is 


{ Ho [K,S0(7!) | dt" = i Hy [KP@)|-sdt", (15) 
vs G 


for 7 on Ty. or T'y2, where it must be remembered that 
[| K,,8o(7’) |_s is zero except on I',+. When this equation 
is multiplied by [K,?(7)]; and integrated over both 
sides of (Ta2t+tIs2) we obtain the following alternate 
form of (13); 


elles 


= (Pin a [K,P() Ho [K,800(7') sd. (16) 


Hy |K,?(#) |_sd¥d0 


When (15) is used (14) becomes 


I= [ fixrolen 


=f oar’ f (RPO Late [KF ear 
S+0 “vo 


H)©[K,?(#) |ed0’a0 


fi fi av’ { [KyS00(#") |eHo [K,80(7) ]_sdP. (177) 
(B S 


By adding and subtracting a common integral this 
equation can be rewritten as 


n= ff (era Lan 
_f ow finse 
_f af tne 


mah { [KPOH) Ho | K8eo(F) |ed0 al 18) 


[KP (#) |sd0'aT 
Ho [Ky8(7’) |zd0 


Hy [Ky8°o(#) |_sd¥ 
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AND PROPAGATION October 

Consider the scattered field for the case of no aperture 
in I, (complete screen) for energy incident in the 2 
direction. In accordance with (2), 


Wii 


Hy [K 807’) |sdV 
4 Jsic 


[Bo = 


= — exp (tkz) for z > 0, (19) 


where it must be remembered that [K,°°(7) |: is zero 
except on I~ 

Multiply both sides by [K,?(#) |_ and integrate over 
both sides of (Ta2+Tb2) to get 


ib exp (jk) [K,P@) Leah 


Se! an’ f [K,? (7) | 2H [K,8(r’) |e, 
S+C Cc . 
LOG Zee Os 


Multiply both sides of (19) by [K,S@(#) ] 
grate over both sides of (ui +I) to get 


(20) 


_s and inte- 


if exp (ikz)[K,S°o(#) |_2dT 
S 


ale av’ f Ik )8C0(7) |_ so [K,SC0(7") Jad D 


for’ > 0. ow) 


Employing (13), (16), (18), (20), and (21) the follow- 
ing expression can be obtained: 


JJ tere 
et [K,P(*) |-sHo? [KSo(#") | 3d 
wy fp a [rel Hy 


4 
yi ead fab cxP (ikz) [K,P(7) |-ar 
C 


Wo 


Hy [K,P(r') |d0'a0 


O[K,8e0(7') Led 


— ff exp (its) [xs aa | 


mh [ tx,seoey lett [K,yS(7) |_2dV'ar, (22) 
S 


If 2 and —% are interchanged this equation will be 
identical to its present form except for the term in the 
brackets on the right. Since the rest of the equation is 
invariant to this change it follows that the term in the 
bracket must also have the same value when 2 and —% 
are interchanged. When this exchange is made the 


bracketed term can be expressed in terms of A(é, 8) in 
accordance with (5). If we further define 
A.(é, 8) 
10 U9? 
= - Sf [le sow kno ix soum Lavan, (23) 
z S 
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then (22) can be written as 


14. ; 
[A (3 g) am A;,(2, ) | 


fe i) J [Ky?() ]-2Ho] Ky? (*) at’ ar 
c 
zs if aq” f [Ky?(7) | 2Ho [Ky8(7’) |sdP 
Ss Cc 


-forf [Ky?(7) |sHo [KyS(r’) | 2d. (24) 
Ss Cc 


Thus we have an expression which yields A (é, 8) and 
hence o if we know the perturbation currents that flow 
on the conducting screens forming the slits first for an 
excitation field traveling in the 2 direction and secondly 
for an excitation field traveling in the -—% direction. 

Eq. (5) which is much simpler also gave an expression 
for A(é, 2) in terms of the currents. The special value of 
(24) over (5) is that the former is stationary with respect 
to changes in K? provided [K”]; and [K?|_; are chosen 
so that (11) and (15) are satisfied. It follows that if 
6A (Z, 2) is made equal to zero then K? must satisfy the 
proper integral equation, that is, in principle, (24) pro- 
vides a method of determining K?. 


EXPERIMENTAL WoRK 


The microwave parallel-plate system built by Row? 
was available. With the plates closed a region is formed 
which is bounded on the top and bottom by plates of 
aluminum 4 feet by 8 feet. These plates are separated 
vertically by a distance less than one-half wavelength 
so that only the TEM mode can propagate if all source 
currents are in the vertical direction. 

The parallel-plate region was excited at one end by 
an open-ended waveguide. This gives, in effect, a line 
excitation and a corresponding cylindrical wave so that 
corrections must be made in the resulting data or the 
diffracting objects and the measurements must be re- 
stricted to a small region in which a plane wave ap- 
proximation is accurate. In this investigation only data 
for slit widths less than 1.5 \ were used to obtain scatter- 
ing cross-section coefficients, except for some phase in- 
formation for slit widths out of 2.0 \ as described later. 
The apparent line source was found to be 106 cm in 
front of the leading slit, which means that for a slit 
width of 1.5\ the amplitude was constant to within 
0.014 per cent across the slit and the phase was con- 
stant to within 3.4 degrees. 

The desired fields were measured by using a small 
vertical probe 101.3 cm in back of the slit nearest the 
source which was insulated from and extended up 
through the bottom plate. The signal from this probe 
excited a waveguide and this signal was then used in 
various standard ways for the determination of phase 
and amplitude of the original signal. 

Construction of the tandem slits used in the parallel- 
plate region presented several difficult problems; the 
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diffracting edges forming the slits should be very thin 
compared to a wavelength, should make good contact 
with the top and bottom plates, and should be quickly 
and accurately adjustable in tandem separation and slit 
width. The requirements appear to be mutually ex- 
clusive, and it is true that compromises must be made 
in all of them. The final design used 0.003-inch thick 
silver foil for the leading edge of the slits but used a 
thicker construction back from the edges to facilitate 
the mechanical adjustment of slit width and tandem 
separation. 

Fig. 2 (next page) is a photograph of one of the dif- 
fracting edges used to form the slits. Fig. 3 shows the tan- 
dem slits in position on the bottom plate of the parallel- 
plate system. In this photograph, the top plate is in 
the open position. In this position major adjustments 
can be made in the diffracting region. 

Five watts of microwave power at a free space wave- 
length of 3.182 cm from an X-21 two-cavity klystron 
was modulated at 1000 cps by a ferrite Faraday rotation 
switch. Amplitude measurements of the field were 
made by feeding the rectified 1000 cps signal into the 
input of a tuned audio-amplifier. The output of this 
amplifier was connected to a Ballantine vacuum-tube 
voltmeter. In making phase measurements, the phase of 
the signal picked up by the probe was compared with 
the phase of the known, but adjustable, reference signal 
selected from a terminated slotted guide. The location 
of the probe in the slotted guide was measured to 0.001 
inch by an Ames gauge. 

Experimentally it would be very difficult to make the 
absolute phase measurements called for by the theory. 
Most of the difficulty regarding the phase reference for 
the far field factor can, however, be dispelled by observ- 
ing for a single slit (@=0) both the way in which the 
phase changes as the slit separation is increased and 
how the resulting scattering cross-section coefficients, 
computed by assuming various zero phase reference 
points, compare with the exact theoretical values.! Ex- 
perimentally it was found that the best fit occurred 
when the far-field factor was taken as imaginary for the 
most extreme slit width used (2A). This procedure was 
adopted for all the calculations made. 

The above procedure for choosing the zero phase 
angle is justified by the agreement with the theoretical 
curve for a single slit. When the experimental and 
theoretical curves for the single slit are compared for 
large slit widths, however, discrepancies begin to appear 
for slit widths greater than 1.3. These discrepancies 
are not of the type that can be corrected by changing 
the zero-phase angle-reference-point, without causing 
major discrepancies for small slit widths. For this reason 
it seems proper to conclude that at these large slit 
widths, the undesirable effects of a cylindrical incident 
wave and of the probe being too close to the slits are 
becoming important in determining the amplitude of 
the signal. Because of this, the final results given here 
do not show data for slit widths greater than 1.4A. 
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Fig. 2—One of diffracting edges used to form slits. 


Z 


Fig. 3—Tandem slits in parallel-plate region with plates in open position. 


The final results are displayed in Fig. 4, where both 
the amplitude of the received signal divided by slit 
width and the scattering cross-section coefficients are 
plotted as functions of the slit width for a variety of 
tandem slit separations. All of the experimental data 
has been normalized to fit the Morse-Rubenstein theo- 
retical values for the single slit at b/A=1. Actually the 


absolute value of the scattering cross section divided by 
two is plotted so as to be comparable with existing 
theoretical data for the transmission coefficient, 
t= —a/2, for the single slit. 

The experimental result for the magnitude of one-half 
the scattering cross section (or the transmission) coef- 
ficient for zero tandem separation (single slit) are seen 
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Fig. 4—Scattering cross-section coefficient « and normalized 
amplitude in far zone for tandem slit diffraction. 


to agree quite well with the Morse-Rubenstein com- 
puted values. As the tandem separation increases, the 
primary peak in the scattered cross-section coefficient 
decreases and moves toward larger slit widths. The 
curve changes quite markedly when the tandem separa- 
tion goes from 1.5 cm to 2.0 cm. This change may be 
associated with the fact that at 2.0 cm separation the 
space between the planes forming the slits can now act 
as a waveguide for the Hip mode. Similar changes seem 
to occur at each new tandem separation that can permit 
the next higher mode to propagate between the con- 
ductors forming the slits. This change in behavior, 
which seems to be associated with the modes that can 
be propagated between the conductors, appears in 
greater relief in the amplitude data. Large peaks in the 
amplitude have developed for tandem separations of 
1.5 cm, 3.0 cm, and 4.5 cm. These peaks are greatly re- 
duced when the slit separations increase respectively to 
2.0 cm, 3.5 cm, and 5.0 cm. It would be interesting to 
investigate this effect in detail to determine how sharp 
the transition is. It might be caused by a large amount of 
energy being coupled into the waveguide formed by the 
slit assembly or it may be that the incident wave sees a 
greatly modified input impedance at the leading slit. 


COMPARISON OF THEORY AND EXPERIMENT 


In general when trial functions are used in a varia- 
tional formula there is no direct way to determine the 
limits of error that will result. Therefore, it seems wise to 
first use an approximate current distribution for the case 
of a single slit, where an exact answer is available and 
hence the error can be determined, before time is spent 
trying to work out the more complicated tandem slit 
problem with similar types of current distribution. 

To do this we start with the variational form of 
A(é, %) given in (24). When this equation is reduced to 
the case of a single slit, we get 


A(é, 8) = — 100” f { [K,?(#) |sHo [Ky?(7’) |edVd0 
sets 


ay 


€ Wie 
ah. 212 W9? (=) f arf [K,P(7) |:Ho dr 
Ho Vo, Pa 


= iotoe f [aooarar. (25) 
Tai 


The simplest approximation is to set the current 
density equal to zero. This is a Kirchoff-type approxi- 
mation which becomes 


A(&, 2) = Ax(4, 8) = — iotwoee ff Hydi'dl. (26) 
Tie 
Referring to Fig. 1, we can write 
b/2 
Im A,(, 2) = — wae f foc 7 — #|)dr'dr. (27) 
—b/2 


The Bessel function Ti(R| 7—7'|) isan even function so 
the absolute sign indication on |*—#’| can be removed 
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since the vectors * and # are collinear. By using a 
change in variable and the integral formula?® 


it Than Bo D Fossa) 


we get the double infinite sum 


Im A,(é, 2) = — 8 py J 9y4-242( kb). 


v,p=0 


It can be shown that for a single slit the transmission 
coefficient ¢ is equal to minus one-half the scattering 
cross-section coefficient, so we have 


oh 4 
9 J 2v4-2p+2( RO). (28) 


SEEN seat 


This double infinite sum is amenable to numerical cal- 
culations for relatively small values of b/d. This nu- 
merical calculation is shown in Fig. 5, where it can be 
readily compared with the exact calculations of Morse 
and Rubenstein and the experimental value. 


O OO EXPERIMENTAL 
EXACT THEORY (MORSE AND RUBENSTEIN) 


———=— APPROX THEORY (KIRCHHOFF-TYPE®* 
2 


Jev+ 2p+ 2") 
=0 


Fig. 5—Scattering cross-section coefficient o for single slit. 
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culation are very striking. The currents on the edges of 
the slit, which have been assumed zero, should play a 
lesser part for larger values of b/ so the approximation 
should be better as b/\ increases to values beyond those 
shown. The currents on the edges should be relatively 
more important as b/d is made very small so it is not at 
all surprising that the Kirchoff-type approximation is 
less accurate for small slit widths. 

The agreement between the Kirchoff approximation 
and the exact theory for the single slit is good enough to 
encourage one to attempt the more tedious calculations 
required for a similar approximation for the case of 
tandem slits. 

Setting K? equal to zero in (24) and employing (23) 
we get for the Kirchoff-type approximation for tandem 
slits 


1 U9” 
A z(8, 8) at 


f { [KS©o(#) |sHo [K yS©o(#’) |_2d TaD. 
S 


When we apply the conditions 
1/2 
[K,So], = — 2[H,jne]; = 2(= *) on Tai 
Mo 


= (0 on I,;t and on both sides of Ty; 


and 


ANE 
ee vies ae PAP Re on 2(*) exp (—ika) on Pei 
Ko 


= 0 on Ty; and on both sides of Tay, 


the approximation becomes 


Ax(é, 2) = — iB? 


-exp (—ika) arf HoMd¥", (29) 
one side of Tg; 


one side of Ty, 


where we have used the equality w2u€) = R?. 

By employing a change in variable (see Fig. 1), inte- 
grating by parts to remove a singularity, and rearrang- 
ing the terms we can finally get 


2 1/2 b 
A,(2,%) = — 12k exp (ito) fo ae El (= — rs) + a eS — rs) 


+m [f ee) ea} (Er) 


[(b/2—r,) 24a] 42 
+2f Hy0(y) [y* = (ba)? }8dy + 
ka 


The Kirchoff-type approximation is seen to give 
reasonable agreement with the exact calculations and 
with the experimental results for b/A greater than 0.4. 
The resonant peaks in o even for the approximate cal- 


5 E, Jahnke and F. Emde, “Tables of Functions with Formulae 
and Curves,” Dover Publications, New York, N. Y.; 4th ed., p. 145, 
1945. 


kE(b/2+r,) +a] M2 


LN iil 


(ha)*)ay . (30) 


k[(b/2—n,)-+02}” 


In this form although the labor involved is very great 
it is amenable to numerical evaluation using a desk 
calculator. This calculation has been made for a slit 
separation of a=0.5 cm (ka=0.986) and for a range of 
b/X from 0 to 1.4 to correspond with one of the experi- 
mental cases shown in Fig. 4. 

A comparison of this theoretical Kirchhoff-type ap- 
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proximation with the experimental values is shown in 
Fig. 6. The agreement is fairly good; it still predicts 
the resonance effects quite well. The presence of the 
term exp (—7ka) in (30) makes it possible for even this 
simple approximation to display at least in a qualitative 
way the resonance effects with changes in slit separation 
as were noted experimentally. 


OOO EXPERIMENTAL 
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SLIT SEPARATION = 0.5 CM 
= 3.185 CM 


1.0 1.2 1.4 
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Fig. 6—Scattering cross-section coefficient o for tandem slits. 


It would be desirable to determine the currents on the 
screen by utilizing the stationary properties of (24), but 
the work appears very formidable and this is left for 
future research. 

Since the Kirchhoff-type approximation in which all 
currents have been assumed zero has been shown to be 
a fairly good approximation it would seem plausible that 
almost any simple trial function which had the general 
characteristics of the actual current should yield an im- 
provement. With this thought in mind several attempts 
were made at using simplified trial functions for the 
current. The fact that each one tried failed to improve 
on the original Kirchhoff-type approximation leads to 
an interesting study in itself. 

Moullin and Phillips® have discussed the current in- 
duced in a conducting ribbon by the incidence of a plane 
electromagnetic wave. They conclude that the current 


6 E. B. Moullin and F. M. Phillips, “On the current induced in a 
conducting ribbon by the incidence of a plane electromagnetic wave,” 
Proc. IEE, vol. 99, pp. 137-150; July, 1952. 
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on one edge is independent of the current on the oppo- 
site edge provided the second parallel edge is more 
distant than \/2. 

They assumed that for purposes of computing diffrac- 
tion patterns they could assume a uniform current 
density extending back from the edge in such a way that 
its amplitude width integral was equal to the true 
amplitude width integral. 

This type of approximation was tried in (25) for the 
single slit. The resulting values of the scattering cross- 
section coefficient were not as close to the Morse- 
Rubenstein values as was the original Kirchhoff-type 
approximation. Furthermore the value seemed to de- 
pend quite critically on the width of the assumed uni- 
form current density where the amplitude was adjusted 
each time to preserve the proper amplitude width 
product. 

The stationary characteristics of A(é, #) as given in 
(24) merely mean that the variations in A(%, 4) are in- 
dependent of the variations in the current density pro- 
vided K” obeys the proper integral equation. In general 
we may write for the variation in A(&, 8), 


5A(8, 2) = ab KP + B(6K?)2 + y(6KP)34+.---. 


The stationary property of (24) means that the coef- 
ficient @ is zero, but in general 8, y and the other ex- 
pansion coefficients are not zero. If a AK? is chosen as. 
a trial function such that 6K? is no longer very small 
then the higher order terms in the above expansion may 
make large contributions to 6A(é, ), thereby resulting 
in poor approximations to A(é, 2). The nature of the 
complex integrals in (24) involving the imaginary part 
of H)™ which goes to infinity for zero argument seems 
at least to make it reasonable that an assumed uniform 
current back from the edge may give a worse approxi- 
mation than when the current is assumed to be zero. 
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Transmission Characteristics of Inclined Wire Gratings* 
O. J. SNOWt 


Summary—Small diameter parallel wires were imbedded in 
thin plastic sheets and located closely before an antenna dish re- 
ceiving plane wave X band energy. Polarization was parallel to the 
wires, and the grating interval was varied between a fifth-wavelength 
and a whole wavelength for different panels. Received intensity 
was measured for varied grating tilt angles about an axis lying parallel 
to the wires and near the center of each panel. Sharp and intense 
transmission dips were observed for tilt angles at which the para- 
sitic reradiation maxima lay in or near the end-fire direction. The 
shapes and angular positions of the transmission vs tilt angle curves 
are approximated by a tentative theory which assumes that the in- 
put impedance of the grating is independent of tilt angle and that the 
apparently absorbed power is proportional to the areas under plots 
of antenna array patterns. A more precise theory which includes 
the effect of varying input impedance was required to predict ap- 
proximate amplitudes as well as sharp transmission dips of smaller 
magnitude. 


EXPERIMENTAL PROCEDURE 


HE principal features of the experimental ar- 
rangement are illustrated by Fig. 1. The X band 
transmitting antenna was placed at a sufficient 
distance from the panel to insure a nearly plane incident 
beam over the grating area, and the receiving dish an- 
tenna was placed a few inches behind the panel. The 
latter was continuously rotatable about an axis lying 


FIBREGLAS LAMINATE PANEL 


INCIDENT 


E-VECTOR 


\ PARASITIC 
RADIATION \ 


Fig. 1—Experimental arrangement. 


perpendicular to the plane of the figure. With this con- 
figuration, the received intensity was measured vs a 
varying panel tilt angle 0 for several different values of 
wire separations d, the incident electric vector being 
maintained parallel to the wires. The effect of interac- 
tion between dish and panel on the results was reduced 
by using as the received intensity the geometric mean of 
the maximum and minimum intensities as the spacing 


* Manuscript received by the PGAP, September 2, 1955; revised 
manuscript received, April 9, 1956. 

t U.S. Naval Air Dev. Ctr., Aeronautical Elec. and Elect. Lab., 
Johnsville, Pa. 


of the panel center relative to the antenna dish was 
varied over approximately a half-wavelength. The panel 
widths as shown in the plane of Fig. 1 were about 4 feet 
and both 10-inch and 18-inch receiving antenna dishes 
were used. In all cases a 70 per cent copper—30 per cent 
zinc brass wire material was used, with sizes in the 
number 26- to 29-gauge range, although neither small 
variations in conductivity nor in diameter appeared to 
affect the results. Except for the case where the wire 
separation was }-inch, where the wires alone were used, 
they were all imbedded near the centers of laminated 
Fiberglas panels (six piles of number 180 Fiberglas cloth 
plus Laminac P-43 resin). Since the fabrication details 
varied somewhat from panel-to-panel, the panel thick- 
nesses varied from 0.044-inch to 0.063-inch as indicated 
in the upper right hand corner of Fig. 1. The wire 
spacings chosen were 4-, #-, 1, and 14 inch. It was ex- 
pected that more or less uniform absorption would occur 
for wire spacings less than about one-half wavelength, 
therefore, except for the single panel with quarter-inch 
spacings, spacings greater than one-half wavelength 
were chosen. 


EXPERIMENTAL RESULTS 


The principal experimental results obtained with such 
an arrangement are given in Fig. 2. Here per cent power 
transmission is plotted vs grating tilt angle for four dif- 
ferent panels, each with different wire spacings. For 
simplicity in presentation, only regions of the curves 
where the variation in transmission is marked are in- 
cluded. The curves are, of course, in all cases approxi- 
mately symmetrical about the zero tilt position. It is 
noted that, except for the grating with {-inch wire 
spacings, sharp and pronounced transmission dips occur 
at tilt angles which increase in value as the spacings are 
decreased. 

The separate points which are present in Fig. 2 indi- 
cate measured transmissions of wire-free panels of ap- 
proximately the same construction as those in which the 
gratings are imbedded. It is observed that these panels 
exhibit no absorption maxima of their own; also, in 
some instances, the presence of the wires appears to 
slightly improve the over-all transmission. Thus the 
transmission minima may be attributed, primarily, to 
the effect of the wires alone. 

For the case of 1.25-inch wire separation, that is, ap- 
proximately one wavelength, it may be observed in Fig. 
2 that a single major dip occurs at normal incidence. As 
the wire separation decreases from one wavelength, this 
single minimum appears to split into two minima as is 
evidenced by the curve for a 1-inch separation. Here, 
there are major dips in the neighborhood of plus and 
minus 15° along with a number of minor variations. 
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Fig. 2—Experimental transmission curves for 
four wire gratings. 


When the wire separation is reduced to 3-inch, the 
principal minima have moved to approximately plus 
and minus 44°. Referring to Fig. 2, it may be seen that 
these major dips occur at tilt angles which are given, 
very nearly, by the formula: 


nN 
sin7! (— _ 1). 


It is immediately apparent that this is also the condition 
at which the array would be expected to reradiate para- 
sitically in a reverse end-fire direction. In other words, 
the apparently absorbed energy is reradiated in an end- 
fire direction when the tilt angle reaches the proper 
value for this to occur. This relationship is graphically 
illustrated in Fig. 3. Slight errors in angular measure- 
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Fig. 3—Observed condition for minimum transmission. 


ment, panel flatness, and wire spacings are present such 
that the measured transmission minima are no more 
accurate than about one-half of one degree. The extreme 
differences between the measured values and values 
computed by the formula above are only slightly higher 
than this. 


THEORY 


First, a simple theory which predicts the angular 
positions and approximate shapes of the transmission 
curves will be briefly discussed, after which the more 
complex and accurate theory will be outlined. The 
simpler tentative theory assumes that the total im- 
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pedance of each wire is invariant with tilt angle and 
that the current peak amplitude is identical and uni- 
form over all of them. With these assumptions, it can be 
stated that the total reradiated power will be propor- 
tional to the integral of the power antenna array pat- 
terns, where the relative phases between the currents of 
successive wires is given by d sin @ (see Fig. 3). Approxi- 
mate patterns computed in the plane of the gratings for 
the case of 1-inch wire spacings (47 wires total), and for 
three significantly different incident beam directions are 
plotted in Fig. 4. Note that only reradiated energy is 
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Fig. 4—Approximate parasitic radiation patterns for varying 
tilt angle (1-inch wire spacing). 


considered here, and that the large effect of the main 
transmitted beam is not shown. It is observed that 
radiation peaks are found in the principal transmission 
or incident beam directions noted by lobes “a” in the 
figure in addition to lobes noted by “b” which are peaked 
in the specular reflection directions. The absolute in- 
tensities of these lobes, as well as their widths, change 
comparatively slowly with changing incident direction, 
still assuming that the wire current peak amplitudes do 
not change with incident direction. In pattern “A” the 
incident or tilt angle is appreciably less than 15°, a con- 
dition for which the measured transmission is well above 
90 per cent, with only minor transmission dips occurring. 
From the antenna array formula it can be easily shown 
that, with this condition, the only significant lobe other 
than those in the main transmission and specular reflec- 
tion directions is lobe “c” which points in an end-fire 
direction. The amplitude of this lobe grows but remains 
relatively small until the critical angle of incident beam 
direction, in this case 15°, is approached. Pattern “B” 
represents the case where the incident beam direction 
is exactly this critical value. Here the end-fire lobe has 
grown rapidly to its maximum value, and while its 
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amplitude is identical to that of the other two lobes, as 
given by the antenna array formula, its width is con- 
siderably greater, thereby representing a decidedly 
greater quantity of power. As the incident beam direc- 
tion angle exceeds the critical angle, the lobe “c” splits 
into two separate lobes, represented by “c,” and “c,” 
in pattern “C.” These two lobes maintain the same 
amplitude as the others, however, their widths are suf- 
ficiently small so that the energy which they represent 
is significantly less than the energy represented by lobe 
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THEORETICAL RESULTS 


Using the simplified assumptions stated above, areas 
under computed antenna patterns were determined for 
the case of one-inch wire separation and plotted as a 
function of grating tilt angle. The detailed results of 
these calculations are omitted here, but it may be stated 
that the shape and angular position of major transmis- 
sion dips were predicted reasonably closely by the com- 
puted curve. Of course, nothing could be computed by 
this simplified method which would give information as 
to the absolute magnitudes of the transmissions. The 
minor maxima and minima which lie between the major 
dips were not predicted at all by the simple theory. 

Both in order to investigate the significance of the 
minor variations and to obtain a theoretical value for 
the absolute magnitudes of the transmissions, a more 
thorough examination of this problem was then made. 
It was assumed that if the absolute magnitudes of the 
wire currents along with their radiation resistances 
could be obtained, then the parasitically radiated power 
could be computed, although omitting the effect of inter- 
ference between the main transmitted beam and the 
parasitic lobe radiation which is pointed in the incident 
beam direction. As a basic reference for this more 
thorough study, the article by Groves,! was used.” Al- 
though Groves was concerned with the case of pairs of 
grids with normal incident radiation, he obtains the im- 
pedance per meter of the central wire of a single grid as 
one step in his derivations, still at normal incidence. 
Following this derivation, it is a relatively simple matter 
to modify the original equations to make them valid for 
an arbitrary grid tilt angle and for a finite, rather than 
an infinite number of wires. 

The equation thus obtained was used for the basic 
computations of the present study. Because of its ortho- 
dox nature it is not given here. The formula is composed 
of two terms which give the self impedance of a single 
wire, plus terms giving the mutual impedance of the 


1 W. E. Groves, “Transmission of electromagnetic waves through 
pairs of parallel wire grids,” J. Appl. Phys., vol. 24, pp. 845-854; July, 
1953. 


* See also: J. Shmoys, “Diffraction of Electromagnetic Waves by 
a Plane Wire Grating,” Res. Rep. EM-18, Mathematical Research 
Group, New York University, New York, N. Y.; March, 1950. 
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central wire due to each pair of wires located equidistant 
from it. For the case of the panel with 1-inch wire sepa- 
rations with 47 wires, a total of 25 terms are obtained. 
This formula also assumes that the wires are of infinite 
length and that all of them carry the same current. It 
was known that these assumptions did not hold rigor- 
ously for the case under study here, and thus prelimi- 
nary studies were made to determine approximately the 
degree to which they did or did not hold. The first as- 
sumption was determined to be sufficiently valid, since 
it was found that changing receiving dish size made no 
appreciable difference in measured results, indicating a 
reasonable current uniformity lengthwise over the wires. 
The second assumption is more critical, since the mutual 
impedances of the wires near the edges of the grating 
are obviously quite different from those near the center 
of the grating. However, the receiving dish sizes were, 
in practically all cases, of diameters of no more than 
half of the grating projected widths, indicating that the 
currents in the wires near to the dish should be rela- 
tively uniform. Thus it was concluded that this assump- 
tion was probably not a serious one either. The im- 
pedances for each of the wires could, of course, be com- 
puted by taking the proper summations for each one 
separately, but this presented a tedious task and it was 
avoided by accepting this second assumption. Even 
with the simplifying assumptions, the computations 
involved are considerable. 

Further preliminary investigations were made to de- 
termine if computational simplifications could be made 
by utilizing infinite summation formulas similar to those 
utilized in Groves’ article. To gain information on this 
subject as well as general knowledge of the phenomenon 
involved, transmission curves were measured using the 
4-foot wide panel with 1-inch wire spacings and with 
1 foot of each edge being covered on its front surface 
with a microwave absorbing material. A 10-inch dish 
size was used so that more than 6 inches of clear grating 
beyond the dish edges was maintained. From the results 
plotted in Fig. 5, with and without the absorber, it ap- 
pears that the wires beyond the dish edges have an ap- 
preciable effect on the measured transmission curve. A 
similar result was obtained by moving the panel 
laterally so that its edge which pointed towards the dish 
was near the dish. Moving the panel laterally in the 
opposite direction near its extreme caused an apprecia- 
ble sharpening and deepening of the absorption dip. 
These results indicate that the grating lengths lying 
beyond the dish edge in the direction of the end-fire 
parasitic pattern have no large effect on the transmis- 
sion curve, while the grating lengths lying beyond the 
edge in the opposite direction have a significant effect. 
It implies that the degree of coupling between the wires 
of one edge to the wires in the central region is appre- 
ciable. This also indicates that computations must be 
based on the impedance formula with the finite number 
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Fig. 5—Measured transmission with and without 1-foot 
absorber over tilted panel edges. 


of terms and that use of infinite convergence formulas 
are not valid in this application. 

The microwave power which is parasitically radiated 
per meter of wire per unit of incident field intensity 
squared is given by the ratio of the resistive component 
of the impedance per meter to the square of the magni- 
tude of the total impedance per meter. In the computa- 
tions for the problem being considered the fraction of 
energy incident upon the grating, which is parasitically 
radiated, is assumed to be given by the ratio of the 
power radiated per meter for the central wire alone to 
the power per meter per grating interval incident upon 
the grating. When this radiated fraction is subtracted 
from unity, the fraction of nonradiated power is ob- 
tained. The latter computed value together with the 
measured transmission percentage is plotted in Fig. 6 
as a function of tilt angle for the case of 1-inch wire 
spacing. Note that both major and minor transmission 
dips given by the measured data are predicted by the 
computed results. The amplitudes agree as to order of 
magnitude, although the discrepancies in this regard are 
not at all negligible. 

The sign and order of magnitude of the discrepancy 
at the transmission minimum which occurs at a tilt 
angle of 15° can be explained as due to interference be- 
tween the transmission beam unaffected by the wires 
and the parasitic lobe which points in the incident beam 
direction. That the phasing is correct for this can be 
seen by an observation of Fig. 7 which gives a plot vs 
tilt angle of the phase delay of the parasitically induced 
currents relative to the incident field intensity. Note 
here that at the angle of 15° where maximum absorption 
occurs, the phase drops suddenly from a value which is 
near quadrature to a zero or “in-phase” value. If then, 
the equivalent current sheet method is applied to the 
incident field intensity between the wires, it will be 
found that this equivalent current has a 180° or “out- 
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Fig. 6—Measured transmission compared with 
computed nonradiated energy. 
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Fig. 7—Resistive component of impedance and phase 
angle of induced current. 


of-phase” value at a 15° tilt angle, which will cause de- 
structive interference with the fields resulting from the 
wire currents. 

The amplitudes of the minor transmission dips are 
computed to be greater than the measured values. The 
phasing between wire fields and incident fields are not 
destructive as they are in the case of the major dip; 
however, it appears that they may be slightly con- 
structive which would explain a part of the discrepancy. 
The remaining part of the differences are probably due. 
to the effect of the assumptions described above, re- 
sidual experimental errors, and to the presence of the 
thin Fiberglas laminated sheet which was neglected. 

Fig. 7 also gives a plot of the radiation resistance per 
meter of the central wire as a function of tilt angle. Note 
that the wiggles in this curve are approximately oppo- 
site in phase to those in the transmission and phase 
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angle curves, as might be expected. The peak of radia- 
tion resistance, however, occurs at a tilt angle which is 
about one degree higher than that for maximum ab- 
sorption. The total impedance curve, not plotted here, 
has a shape approximating that of the resistance curve, 
which indicates that the assumption made previously, 
that impedances and currents are invariant with tilt 
angle, are completely invalid. The degree of agreement 
between theory and experiment obtained when using 
that assumption may be fortuitous. 
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CONCLUSION 


In conclusion, the transmission minima introduced by 
thin wire gratings imbedded in thin plastic sheets with 
wire spacings lying between a half and one wavelength, 
when placed close to a receiving antenna dish whose 
dimensions are smaller than the grating size, have been 
predicted theoretically with fair results. Sharp trans- 
mission dips, accentuated by interference effects, occur 
at tilt angles where parasitic reradiation is in or near a 
reverse end-fire direction. \ 
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Summary—tThe diffraction by a grating is examined (for spacing 
large compared to wavelength and dimension of grating element) 
for wavelengths in the neighborhood of the ‘‘Rayleigh” wavelengths. 
The shape of the elements, and their size in wavelengths is un- 
restricted. The results, including the effect of interaction, are ex- 
pressed in terms of quantities relating to single scattering. Some 
properties of certain determinants formed from single scattered 
amplitudes are derived. The results are compared with those obtained 
by other authors, using various restrictions on the parameters. 


I. INTRODUCTION 


HE PURPOSE of the present work is to present 
Oe: results of an investigation of the behavior of 

diffraction gratings of cylinders for wavelengths 
neighboring certain critical wavelengths, the so-called 
“Rayleigh Wavelengths.” The analysis is based on the 
results of a previous investigation,!” in which expres- 
sions were derived for the Fourier amplitude of the dif- 
fracted field of a grating of arbitrary cylinders, subject 
to the restrictions that the spacing is large as compared 
to cylinder dimension and wavelength, and that the 
Rayleigh wavelengths are excluded. The expressions 
there obtained for the Fourier amplitudes only involved 
quantities associated with a single scatterer, 7m spite of 
the fact that interactions were taken into account. The 
principal result of the present paper is the extension of 
the conclusions of Karp? to the case of critical wave- 
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lengths, in which the significance of interaction is maxi- 
mal. The required results of Karp! and Karp? are re- 
capitulated in Section II of the present work. They form 
a special case of a general theory® for the calculation of 
the field scattered by a combination of obstacles in 
terms of the scattering patterns of the obstacles in isola- 
tion. Experience with the case of diffraction by a pair of 
half planes (7.e., a slit)* indicates that the accuracy of 
this kind of calculation is good down to spacings of one 
wavelength. 

In Section III of the present paper we discuss the be- 
havior of the Fourier amplitude of the grating in the 
limit in which the critical or Rayleigh wavelengths are 
approached. In this limit the “single scattering” result 
diverges, while the result given in Section II is inde- 
terminate. However, the limit is evaluated explicitly in 
the present work. We also include in this section a dis- 
cussion of certain determinants relating the scattering 
amplitudes of a single cylinder; the behavior of these 
determinants is crucial for the limiting process in ques- 
tion. 

Section IV is concerned with a brief discussion of 
energy conservation. It is shown that the “single 
scattering” result conserves energy with an error which, 
while small for general large spacings, contains a factor 
which diverges for critical spacings. On the other hand 
the results deduced from Section II are free of this 
defect. 

In Section V we compare our results with those of 
previous investigators. Most previous work has been 
based on other types of restrictions than our own. How- 
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ever, Twersky,*:® who treated the case of circular cylin- 
ders for noncritical spacings obtained equivalent results 
for the contributions of the individual cylinders in that 
case, for certain angles of incidence and these results 
were extended by Magnus (see Appendix)*. to arbitrary 
angles of incidence, subsequent to the developments by 
Karp:! 

A more interesting comparison is made possible by 
considering the application of our results to the reflec- 
tion grating, which is the seat of the celebrated “Wood's 
Anomalies” discussed by Wood,’ Rayleigh,® Fano,® and 
Strong.!® A definitive theory has been given by Art- 
mann." Since the only restriction of Artmann’s analysis 
is that the corrugations be small compared to wave- 
length, a comparison with our theory is possible in the 
case in which corrugations are small and the spacing is 
large. We specialize our results to small corrugations in 
Section 5, and, after giving explicit approximate for- 
mulas for the scattering amplitudes of a thin strip in the 
case of when the magnetic vector is parallel to the 
generators of the grating, we are able to show agreement 
with the results of Artmann. 

On the other hand it is pointed out that both for the 
resonance and nonresonance cases our results are in dis- 
agreement with results given in a new report by 
Twersky.” The nature of the discrepancy is explained, 
and possibilities for its removal are suggested. 


II. PRELIMINARY RESULTS 


It is well known that the far field form of every out- 
going wave function! takes the form ~- 


ay T 
H (kr) f(@) = \/— exp (i — i7) 10 (1) 


where 7 and @ are polar coordinates. Also, in a problem 
dealing with diffraction by a combination of obstacles 
Green’s theorem shows that we can regard the scat- 
tered field as composed of a sum of outgoing wave func- 
tions, each regular in the exterior of one of the scattering 
bodies. Suppose the bodies in question are far apart, so 
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that all dimensions of each body, and the wavelength as 
well, are small compared to the distance to any other 
body in the array. Then it can be shown from (1) that 
in the neighborhood of any body in the configuration the 
total field is approximately equal to the incident field, 
plus the field scattered by the body in question, plus a 
sum of plane waves arriving from the directions of the 
other bodies. (The complex amplitudes of the latter 
plane waves are at first unknown.) Consequently the 
actual field contributed by each body is approximately 
an unknown linear combination of the fields which that 
body would scatter when responding to the plane waves 
in question, in isolation. The necessary coefficients can 
then be calculated by a consistency argument. For de- 
tails in the case of E\, see Karp;?* the argument is 
easily seen to apply for arbitrary boundary conditions, 
however. 

The theory in question can be applied to an infinite 
grating of cylinders. Before stating the results we recall 
some formulas from the general theory of infinite 
gratings. 

We consider a grating of congruent cylinders lying in 
the plane y=0. If a plane wave is incident at the angle 


¢, so that 
Dose = eik(e2 cos #+y sin ¢) (2) 


and if the distance between corresponding points of the 
grating is d, then the fields contributed by the various 
cylinders differ among each other only by a translation 
and by a phase factor exp (zknd cos ¢), so that the total 
field may be written 


U = tine + >. exp (iknd cos $)F(x — nd, y). (3) 


Suppose further that the far field form of F(x, y) is given 
by 


Fy ay ry eae (ier ay *) 0), (4) 
a kr 4 

where x =r cos 0, y=r sin @. Then in consequence of the 

regularity of the field in the exterior of a horizontal strip 

containing the grating elements, the total field can be 

written as follows 


2 y 
= Wine ta 


-f(46,) exp ikr cos (6 = 6,) (5) 
_~ kd sin 6, " p 


except for values of y in horizontal strips containing the 
grating elements. The upper sign is to be taken when y 
is above the strip. In this exact formula the spectral 
angles 6, are given by: 


k cos 6, = k cos 6 + 2nv/d 
k sin 6, = Wk? — (Roos ¢ + 2nv/d)?. (6) 


The infinite grating problem is thus equivalent to the 
determination of the values of the function f(@) for 6=8,. 
A finite number of the quantities 6, are real; these lead 
to propagating waves in (5);17.e., orders of the spectrum. 
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The function f(@) is expressible in terms of the boundary 
values of F(x, y) by quadratures which are meaningful 
also for complex 6. We can now state our preliminary 
formulas, with which we operate in the balance of this 
paper. We have found (cf. Appendix) using the type 
of analysis described in the foregoing paragraph, that 

F(8) = fr) + AxfO)fr(@) + Arf(r)frl@), (7) 

F(x, y) = Fy(x, y) we Axf(O)F 1(2, y) 

a Arf(m)Fr(«, y) (7a) 
where the functions f7(@), fr(@), and fr(@) are those de- 
scribing the far field of a single cylinder in the grating, 
when excited respectively by the incident field, a plane 
wave coming from x= — «, and a plane wave coming 
from x«=-+ 0. The function F7(x, y), etc., are the cor- 
responding scattered fields defined for all « and y. The 
constants Az; and Ap are given: 


‘xe wT ee 
A, = = FiCX Pe 
e A) = p( ae 


-exp inkd (1 — cos ¢)/./n 


DS T eS 
SSRN S.6 == 
Vie p( ae 


-exp inkd(1 + cos )//n. (8) 


The quantities f(0), and f(z), which are obtained from 
(7) by setting @=0 and 0=7 respectively, and solving 
the resulting pair of simultaneous equations, are given 
by the following: 


f(0) = {fr(0) — Aror}/{1 — [ALF L(0) + Arfe(n)] 

+ Ard ras} 
[fr(r) — Azae|/{1 — [4rf1(0) + Arfe(r) | 

+ ApAras}, (9) 


Ar = 


fa) = 


in which ay, dQ, a3 are given: 


_ | fr(0)  fr(0) _ | f2(0) fr(0) 
~ | fre) fel) SE SL GOY Ce) 
vie fx(0) fr(m) 
| f20) fer) | 


This result which is obtained by Karp,? holds with the 
proviso that k, d, and ¢ are such that there is no positive 
or negative integer v for which sin 6, vanishes. This im- 
plies that both A; and Ar converge. Then Az and Ar 
are small quantities for sufficiently large kd, so that the 
leading term of (7) is fr(0). The approximation f(6) ~f1(0) 
thus expresses the natural assumption that each grating 
element responds as if in isolation, when the spacing is 
very large. 


Oy} ’ ) 


(10) 


Ill. THe CAsE oF RESONANCE 


Our task in the present work is to discuss the behavior 
of the solution in the neighborhood of such frequencies 
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or angles of incidence as make sin 6,=0. We shall refer 
to such a situation as a “resonance,” and shall discuss 
the limiting value, of the result of (7), as well as the 
behavior in the neighborhood of the limit. (We shall not 
discuss the mathematical question of the existence and 
properties of the solution at precisely critical frequen- 
cies. Such a solution may well differ from the limit of 
the noncritical solution which it seems to us is the ob- 
servable quantity here.) 

The first thing we shall discuss is the limiting value 
of the solution given by (5) and (7). 

To do so we note that if v is such that sin 6, tends to 
zero, then the denominator of the vth term of (5) goes 
to zero. The corresponding exponential factor in this 
term becomes exp (+ikx). This analytical phenomenon 
is well known since Rayleigh, who pointed out the con- 
nection of grating resonances with the case when a 
spectral order passes off tangentially. Of course, the 
approximation f(@)~f1(@) breaks down completely in 
this case. However, we shall show that the use of (7) 
leads to an intelligible result. Our task will be to com- 
pute the limit of f(@,)/sin 6, as 6,0, 7. In addition it 
will be noted that all other terms of the series (5) require 
re-examination in this case. For one or both of the 
quantities A; and Arg diverge, and the quantities 
Axf(0) and Arf(ar) become indeterminate. 

In order to proceed it is necessary to scrutinize the 
manner of divergence of A; and Ar. This divergence 
can occur in several ways [cf. (8) ], 1) Rd(1—cos ¢) ap- 
proaches an integral multiple of 27, 2) kd(1-+cos ¢) ap- 
proaches such a value, 3) both case 1) and case 2) occur 
simultaneously. (For the last mentioned case to occur 
tan? (@/2) must be rational.) Suppose therefore that 
there exist positive integers, pa, p» such that either 


kd(1 — cos ¢) = 2pa7 + 6, (11a) 


og UO; 
or 
kd(1 + cos ¢) = + 2pom + 5, 


or both. In case 1) the denominator of the term vy = pu, 
vanishes in (5) and also Az diverges. In case 2) it is the 
term v= — f» which is critical while A diverges. It can 
be shown by comparison with the binomial expansion of 
(1—2)' that for small 64, we have 


ALR == 2/(2Rdbq,4)1!? 


5, > 0 (11b) 


(12) 


where the root is positive, or positive imaginary. In case 
3) (11) and (12) show that 


lim Axr/Ar = lim V60/ba. (13) 


If we keep @ fixed and let k approach a critical value, 
then the limit is cot ¢/2. If k is a suitably fixed value 
and @ approaches a correspondingly critical value the 
limit is +7. On the other hand, (6) allows us to express 
k sin 6, approximately in terms of 6, and 6, when y=p 
and y= —gq respectively. Introducing k,=k sin ¢,, we 
find 

Rpg — (2kd6,)1/2/d; 


Roy — (2hd8y)12/d (14) 
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so that 
lim Rogar == liven k_pAr as 
: y 
lim ky Ar — vat lim (Ap/A L) 
lim R_p Ay == lim (A L/Ar). (15) 


d 


Using (12) to (15) we can pass to the limit in (7). In 
case 1) we find 


{1(0) — Aray 
[@) —fr(0). + ——_—f@ 
f(8) flO): + eG f1(8) 
gitia @ 16 
Brea & Fit) 1 Ne (16) 


We noticed earlier that for large and noncritical kd the 
single scattering hypothesis f(@) =fr(0) supplies a first 
approximation to (7). This is no longer the case as a 
critical value is approached. Indeed, neglecting the cor- 
rection supplied by terms of order Ar =(1/~W/kd) in (16) 
we find 

fr(0) 


1G) > FO) Fi 1 (17) 


which differs from f;(@) by a finite amount. In case 2) a 
similar result is obtained. In this formula it should be 
noted that the denominator is the forward scattering 
amplitude of an isolated cylinder. Its negative real part 
gives the scattering cross section, by a well-known 
theorem. The vanishing of the cross section implies that 


if "| f2(6) 


This in turn implies f;,(@) =0 for all 6, whence also the 
vanishing (everywhere) of the field scattered by the 
cylinder when excited by the incident field exp zkx. The 
incident plus scattered field could not then combine so 
as to fulfill homogeneous boundary conditions. Conse- 
quently fz(0) cannot vanish. 

We now consider case 3). In this case we have to con- 
sider two possibilities: either a; vanishes, or it does not. 
Both cases are possible for it is easily verified, by calcu- 
lating the cross section explicitly, that for the case of a 
circular cylinder a3 cannot vanish. On the other hand 
the symmetry of the scattered field for a strip of zero 
thickness along the y-axis shows that in this case a3 does 
vanish. 

It is of interest to discuss the determinants a1, a and 
as. These are properties of a single cylinder. It can be 
shown that when a3 vanishes so do a and a. To see this 
define v=f1(0) Fr(x, y) —fr(0) Fx(x, y) which is an out- 
going wave function. Then form the quantity 
E(v) =Im fv(0d/dn)ds over the body. Then using the 


2dd = 0. 


Karp and Radlow: On Resonance in Infinite Gratings of Cylinders 


657 


boundary condition“ for fF, and Fr, and the well- 
known expression for the far field amplitude as a surface 
integral, we find E[v]}=+4 Re {F.(O)as 2 Consequently 
the vanishing of a3; implies E[v]=0. Since v is an out- 
going wave function, it now follows that v is identically 
zero. The argument is that used in the proof of the 
uniqueness theorem for diffracted fields, which actually 
requires only the vanishing of the energy flux E(v) over 
the surface. Having now shown that a3 =0 implies v=0, 
consideration of the far field form of v shows that 
fx (0) fe(@) — fr(0)fr(0) = 0 
for all @. Now by the reciprocity theorem for plane waves, 
we easily see that f7(0) =fr(@+7), fr(r) =f1(d+7), and 
fr(0) =fr(ar). Recall, for example, that f;(0) is the far 
field amplitude in the direction 0 produced by a plane 
wave incident in direction @ while fr(é+7) is the far 
field amplitude resulting in the direction @+7, when a 
plane wave is incident from the direction @=zero. 
Making these substitutions it follows directly that 
a,=0. The vanishing of a; implies [by direct substitu- 
tion of fr()fr(0)/fr(0) for fr(a)] that 
fr(O)as 
Oe = eee arn ty be 
fr(0) 
Consequently a;=0 also implies a2=0, as was to be 
shown. It seems worthwhile to inquire as to when a3 can 
vanish. The necessary and sufficient condition (for per- 
fectly conducting bodies) is that the boundary is com- 
posed of one or more vertical zero thickness strips, in 
the case of E||. In the other polarization any number of 
horizontal strips may be added. The spacing of the 
vertical strips must be according to half integral multi- 
ples of the wavelength. Thus we can have a rectangle of 
suitable width, in the H| case. In all these cases the back 
and forward scattered amplitudes agree, for E|l, and 
differ only in sign, for Hil. The result generalizes well- 
known symmetry properties of a screen of zero thick- 
ness. If a general homogeneous boundary condition is 
to be satisfied, one can derive a differential equation re- 
lating the equation of the curve and the coefficients in 
the boundary condition. These results all follow easily 
by applying the boundary condition. For the E|| case, 
for example, we have 


v = fr(0)Fe — fa(O)Fi = — [fx(O)e-** — fa(O)e**] 


on the boundary. If a;=0, v vanishes identically, hence 
x must be constant on the boundary. Taking the origin 
of x at a boundary point, we see fz(0) =fr(0). Hence at 
other boundary points sin kx =0, whence the spacing 
can be deduced. The sufficiency is shown by defining 


WV= Fy, = Fr 
in the case E||. Then W will be seen to vanish on the 
boundary and hence W=0, so that F,=Fr. Conse- 


quently f:(@) =fr(@) for all 6, whence it is easily shown 
that ay, a and a3 vanish. 


14 Thus, for Fl], Fz = —e7*** on the body. 
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Returning now to case 3) we first note the result when 
a3;%0. In this case, we obtain, in the limit 


f(0) = frl@) — (arfz(@) + a2fr() \/as. 


If a3 vanishes, then so do a; and a», and we obtain in 
general, using the relation of f,(@) and fr(@) 


[41 + Ar] [fr(0)f2(9) | 


(18) 


6) =e ir 19 
Oe apt ataty (19) 
and therefore in the limit 
fr(0) 
0 SSA UE 7] on i 7] e 20 
f(0) = fr(®) fio) 2 (20) 


We have now calculated the limit which f(@) ap- 
proaches." This allows immediate calculation of the 
complex amplitudes of the noncritical spectra, 7.e., the 
terms in (5) for which the denominator does not tend to 
zero. We still have to examine the terms vy =,, and/or 
y= —ps, where pa and fy are given by (11). 

In case 1) we find, using (15), 


2f (Oa) 


9 | 
Lit ee {in(o shin = (fr'(0)fx(0) — fr(0)fr"(0) — Arb coy} i, [Ara; — fr(0)] 


m . 
kd sin 0 ya 


where — (0) =asf1(0) —aifi(@)—a2fe(@) (note that 
(0) =0). Since kd is large in any case, this simplifies to 


_ 2 fre) Fr(0) 

eee ee 
kd sin 05, f(0) 

the result of (21) furnishes a correction to this formula. 

Thus the fotal coefficient in the last spectrum tends 


to a finite quantity as this spectrum tends to the 


horizontal. 
In case 2) we find, for the critical term v= — pp, the 
corresponding approximate value 


—fr(r)/fr(m) = — fr(w)/fr(0). 


Turning now to case 3), consider first the general 
situation, 7.€., a340. We find that the coefficient of the 
term vy =p, tends to 


(22) 


2f (Ona) a 2 B'(0 
(cee eee ee ee Os (23) 
kd sin Oia a3 kd a3 
Similarly, for »= — po 
se 2s Oy) ae 2 p(x) 
lim — Ra oe ae 
kd sin 0, as kd as 


which is again a finite result. 
As resonance is approached, in the case a3 =a, =a; =0 
we get instead, for y= pa 


edhe wa pas 
f1.(0)| 1 +%e 


1 In cases 1), 2), and 3). 
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plus a correction term in (kd)~!. Here a is the limit ap- 
proached by the ratio Ar/Ax, [cf. (13) and the subse- 
quent text |. 


IV. REMARKS ON ENERGY CONSERVATION 


For any wave function possessing the periodicity suit- 
able to the grating , and regular in the exterior of its 
elements, the total energy passing through the bound- 
aries of a typical cell will be equal to the energy passing 
through the surface of the typical cylinder contained in 
that cell, as can be seen from direct applications of 
Green’s theorem. Furthermore, due to the periodicity, 
the energy passing through the vertical boundaries of 
the cell is zero. 

Consequently, energy is conserved in the correct solu- 
tion just as well as in the incident field. For an incorrect 
solution, however, the degree of conservation of energy 
will be measured by the value of the surface integral of 
Poynting’s vector, provided we are dealing with a suita- 
bly periodic wave function, which is regular in the ex- 
terior of the cylinders of the grating. Now, with regard 


(21) 


to the single scattering solution, according to which each 
cylinder contributes to the field as tf in tsolation, it is easy 
to show that the boundary values are adopted with an 
error of order Az, say. Energy is conserved to order 
zero, except in the resonance case. In that case Ay 
diverges so that energy is not conserved at all. To see 
this, we simply note for example that a cylinder at a 
distance nd to the right of a given “central” cylinder 
contributes a field of the order of 


Ho(nkd) fr(m) 


to the boundary values on the cylinder selected as 
“central.” On the other hand, in our analysis, the con- 
tribution of a single cylinder 1s given by (7a). It is easy to 
show that in this case, the surface energy integral 
vanishes to order (kd)~'? and that the error is not only 
of higher order in (kd)~! for a noncritical case, but it is 
also uniformly bounded, 7.e., does not diverge as reso- 
nance is approached. Consequently in our result energy 
is conserved to the proper order, and this is independent 
of the “resonance.” 

It should be added that conservation of energy does 
not seem to guarantee that the boundary conditions are 
approximately fulfilled. Far from the vanishing of 
Im /uda/dnds over a surface, only the orthogonality of 
the function and its normal derivative follows. For 
example the exact solution for a grating of perfect di- 
electric cylinders is useless as a substitute for the solu- 
tion of the corresponding problem of conducting 
cylinders. Yet it certainly fulfills conservation of energy. 
While the energy law is useful as a check, proposed ap- 
proximate solutions always require additional support. 
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V. COMPARISON WITH OTHER THEORIES 


Since, as is well known, the theory of a grating of 
conducting cylinders implies a corresponding theory for 
the conducting reflection grating, it seems useful to com- 
pare the results obtainable from our solution with those 
by Artmann. Artmann’s results are obtained under the 
assumption of protuberances small compared to wave- 
lengths, without restriction as to spacing. By specializ- 
ing his results to the case of large spacing, while 
specializing our results to the case of thin cylinders, a 
valid comparison is possible. However, we require ex- 
plicit expressions for the quantities fz(0), etc., since 
Artmann’s results are expressed in terms of the geome- 
try of the surface. To this end, we have calculated, by 
the method of small perturbations, the pattern produced 
by diffraction by a thin strip of thickness y=2ef(x) 
where ke is small. The procedure is omitted here for 
brevity; only the results are given as follows. 

We consider a symmetrical strip y= +ef(x), for 
| «| </. The incident magnetic field has a single com- 
ponent, which is perpendicular to the xy-plane and is 
given by Fine =exp tkx. Then if we write 


H(x«, y) = exp tka + F(x, y) 


Fila, y) = F(x, y) + OF, (x, y) +--+ (24) 
we find 
k y 
F(a, 9) => f 
-e*®20f" (49) Hy (Rv/ (x — a0)? + y2)dxo (25) 


i? OF (x0 
NE 2 eres {i"(o) Bea 
2 Ey 0X0 
OF 1, (x0,0) 
— f(%o) 5 
Oy? 


{Ho (kV (x — 0)? + y?} dao. (26) 


These results are obtained by expanding the boundary 
condition for /; in a power series ine. It is easily verified 
that the function exp ikx+eF,%+eF,™ fulfills the 
boundary condition to order e?. From (25) and (26) we 
find 

ke f°! 
ful) == ff (aeitoeer Md + Of) (27) 
Ld_j 
The leading term vanishes for 6=0, since f(+/) =0. 
This is natural since Re f,(0) is quadratic in e€ by the 
forward amplitude theorem. To obtain fz(0) we use 
(26) together with (25). After various integrations by 
parts we obtain 


ae, he? ee) | (1) 2 
ia (0) i dt O( Re)? 28 
1(0) 7a ie ip + O(Re) (28) 
where 
1 
g(t) = f' (xo)e®O—) dry, (29) 
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Here we have used the expression 
H™(|x])=— f° a (30) 
tid. Vt—1 
where 
Vi—1 = = iVl—F- for. i] <1. 


It is easy to verify that Re fz(0) is negative, while 
Im fz(0) is positive. 

Let us now briefly discuss the solution of the reflection 
problem. It is well known that to obtain this solution 
from the cylinder solution we have to consider a grating 
of cylinders symmetric with respect to the plane y=0, 
and construct the function u(x, y; ¢)+u(x, y; —@), 
where u(x, y; @) is the solution corresponding to the 
incident wave exp 7k(x cos ¢+y sin @). The upper sign 
is taken when we have H , while the lower sign refers to 
the case of E\|. In calculating by means of (7), (8), (9), 
(10) we have to consider the subscript J to refer to both 
the directions of incidence ¢ and —@. If ¢ is replaced by 
—d¢, then A, and Apr are unaffected. Because of the 
symmetry, the quantity f7(@) for the direction of inci- 
dence —¢ is equal to fr(—8) as calculated for the direc- 
tion of incidence +¢. Consequently also f7(0) and fr(7) 
remain unaltered under this transformation, whence the 
same is true of a; and a of (10). It follows that in calcu- 
lating with (7) for the direction of incidence —@, the 
first term, f7(0) is replaced by fr(—#@) and the remaining 
terms are unaltered. Consequently, for E||, the neces- 
sary subtraction leads to a simple value for the total 
f(@) in the reflection problem, 7.e., 


Feotar() > fr(@) aa fr(—8): 


Since fr(@) is an analytic function of cos 6, it is easily 
verified that there is no divergence of the terms 
f(0,)/sin 6, as 8,30, in (5), so that there is no resonance 
phenomenon. 

On the other hand, in the case of H]|, we have to add 
the terms f(@)+/(—@). It is interesting to note in this 
connection that, if the cylinders are taken to be of rec- 
tangular cross section and of unit or integral wave- 
length thickness, then it follows from our discussion of 
the determinants a in a previous section that fr(@) 
+fr(—6) =0, if we choose ¢=7/2. In this case our re- 
sult for the reflection grating is the trivial result 
exp iky+exp (—zky), regardless of whether the spacing 
is critical. However this result is easily verified for any 
collection of rectangular proturberances, for any spacing 
whatever, if ¢=a/2 and the thickness is half a wave- 
length. A protuberance of half a wavelength corre- 
sponds to a cylinder of thickness equal to a wavelength. 

We continue our discussion of the H]|| case by special- 
izing to the case of thin cylinders. For strips of zero 
thickness the far field amplitude is well known to be 
odd in 6. Consequently the term fr(@)+fr(—@) will be 
of order ke where ¢ is the thickness. If we turn to the 
next term, however, in the case of resonance [cf. (17) | 
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we find a different order of magnitude. For the reflection 
grating this term becomes in fact 


[f.(0) + fr(—8) |. 


Now f7(0) =fr(@+7) by reciprocity. For zero thickness 
this quantity is identically zero and for a thin strip it is 
of order ke. The quantity fz (6) +fz(—8@) is also of order 
ke. But the denominator f;(0) is of order (Re)? [cf. (27) ]. 
Therefore the total amplitude is a finite number, no 
matter how small ke is, provided we have a critical com- 
bination of incidence angle and wavelength. This is in 
contrast to the situation for noncritical wavelengths. 
Artmann has shown also that the maximum amplitude 
occurs at slightly larger wavelengths than critical, and 
that the shift is independent of the order of the spec- 
trum. To see this from our result we return to the for- 
mula for f(@,) at a noncritical wavelength. Consulting 
(7) and (9), and neglecting terms of order Ar, we find, 
for a case when A, is large but not infinite, 


fr(0) fr.) 
(yA 7 ener O) 


From (12), A,~2: (kd6,)—!/?, where, if 6 is negative, the 
positive imaginary root is taken. Since —Re fz(0) is 
positive by the forward amplitude theorem there can be 
no cancellation between 1/A,z and f,(0) if 6 is positive. 
On the other hand, if 6 is negative, then 1/A z is positive 
imaginary. Comparing (28) we see that Im fz(0) is also 
positive. Thus the imaginary part of the denominator 
will vanish at a certain negative value of 6., given by 


V2kd | 6a | = —f- | o(0) —— Ki 


S(O») + f(—8) = fr) + fr(—%) ++ - (31) 


(32) 


Consequently [cf. (11a)], a maximum occurs for kd 
slightly less than critical, in general, since the numera- 
tor will only vary slightly with such a small variation in 
k. Eq. (32) shows that 6, is independent of the order », 
in which the maximum is being observed, and is also!” 
independent of the index p,. of the order which tends 
toward grazing. Furthermore it is easily seen that if f. 
is the amplitude at the critical value, while f,, is the 
value at the maximum then (f.—fm)/6 is of the order of 
f./6, 1.e., inversely proportional to 6. Consequently the 
maximum is very steep. 

It is possible to show that all these results are in 
agreement with the results of Artmann [compare for 
example his equation (3.32)], provided we replace 
Artmann’s sums by integrals, in accordance with the 
fact that kd is to be assumed large. Using our notation, 
Artmann’s result can be rewritten 


16 Then the leading term is f7(@) +f1(—6), which tends to zero with 
the thickness. 

17 Tf we keep kd fixed, and vary the angle of incidence in order to 
vary pa. 
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0 = V5.(2kd — 5.) + ilm ») “(= : “) (2) 


VAPa d 


f(%) exp ix(* + kcos6, — ) a: ES Kamar 


=) 


Using the fact that d is large we can replace the above 
sum by an integral, via the substitution k cos 6,=Rt, 
Ir/kd=dt, —» <t<o. Since 6, is small we have 
In (v—pa)/d=kt—k, while the term 6,/d in the ex- 
ponential can be neglected. This results ‘precisely in our 
(32). On the other hand Artmann’s theory is also valid 
for small thickness of protuberances, even if the spacing 
is not large, while our results, although restricted to 
large spacing (kd>>27r), do not require that the thickness 
is small compared to wavelength. 

Having briefly compared. our theory with that of 
Artmann we now turn to a brief comparison with new 
results of Twersky,!? which appear to be in disagree- 
ment with ours. To summarize this disagreement, we 
note that a functional equation is obtained for the 
amplitude in the reflection grating. This functional 
equation is solved approximately by iteration, the 
justification for retention of the leading term (alone), 
being the assumption of large spacing.!8 Then the result 
is a quantity f,=/(0,0, * —0) which is in turn defined as 
the solution of a new integral equation [equations (60), 


(56) ] 


foo = F@, 7 — %) = fO,, 7 — 60) 
r/2 A 
eo Bas > f(8, « — $) f(b, + — A) dd, 
dad _x/2 


in which the quantity f is the result based on single 
scattering theory, and f is to be solved for, with f as 
given data.!® It is shown that f is, under certain circum- 
stances, purely imaginary, and that consequently if one 
retains terms of higher order of magnitude in the 
specular amplitudes, while neglecting them in the other 
spectra, then energy is conserved exactly.?° It is then 
emphasized”! that the results (which according to page 
29 and (70) are obtained for large kd) “differ significantly 
from a single scattering approximation in that they con- 
tain f instead of f.” Arguments are given to show that 
the presence of f rather than f would be incorrect 
physically. The results insofar as they are explicit seem 
to be in disagreement with our own even before the 
passage to the limit of resonance, since our naive leading 
term is the single scattering term, in the noncritical 
case, and our principal result is the correction to this 
term so that in the limit of infinite spacing our results 
reduce to single scattering. We reserve further comment 
until Twersky’s results have appeared in final form, but 
note that the discrepancy would be removed if the 
validity of Twersky’s result were additionally restricted 

18 Twersky, op. cit., cf. (70) and p. 11. 

19 In this formula the superscribed bar does not mean complex 
conjugate. 


20 Twersky, nee Clie 


pp. 29-30. 
4 Tbid., p. 30. 
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to smali scatterers, or at any rate to scatterers whose 
scattered amplitude is small, for then his integral 
equation shows that ff. 


APPENDIX 


In this appendix we give a brief intuitive derivation 
of (7) and (7a). 

The mth element of the grating contributes an out- 
going wave function F(x—mnd, y) to the total field 
[cf. (3) of the text]. The function F(« —nd, y) is regular 
in the exterior of the mth cylinder. If 7,, 6, are polar co- 
ordinates measured from the origin (nd, 0), then we 
know that for large r, the asymptotic form of F is 


F(x — nd, y) > (2/rkr,)1!? exp tk @ —1 =) ie (33) 
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where f(6,,) is unknown [cf. also (4) of the text ]. Now if 
the spacing is large compared to both wavelength and 
dimensions of the cylinders, then the zeroth cylinder 
subtends an infinitesimal angle at the other cylinders. 
Hence, for cylinders to the left of the zeroth cylinder, 
(1.€., m negative) 0,—0 in the neighborhood of the zeroth 
cylinder, while for 1 positive 0,—7. On the other hand 
if («, y) is a point in the neighborhood of the zeroth 
cylinder, then 7,(x, y)—|« —nd] . Using these approxi- 
mations in (33) and then evaluating (3) of the text in 
the neighborhood of the zeroth cylinder, we find that 


u = tine + Arf(O)et*? + Anflme *** + F(x, y). 


Thus the effective excitation contains two plane waves 
in addition to Wine). Now F(x, y) is the sum of the re- 
sponses to all these excitations, whence (7a) of the text 
follows. Eq. (7) is the far field form of (7a). 


Line-of-Sight Wave Propagation ina Randomly 
Inhomogeneous Medium* 


BOB M. FANNIN{T 


Summary—Theoretical calculations have been made, using the 
single-scattering approximation, for propagation in a randomly in- 
homogeneous medium in which the deviations of refractive index 
from the mean are small. The statistical quantities considered were 
the variance, correlation function, and power spectrum for the phase 
and relative amplitude of the field at a point and their differences 
at two points. The emphasis in this paper is in indicating the transi- 
tion from the ray treatment results to the scattering cross section 
results. The correlation function for the refractive index is taken to 
be time as well as space dependent so that tiie power spectrum can 
be computed from the original formulation. 


INTRODUCTION 


HE DEGREE of inhomogeneity of the atmos- 
apse constitutes the limiting factor in the ulti- 

mate accuracy of radar systems. The quantity 
which describes the propagation characteristics of the 
atmosphere at a given point is its refractive index 
(complex if absorption is considered), which is a func- 
tion of the temperature and pressure of the air as well 
as the concentration of such constituents as water 
vapor, precipitation, and free electrons. The refractive- 
index distribution in the atmosphere is such a compli- 
cated function of space and time that, when making a 
theoretical study of atmospheric propagation phe- 


* Manuscript received by the PGAP, February 13, 1956. The re- 
search reported here has been sponsored by the AF Cambridge Res. 
Ctr., Air Res. and Dev. Com., under Contract AF 19(604)-494 and 
St. Louis Ordnance District, under Contract DA 23-072-ORD-763. 

+ Elect. Engrg. Res. Lab., University of Texas, Austin, Tex. 


nomena, it is generally necessary to consider one at a 
time the different characteristic types of variations 
likely to be encountered. This study is concerned with 
the relatively small scale, random-type refractive index 
variations that are present in the atmosphere primarily 
due to its turbulent character. The stratified nature of 
the atmosphere and such phenomena as inversion layers 
are not considered in this paper. 

The situation being considered is one in which the 
propagation path is line-of-sight and such that the 
earth and all other obstacles are far enough removed so 
that their effect can be neglected or considered sepa- 
rately. Computations of quantities associated with a 
wave traveling through a randomly inhomogeneous 
medium such as the atmosphere (in which the root- 
mean-square of the refractive-index deviations from the 
mean is very small) are usually carried out on the basis 
of a scattering mechanism! or a ray concept.” Not only 
are the fundamental physical mechanisms visualized 


1C, L. Perkeris, “Note on scattering in an inhomogeneous me- 
dium,” Phys. Rev., vol. 71, p. 268; February, 1947, H. G. Booker 
and W. E. Gordon, “A theory of radio scattering in the troposphere,” 
Proc. IRE, vol. 38, pp. 401-412; April, 1950, are references on basic 
scattering theory, while A. D. Wheelon and R. B. Muchmore, 
“Line-of-sight propagation phenomena—II. Scattered components,” 
Proc. IRE, vol. 43, pp. 1450-1466; October, 1955, deals specifically 
with the line-of-sight problem. 

2 See, for example, R. B. Muchmore and A. D. Wheelon, “Line- 
of-sight propagation phenomena—I, Ray treatment,” Proc. IRE, 
vol. 43, pp. 1437-1449; October, 1955. 
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quite differently in connection with these two ap- 
proaches but the mathematical formulations involve 
different approximations which are valid for contrast- 
ing situations. In short, the ray concept is generally 
considered to be valid when the atmospheric “blobs” 
are relatively quite large so that the receiver is in the 
near (Fresnel) field of the blobs while the scattering 
cross section approach is appropriate when the blobs 
are small so that the receiver is in their far (Fraunhofer) 
field. The treatment employed here is that obtained by 
applying first order perturbation theory to the wave 
equation; 1.e., the single scattering approximation.’ The 
scattering cross section concept generally employed 
evolves when the additional assumption is made that 
the blobs are sufficiently small compared with the range 
and wavelength. This last restriction will not be imposed 
in this development so that the results are equally valid 
for large as well as small blobs.* This approach thus has 
the advantage of presenting a unified treatment which 
furnishes results that, for the limiting cases of very 
large and very small blobs, coincide with those obtained 
by applying the usual ray or scattering cross section 
methods, respectively. 


REFRACTIVE-INDEX CORRELATION FUNCTION 


The following development incorporates a correlation 
function for the refractive index which is both time and 
space dependent. This makes it possible for power 
spectra calculations to be made directly from the 
original formulation. 

The medium is taken to have small random variations 
about a mean, the statistical character of the fluctua- 
tions being homogeneous and isotropic except for a mean 
wind. Thus the refractive index (7) can be expressed as 


n(x, t) = no + An(x, t), An<<<m, m1 (1) 


in which 7p is the mean value of refractive index, x de- 
notes a position in space, and ¢ denotes time. Also, using 
boldface symbols to denote vectors, 


An(x1, t)An(a2,t-+ 7) = An?C( | CT a) (2) 


where r is the vector distance between x; and xs, v is the 
mean wind velocity, C is the normalized time and space 
dependent correlation function for the refractive index, 
and the bars denote time averages. The effect of absorp- 
tion is neglected so that An? is real. 

In what follows the form of the correlation function 
will be taken as® 


’ D. Mintzer, “Wave propagation in a randomly inhomogeneous 
medium.I,” J. Acous. Soc. Amer., vol. 25, pp. 922-927; September, 
1953, or A. D. Wheelon, “Near-field corrections to line-of-sight prop- 
agation,” Proc. IRE, vol. 43, pp. 1459-1466; October, 1955. 

4L. A. Chernov, “Correlation of the amplitudes and phases in 
wave propagation in statistically-inhomogeneous media,” Doklady, 
Akad. Nauk USSR, vol. 98, pp. 953-956; June, 1954, treats a very 
similar problem by almost the same methods presented in this 
paper. 

5 The author is indebted to K. A. Norton for pointing out the 
equivalence of the two expressions in (2a). 
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C = exp [—| r — wr}?/P] 


=exp {—[wr?+ |r—- vr || /12} 


in which w denotes the wind velocity composed of a 
mean component v and a component of magnitude w 
with random direction while / is a measure of the scale 
of turbulence and is frequently referred to as the mean 
blob size. Thus incorporated in (2a) is a Gaussian de- 
crease with separation for 7 =0, a drift with velocity v, 
and a Gaussian decrease with 7 to take into account the 
random movement and dispersion of ‘the blobs. Un- 
doubtedly (2a) is not an entirely adequate representa- 
tion of turbulent atmospheres. There is, however, little 
direct data upon which to base a more appropriate cor- 
relation function and, in at least some respects, the re- 
sults for the line-of-sight case are relatively insensitive 
to the exact form of C. The-additional assumption that 
I> is also made. 


l| 


(2a) 


MATHEMATICAL FORMULATION 


Let us consider an isotropic radiator (with wave- 
length \) at point (a) in space and study the field at 
point (b), a distance R from (a). The wave equation for 
which a solution is sought is 


(V2 + kn) E = 0 Ie 


where k=(27/X). Neglecting terms in (Az)?, defining 
E, by E=(E)+E;) where Ej is the field that would exist 
if An =0, and taking mp» as unity, (3) then becomes 


(V+ B)E, = —2PVnE (4) 


In this form the equation is immediately interpretable 
as expressing E, as being composed of contributions 
from each differential volume (dv) in space radiating like 
an elementary dipole whose moment and orientation is 
given by —2AnEdv. Eis not known so we take the single- 
scattering approximation by using Ep instead of the 
correct E for the field at the differential scattering 
volume. When the intuitive assumption is made that the 
important scattering volumes are relatively close to the 
straight-ray path (since />\) the sine of the angle be- 
tween the direction of Eo at the scattering volume and 
the direction to the receiving point can be taken as 
unity. Also, only the component of the field at the re- 
ceiver that is in the direction of Ey need be considered 
so that the vector notation can be dropped. Considering 
the E’s (not boldface) to denote complex quantities in 
the customary convention for sinusoidally varying 
fields, the single-scattering expression for (E,/E,) be- 
comes 


(Ex/ Eo) = (?R/2n) if A(n/RaRs) 


-exp [—jk(Ra + Ry — R)|dv (5) 


in which R, and R, are the distances from dv to a and 8, 
respectively. In keeping with the small perturbation 
approach employed throughout, the real and imaginary 
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parts of (5) are taken as the deviations in magnitude 
and phase of E relative to Eo. Since An is defined to be 
zero the average value of (5) is also zero. 

The fundamental quantities to be determined are 


Pm,p(D, T) = Om,p(b1, t)Qm,p( be, y = T) (6) 


in which a,(6, t) and an(b, tf) denote the phase of E/E, 
and the magnitude of (E—E»)/Eo, respectively, at 
point 6 and at time ¢. Eq. (6) define two quantities, the 
comma between the subscripts indicating that either the 
first or the second subscripts may be taken. Both 6, and 
bz are assumed to be a distance R from the transmitter 
at (a) and displaced one from the other a distance D, the 
assumption being made that R>>D. The p’s are thus the 
cross-correlation functions for the phase or relative 
amplitude variations of the field at point 0}; and bo, D 
being a parameter. Of course, when D =0 the p’s become 
autocorrelation functions. 

From (5) it follows that 

An(«, t1)An(x’, te) 


Dy a 1/2027 RIX ale 
een PIPER RER YR 


: [cos R(Rg a R, =< Re = R,’) 
+ cos R(Ra + Rp + Ra’ + Rs! — 2R)|dvdv’ (7) 


the upper sign being taken for p, and the lower sign 
for pp, the primed and unprimed symbols referring to 
two different points in space. To evaluate this double 
volume integral, spheroidal coordinates were used and 
(2) and (2a) employed. Again presuming that the im- 
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The component of v along the propagation paths has 
been taken as zero and ¢ denotes the angle between v 
and the direction from 0; to bs. The parameter g is 


g = R/4nl? = 1y2/xP? (11) 


in which /) denotes the half-width of the first Fresnel 
zone at the midpoint of the path. 
The variable of integration, y, in (9a) and (9b) is 


y = (Ra + RR.’ — Ry — Ry’)/2R (12) 


so that these integrals denote the summing up of the 
contributions of different layers sliced normal to the 
propagation direction. Of course, in these spheroidal 
coordinates the layers are hyperbolic but they are ap- 
proximately plane over the important region. Therefore, 
(9a) and (9b) are appropriate for considering the effect 
of a turbulent layer. Assumed variations along the path 
in the intensity or scale of the refractive-index fluctua- 
tions can be inserted into these equations. Also, the 
transition to the case of a plane wave incident upon a 
turbulent region (R-—> ~~) is easily effected. 

Continuing with the assumption that the statistical 
character of the refractive-index fluctuations is homoge- 
neous throughout the whole region, (9a) simplifies to 


P(D, 1) = (I/D) exp [—(w? + v sin? $)7?/P | 
: { erf [(D — or cos ¢)/l — erf [v7 cos ¢/1]}. (9a’ ° 


The autocorrelation functions, Ym,»(D, 7), for the dif- 
ferences in phase or relative amplitude at points 0; and 
bs, are given by 


Ym, p(D, Tt) = Lom, p(d1, t) 


portant scattering volumes are relatively near the 
straight-ray path and that R>>l, it was assumed that for 
the correlation function in the integrand, approxima- 
tions could be introduced that amount to taking the 
spheroidal coordinates as cylindrical over a volume of 
dimensions of the order of J. With these simplifying 
modifications (7) reduces to 


pm,p(D, tT) = An®x!21R\—[P,(D, 7) F P2(D,7)] (8) 
in which 


P\(D, 7) = exp [—(w? + 0? sin? ¢)r?/??] 


Lf exp {=[01 + 9/20 = or/t) cos 6] dy (90 


1 
¥ 
P,(D, 7) = exp (—#7?/P) Re | Q-texp {—[1+ y)?D?/4 
=i! 
— (1 + y)Dor cos ¢ + v?r?]/P?O} dy (9b) 
and 


Q = [1 — 2¢(1 — yl: (10) 


as Am, p(b2, t) | [om p(b1, t+ T) 
2pm p(0; T) mit PmipD, T) op pm,p\D; —r)]. 


ae Om, p(b2, t oe 7) (e 
(13) 


Also of interest are the power spectra, S»,,(D, w), of the 
phase or relative amplitude differences, 


ico) 


Sm,p(D, w) = (in) f Ym,p(D, 7) cos (rr)dr (14) 


oO 


in which S»,,(D, w) has been defined such that 
[ Smp(D, 0d = Ym, 9(D, 0) 
0 


The portion of (8) that involves P; is seen to differ 
from the expression for p,(D, 7) derived using a ray con- 
cept only by a factor of one half.6 Thus the functional 
form of P; is independent of g while P» is a function of g. 
A brief inspection of (9b) and (10) reveals that P,P; 
as g—0 and P,—0 as g~. Thus for very large blobs 
Pm~0 while p, has the value that would be calculated 
using a ray treatment. On the other hand, for very small 
blobs, pm», the functional dependence upon /, Rk, and 
d differing from that for p, for large blobs only by a 
factor of one half. Since y(D, 7) and S(D, mw) are defined 


6 Muchmore and Wheelon, op. cit. 
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in terms of the p’s, the statements concerning the p’s in 
this paragraph apply equally well to these quantities. 

Also in (9a) and (9b) it is seen that the “dispersion” 
velocity, “, enters only in the form of the same Gaussian 
multiplying factor that appears in (2a), since this factor 
is not involved in the space integrations. Thus any other 
functional form involving u could be directly substi- 
tuted. Likewise, if it was desired to let / in (2a) be a 
function of 7 this time dependent could be substituted 
throughout. 

Since pm,p( ©, 7) =0, (13) shows that y(©, 7) =2p(0, 7) 
as would be expected. This simply indicates that the 
signals at 6; and by are uncorrelated when D is very 
large. Thus the spectra for the fluctuations of a single 
signal are given by $Smn,»(%, w). 


NUMERICAL RESULTS 


Space limitations preclude the presentation of the 
numerical results for more than a few combinations of 
the values of the many parameters. 


Single Path Results 


The variance of the phase and relative amplitude of 
the field at a single point is just pm,»(0, 0). From (8), 
(9a’), and (9b) this is seen to be 


J 
pm,p(0,0) = sR 1/2 + 1/2Re i) oay| (15) 
0 


in which Q is defined by (10). The integral in (15) can 
be evaluated and the real part taken but the resulting 
expression is so lengthy that it will not be given here. 
The quantity in the brackets of (15) is plotted vs qg (the 
q scale is logarithmic) in Fig. 1. The coefficient of the 
bracketed expression in (15) is the value of p,(0, 0) that 
is obtained using a ray approach. It is seen that p,(0, 0) 
agrees with the ray concept value for big blobs but 
differs from it by a factor of one half for small blobs. 
Pm(0, 0) goes from zero to the same value as p,(0, 0) as 
the blob size goes from large to small as mentioned in 
the previous section. Fig. 1 shows the transitions as q 
goes from a small to a large value. 

The square root of p,(0, 0), (@p)rms, thus varies di- 
rectly as the square root of / and R inversely as \ when 
g is either much larger or much smaller than unity, with 
the dependence on these parameters shifting to /~°-°78, 
DPriS? AN Na On Gg 1, “Line Value, Of) (Qn) smart alsO 
varies as the square root of / and R and inversely as 
for large g but approaches (8./7An?/15)°*(R/1)15, in- 
dependent of A, for small g. That is, when a scattering 
mechanism predominates then (@p)rms and (Qm)rms have 
the same value and dependence upon J, R, and but 
when a ray concept suffices then the value of (m)rms iS 
much less than that for (@»):mg and its functional de- 
pendence upon /, R, and d is considerably different. 

Eqs. (8), (9a), and (9b), with D taken to be zero, give 
the expression for the autocorrelation functions for 
phase and relative amplitude variations at a single re- 
ceiving point. As mentioned above, the dispersion ve- 
locity, u, enters solely in the form of a Gaussian multi- 
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plier while the drift velocity, v, enters in a more complex 
manner which depends upon the value of g. This de- 
pendence of pm,p(0, 7) upon v is shown in Fig. 2 for 
qg=0, 1, and «. The curves in this figure were normal- 
ized by dividing by pm,p»(0, 0). 

The effect of the parameter g upon the power spectra 
for the phase and relative amplitude variations at a 
single point can be determined qualitatively by inspect- 
ing Fig. 2. The correlation functions for the phase for 
g=0 and - and for the amplitude for g= © are Gaus- 
sian so the corresponding power spectra will also be 
Gaussian. The correlation function for the phase for 
q=1 falls off less rapidly so that a larger portion of the 
power is contained in the lower frequencies. Likewise, 
as g decreases the correlation function for the amplitude 
variations falls off more and more rapidly so that a 
larger portion of the power is contained in the higher 
frequencies. 


Cross-Correlation Functions of Signals at Two Receiving 
Points 


Fig. 3 shows pm,p(D, 0) (normalized) vs (D/l), 7.e., 
a plot of the normalized cross-correlation functions for 
the phase and relative amplitude variation of the field 
at 6; and bs, for zero time difference, as a function of the 
separation between receiving antennas. 


Signal Differences for Two Receiving Points 


The power spectra for the phase and relative ampli- 
tude differences are plotted in Fig. 4 for a number of 
values of (D/l) for gq= ©, u=0, and ¢=0. That is, the 
curves are for the case when only the effect of a mean 
wind across the paths is considered. Naturally these 
curves apply equally well for phase differences when 
g=0. It is seen that decreasing the antenna separation 
reduces the rms of the fluctuations but cuts out the 
lower frequencies more rapidly than the higher fre- 
quencies. 

Fig. 5 indicates the effect of g upon the power spectra 
for the difference variations. These curves have been 
normalized and are for (D/l) =1, u=0, and ¢=0. 
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Partially Reflecting Sheet Arrays* 
GISWALT VON TRENTINI{ 


Summary—Multiple reflections of electromagnetic waves be- 
tween two planes are studied, and the increase in directivity that re- 
sults by placing a partially reflecting sheet in front of an antenna with 
a reflecting screen is investigated at a wavelength of 3.2 cm. The 
construction and performance of various models of such arrays is 
discussed. Thus, for example, a ‘‘reflex-cavity antenna’? with an 
outer diameter of 1.88 \ and an over-all length of only 0.65 » is de- 
scribed which has half-power beamwidths of 34° and 41° in the E 
and H planes, respectively, and a gain of approximately 14 db. It is 
shown that larger systems produce considerably greater directivity 
but that their efficiency is poor. 


INTRODUCTION 


A LARGE, plane, conducting screen placed behind 


an antenna serves as a shield against backward 

radiation and also affects, depending on its 
spacing, the forward pattern.! Improved directivity in 
the normal direction is possible only to a limited degree 
because the illumination of the reflecting screen is not 
optimum and the rays reflected from more distant zones 
produce, in part, out-of-phase contributions. 

Improved illumination and an increase in directivity, 
and therefore in gain, can be obtained by adding a par- 
tially reflecting sheet in front of the antenna and parallel 
to the reflecting screen, causing multiple reflections be- 
tween the sheet and screen. The distance between the 
sheet and screen must be such that the partial rays pro- 
jected through the sheet into space have equal phases in 
the normal direction. 

This paper presents the results of tests with various 
types of sheets and discusses a few simple forms of 
“reflex antennas.” In order to work with small systems 
that are easy to construct and convenient to measure, 
the investigation was carried out at a wavelength of 
3.2 cm. If adequately scaled, such arrays may also be 
used for longer wavelengths. 


DERIVATION OF THE RADIATION PATTERN OF AN AN- 
TENNA BETWEEN SCREEN AND PARALLEL SHEET 


If an antenna is located in front of, or in the plane of, 
a conducting screen, its radiation can be imagined as 
originating in a point P of the screen with element pat- 
tern f(a) (Fig. 1). 

A partially reflecting sheet placed at a distance / from 
the completely reflecting screen introduces multiple re- 
flections with decreasing amplitudes between these two 
planes. Let the reflection coefficient of the sheet be 


* Manuscript received by the PGAP, April 18, 1955; revised 
manuscript received, October 3, 1956. 

t Starnberg, Oberbayern, Luitpoldweg, Germany. 

1 John D. Kraus, “Antennas,” McGraw-Hill Book Co., Inc., 
New York, N. Y., ch. 12-2; 1950. 


Fig. 1—Multiple reflections between screen and sheet. 


pe’¥. Assuming no transmission losses, the amplitude of 
the direct ray 0 is proportional to /1— ?; the ampli- 
tude of the once-reflected ray 1 is proportional to 
pV1—p*; the amplitude of the twice-reflected ray 2 is 
proportional to p?./1— 2; etc. The electric field in- 
tensity in the Fraunhofer zone consists of the vector 
sum of these partial rays, and for an infinite screen and 
sheet we may write 


E = DS fla)Eop"V1 — prer®r. (1) 


n=0 


The phase angle ©, is composed of the phase variations 
during reflections from the completely reflecting screen 
and partially reflecting sheet, and also of the path dif- 
ferences of the partial rays. An additional phase shift 
occurs when a partial ray passes through the sheet, but 
since this is the same for all rays it need not be taken 
into account. From Fig. 1 we may derive the phase dif. 
ference between ray 1 and ray 0 as 


1956 


7 : Mage Dah 
QO, = —2/tanasina — — = mt = W, 
nN A cosa 
and between ray 2 and ray 0 as 
Ue : 2x «Al 
@Q. = —4ltanasina — — —~—27+ 2p 
r A COS a 
which, somewhat transformed, gives 
Ar 
@, = n= | — te0sa— 2 +y | (2) 
Since p <1, we obtain 
( ei?) $ {= Se RE . 


Inserting this in (1), the absolute value of field strength 
becomes 


alee |. sla) 4/ ete ee 
1+ p?— 2pcos ® 


and the power pattern is therefore 


1 — 7 


S= 


Cale) 
Aor 
Wis eee 2pcos (y= == cos a) 


We must take into account, however, that the amplitude 
pb and the phase y of the sheet reflection coefficient are 
a function of the angle of incidence a. Maximum power 
in the direction of ~=0° is obtained when 


j Ly ; 
Sp! ES P= '0! 
r 


and hence the equation determining the resonance dis- 
tance J, of the sheet is 


Wo ON r 
eg ar eee 3, 4 
; ( \o+ . (4) 


with Wo expressed in degrees and N=0, 1, 2, 3, etc. The 
phase angle for a sheet with inductive field impedance 
is in the second quadrant, and for a capacitive sheet, in 
the third quadrant. 


RESULTS OF TESTS AND COMPARISON 
WITH CALCULATIONS 


An exact evaluation of (3) and (4) is difficult since the 
reflection factor of sheets is known for the case of a 
plane wave only. Using the plane wave value, however, 
an approximate, qualitative solution can be obtained 
that should be satisfactory in the Fraunhofer zone. 

If the sheet consists of a grid of parallel equidistant 
wires running in the direction of the electric vector, the 
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field impedance will be: 


iV = j120 “cos a| in + F | (5) 
oN 2p 
The factor F is itself again a function of a, but may be 
neglected if the wire distance d (grid constant) is small 
with respect to A, and the wire radius p is very small 
with respect to d. 

From the relation 


—1 
iy = 6 
pe ae = (6) 
J 607 
we get 
ee TTR AV ANE (7) 
V+ (a) 
607 
and 
V 
y = arctan (- =<) (8) 
Oa 


For a grid with 20 =0.5 mm and d=9 mm, and with 
\=3.2 cm, we obtain po=0.707 and ~).=135° in the 
normal direction (a#=0°). In other words, at the first 
reflection 50 per cent of the energy is reflected back to 
the screen. From (4) we find that the resonance distance 
l,is 14mm, and from (3) we obtain the radiation pattern 
in the H plane, as shown in curve 1 of Fig. 2. The ele- 
ment pattern of the antenna, fia”, is taken from the 
measured values for a waveguide aperture with a large 
reflecting screen, as given in curve 1 of Fig. 3. The total 
pattern obtained in this way has no side lobes, a 3 db 
beamwidth of 35°, and a 10 db width of 58°. The calcula- 
tion for a corresponding capacitive sheet gives a reso- 
nance distance of 18 mm and a slightly more directive 
radiation pattern. A much more considerable increase 
in directivity is obtained with larger reflection coeffi- 
cients. Curve 2 of Fig. 2 was calculated for a pp) =0.837 
(z.e., 70 per cent energy reflection in the normal direc- 
tion) and a decrease of capacitive field impedance pro- 
portional to cos a. Since fy) =213.2°, the resonance dis- 
tance 1, becomes 17.47 mm. Curves 3 and 4 were corre- 
spondingly calculated for po equal to 0.9 and 0.95, re- 
spectively. 

Tests confirmed these theoretical results. A rectangu- 
lar waveguide aperture (1X4 inch) in a plane conduct- 
ing screen (32X32 cm) was used as an antenna. The 
measured radiation pattern in the H plane is shown in 
curve 1 of Fig. 3. When a movable frame with a wire 


2G. G. MacFarlane, “Surface impedance of an infinite parallel- 
wire grid at oblique angles of incidence,” JIEE, Part IIIA, vol. 93, 
p. 1523; March—May, 1946. 
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Fig. 2—Calculated radiation pattern of a sheet with a reflection co- 
efficient of 1) po>=0.707, 2) po=0.837, 3) pPo=0.9, 4) po =0.95. 
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Fig. 3—Measured radiation pattern of a waveguide aperture with 
reflecting screen (1) and awire grid (2) orstrip-grating (3) disposed 
in front of it. 


grid similar to the one calculated above was placed in 
front of the aperture, the resonance distance /, was 
found to be 13.5 mm and a considerable increase in di- 
rectivity was noted. The measured pattern in the H 
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plane is shown by curve 2 of Fig. 3. This curve follows to 
a great extent the calculated curve (curve 1 of Fig. 2) 
but the directivity is somewhat greater. 

If a homogeneous dielectric plateglass (32 X32 X0.3 
cm) is used in place of the wire grid, a resonance distance 
of 15.9 mm results and the radiation pattern in the 7 
plane corresponds almost exactly to that of a grid. 
From the resonance distance we obtain ~o>~180°, which 
corresponds to a plate thickness of one-fourth of the 
wavelength in the dielectric. Therefore, the relative di- 
electric constant would be approximately 7 and 
por0.75. 

Patterns in the £& plane are considerably different. In 
that place, the array with the dielectric plate has a 3 db 
beamwidth of 35°, and a 10 db width of 73°, and at an 
angle of +72° there are side lobes of approximately 
—8 db. Apart from the considerably wider and more 
irregular individual pattern of the waveguide aperture 
with screen, the reflection factor of this sheet decreases 
with the increase in angle and passes through a mini- 
mum (Brewster angle). 

The radiation pattern of the array with the wire grid 
is strongly lobed in the & plane. Side lobes are as high as 
80 per cent at an angle of +13°, 40 per cent at +26°, 
and so on. 

As the angle of incidence increases, the wires of the 
grid are no longer parallel to the electric vector and re- 
flection decreases. Furthermore, in the finite grid, re- 
flections of induced currents are noticeable at the ends 
of the wires. In order to calculate these, a radiation 
originating from the edges may be added in a manner 
similar to the case of a slot-excited ground-plane an- 
tenna.*’ Due to the great width of the screen (10 \), a- 
rapidly varying field is superposed which, when multi- 
plied by (3), results in an irregular radiation pattern 
with pronounced ripples. | 

It is possible to obtain a considerably stronger reflec- 
tion, é.g., by means of a capacitive strip-grating placed 
vertically to the electric vector. The test arrangement 
consists of aluminum foil 12 mm in width, glued to the 
dielectric plate at intervals of 1 mm. Such a strip-grating 
in air gives a reflection factor of pp =0.89 for a plane 
wave.’ Due to the dielectric plate, the capacity and, 
consequently, the reflection factor are somewhat in- 
creased. With a resonance distance of J.=16.3 mm 
curve 3 (Fig. 3) results. Directivity in the E plane is 
somewhat less and the main beam is slightly irregular. 
The absolute gain amounts to 19.6 db. The adjustment 
of the resonance distance is very critical, and the sheet | 
as well as the screen must be perfectly plane and situ- 
ated parallel to each other. 

The influence of the edges on the main beam de- 
creases when both sheet and screen are smaller. A re- 

8 A. Dorne and D. Lazarus, “Very High Frequency Techniques,” 
Radio Res. Lab. Staff, McGraw-Hill Book Co. Inc., New York, N. Y., 
ch. 7-3 and 4, 1947. 


4N. Marcuvitz, “Waveguide Handbook,” Rad. Lab. Ser., Mc- 
Graw-Hill Book Co., Inc., New York, N. Y., vol. 10, ch. 5-18; 1951. 


Fig. 4—Arrangement of the perforated grid in front of the 
waveguide aperture and screen. 


flecting sheet consisting of a conducting plate with 
equal holes uniformly spaced (13 mm diam), as shown 
in Fig. 4, gives a power gain of 18.5 db and a half-width 
of approximately 18° in both planes. Fig. 5 indicates the 
radiation pattern in polar coordinates. Sheet and screen 
measure 23X24 cm and have a resonance distance of 
[,=14 mm. The measured radiation pattern (main lobe) 
is similar to the calculated curve 3 in Fig. 2, and also cor- 
responds approximately to the calculated reflection fac- 
tor of such a perforated grid.°® 

As regards the second resonance distance /,=30 mm, 
directivity still increases but the side lobes become 
greater, as may be determined from (3). 


FoRMS OF CONSTRUCTION 


In practice all forms of partially-reflecting sheets may 
be used. The choice depends to a great extent on the 
wavelength and use of the array. For decimeter and 
meter waves, reflecting sheets of perforated or wire 
grids with low wind resistance, may be used, whereas for 
centimeter waves closed sheets are preferable. These 
sheets can be made of low-loss insulating material and a 
sprayed, homogeneous, semipermeable metal coating, or 
with glued strips or disks of metal foil. 

The reflection coefficient should lie between about 
po=0.7 and 0.9. Smaller values show little effect and 
larger values often produce the disadvantages of high 
loss, excessive tolerance, and great dependence on fre- 


5 J. C. Simon, “Etudes de la diffraction des ecrans plans et applica- 
tion aux lentilles hertziennes,” Ann. de Radioelect., vol. 6, p. 205; 
July, 1951. 
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Fig. 5—Radiation pattern of the array of Fig. 4 in 
the H and E£ plane. | 


quency. The reflection coefficient should increase for 
angles of incidence that deviate from the normal direc- 
tion. 

There are two possibilities of further improving array 
characteristics. Primarily, the amplitude distribution of 
the radiating sheet may be influenced by different re- 
flection factors. Illumination is improved when the sheet 
produces strong reflection in its center and weaker re- 
flection toward the margin. The choice of an appropriate 
taper depends on the nature of the sheet, its over-all 
size, and the excitation. The composition of inhomo- 
geneous sheets is also important and diffraction can be 
used to improve directivity, for instance, on the capaci- 
tive strip-grating.® 

Secondly, lateral radiation, which occurs principally 
in smaller systems, may be reduced by completely or 
partly covering the sides with metal. This produces a 
cavity which is excited by the antenna in its interior and 
which is coupled to free space by the partially reflecting 
front sheet. The number of modes excited in the en- 
closed space, and their configuration, depends on the 
nature and location of the exciting antenna (or anten- 
nas), and on the dimensions of the cavity. 

Fig. 6 (a) is a sketch of a rectangular cavity excited 
in the principal TEi9 mode by a waveguide aperture in 
the back screen. Fig. 6 (b) is a sectional view of a simpli- 
fied distribution of the electric field corresponding to the 
second resonance distance of the sheet; the standing- 
wave distribution formed between back screen and 
front sheet is shown with a zero plane in the middle. 
Due to the presence of the metal side screens, the wave- 


6 G. von Trentini, “Btindelung elektrischer wellen durch leitschei- 
ben,” Z. angew. Phys., vol. 6, p. 462; October, 1954. 
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(a) 


length within the cavity is somewhat larger than with- 
out walls, and the resonance distance /, increases corre- 
spondingly. However, this influence is noticeable only 
with a small cavity in which the coupling of the wave- 
guide also modifies the resonance distance to a slight 
extent. 

As long as higher modes remain suppressed, the di- 
rectivity of such a “reflex-cavity antenna” will increase 
with increasing distances of the side walls and with a 
consequent enlargement of the sheet and radiating 
aperture. Depending on the nature of excitation, other 
modes may however be excited along the dimension B. 
The TE) mode distribution is to be preferred over the 
TE3 9 or TMy: mode. Fig. 6 (c) shows the approximate 
distribution of such a mode, the resonance of which has 
to be prevented for this affiliation. The resonance may 
be expected to occur at a side wall distance of 


r 


oN 2 
Vo) 

2A 
Depending on the field distribution, the exciting center 
is displaced a certain distance s (in this case about 
1—2 mm) from the center of the waveguide. It should be 
noted also that one of the reflecting walls for a mode 
along B is the partially reflecting front screen, so 
that (9) is only approximately valid; fortunately, the 
reflection coefficient may be expected to lie close to —1 
near grazing incidence. 

In addition to affecting the resonance distance, B also 
affects the impedance at the waveguide aperture and 
thus the matching of the cavity. Consequently, radia- 
tion of the array is subject to a periodically-varying in- 
fluence as B is increased and a favorable length must be 
found empirically. 
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Fig. 6—Reflex-cavity antenna with approximated E-field distributions. 


Test MODELS AND RESULTS 


In the test arrangements described in the beginning, 
sizes are large in comparison with wavelength and the 
gain is relatively small compared to geometric area. 
Smaller models that have higher efficiency are more im- 
portant for practical use. 

Fig. 7 (next page) presents three methods of construc- 
tion which give an effective area of over fifty per cent. 

The reflex antenna shown in the middle of Fig. 7 con- 
sists of a waveguide aperture (22.9 10.2 mm) in a re- 
flecting screen (63X63 mm) and a grid of 6 parallel 
wires (58 mm long; 20 =1.5 mm; d=12 mm). The reso- 
nance distance is /,=14.5 mm and two of the sides are 
closed by two adjacent wires (A =74 mm; 29 =1.5 mm; 
d=6 mm). The gain is approximately 14.3 db and 
the beamwidth in the H plane is 31.6° (3 db) and 58.5° 
(10 db), and in the £ plane 37° and 63°, respectively. 
The side lobes in the two planes are —14 and —13 db. 
This system is also convenient for decimeter waves, ex- 
citation being produced by a dipole and the reflecting 
screen being replaced by slightly separated parallel bars, 
or by a fine mesh wire screen. 

In the reflex antenna shown in the lower part of 
Fig. 7 reflection of the capacitive sheet is tapered. The 
antenna consists of a waveguide aperture, a reflecting 
screen (96X82 mm), and a polystyrene plate 1.8 mm 
thick with thin brass disks, 12 mm in diameter, glued 
to it. Their separation is around 15 mm in the center 
and 18 mm at the edge. Two side walls serve to hold the 
plate with the disks at a resonance distance of l,=18.5 
mm. This reflex antenna has a gain of approximately 
17.2 db. In the H plane the 3-db width is 20.5°, and the 
10-db width is 39.9°. In the E plane the 3-db width is 
22°, and the 10-db width, 42°. The side lobes amount to 
less than —14 db. 
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Fig. 7—Reflex antennas for \=3.2 cm. 


The upper part of Fig. 7 shows a closed reflex-cavity 
antenna with tapered reflection. It consists of two tele- 
scoped brass tubes. The tube having the larger diameter 
(55/60 mm) is closed by a metal screen in the middle of 
which the waveguide ends. The tube having the lesser 
diameter (50/55 mm) is closed in front by four brass 
strips, 8 mm wide, placed perpendicularly to the electric 
vector; these brass strips are 0.8 mm thick and are 
spaced apart 1 mm and 2.8 mm, respectively. The 
(adjustable) resonance distance is /,=18 mm. This 
antenna has a gain of about 14 db and its radiation 
pattern is shown in polar coordinates in Fig. 8. 


Von Trentini: Partially Reflecting Sheet Arrays 


671 


aH 


na 


Jf 


NO 


Fig. 8—Measured radiation pattern of the reflex-cavity 
antenna in the H and £ plane. 


To improve the match of the antenna to the wave- 
guide, all models can be provided with an inductive 
diaphragm in the opening plane of the exciting wave- 
guide. In the antenna just described, a diaphragm 
14.8 mm wide gives a vswr of 1.05:1, and the antenna 
efficiency is comparable to that of a horn-type antenna. 
However, with an approximately equal aperture, the 
length of the reflex antenna is appreciably shorter than 
that of the horn antenna. 

Of course, a group of these arrays may be fed from a 
common transmission line, or several distributed an- 
tennas may excite a single reflecting sheet. If radiation 
is to be effected in different directions, various partially- 
reflecting sheets may surround the antenna. Finally, the 
possibility exists of exciting other modes, thus producing 
additional radiation and polarization properties of the 
reflex-cavity antenna. 
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staff of the Hughes Aircraft Company, 
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At Hughes, Dr. Elliott specialized in surface 
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University of Texas 
in 1944 and 1947, re- 
spectively. From 
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Navy. 
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came a research as- 
sociate in the Elec- 
trical Engineering 
School of Cornell 
University. From 1951 until July, 1956, he 
was with the Electrical Engineering Research 
Laboratory, University of Texas, as a re- 
search engineer. At both Cornell University 
and the University of Texas, he was engaged 
primarily in tropospheric propagation 
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In 1956, he received the Ph.D. degree 
from the University of Texas, and became a 
member of the teaching and research staff of 
the electrical engineering department of the 
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for ballistic measurements and the evalua- 
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serving as staff assistant to the chief engi- 
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Engineering and Ap- 
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problems under the direction of Prof. 
Ronold King, and has published several 
reports. 
Miss Kennedy is a member of the Fed- 
eration of American Scientists. 
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He received the B.S. degree in 1945 and the 
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from the University 
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ing assistant at the 
University of Cali- 
fornia from 1947 to 
1949, as well as a re- 
search engineer for 
the U. S. Navy An- 
tenna Project at the 
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Company, where he was engaged in the de- 
velopment of linear and two-dimensional 
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and other microwave antennas and com- 
ponents. Mr. Kurtz is presently head of the 
microwave Antenna Department of Rantec 
Corporation. 

He is an associate member of Sigma Xi. 
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Arsenal and the Naval Research Labora- 
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the Naval Air Development Center, engag- 
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missile systems, in particular, with contribu- 
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tions in the electromagnetic theory of ra- 
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to the Air Warfare Research Department 
of the same Center, where he is now par- 
ticipating in operational research. 
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Peru, Ill., in 1921. He received the B.S. de- 
gree in electrical engineering from the Uni- 
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development, and production of microwave 
antennas and components, chiefly in the 
automatic tracking radar field. 
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ety. 


R. J. STEGEN 


COL 


October 


Giswalt von Trentini was born on April 
21, 1913, in Munich, Germany. He gradu- 
ated as Engineer-Physicist from the Univer- 
sity and Institute of 
Technology in Mu- 
nich in 1939. From 
1940-1945 he was 
with the “Drahtlose 
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was on the staff of the “Centre National 
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ciones Scientificas Y Tecnicas de Las 
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There he works on radar antennas, wave- 
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Charles P. Wells was born January 23, 
1909, in Tabor, Iowa. He received the B.A. 
degree with a major in mathematics from 
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and Spin Wave Equations: AP261 
Phase Centers of Microwave Antennas: 
AP323 
Phase Shifters, Scanning Antenna Errors: 
AP327 
Plasma Oscillations AP262 
Prolate Spheroid: AP192, AP230 
Antenna, Radiation from: AP192 
Radar Cross Section: AP230 


R 


Radar Cross Section of a Prolate Spheroid: 
AP230 
Radiation: AP192, AP207, AP225, AP245, 
AP246, AP254, AP330 
Cerenkov and Undulator: AP245 
by Disks and Conical Structures: AP330 
Electromagnetic, Patterns and Sources: 
AP225 


Nonreflecting Absorbers: AP246 
from Prolate Spheroidal Antennas: AP192, 
AP230 
from Ring Quasi-Arrays: AP254 
of Spherical Luneberg Lens: AP207 
Reflections: AP195, AP212 
F-Layer, and Long Range Meteoric 
Echos: AP195 
Microwave, from the Ocean: AP212 
Refraction, Superrefractive Layers, Effect 
on Nonoptical Fields: AP214 
Resonators, Electromagnetic, Variational 
Principles: AP201 


iS) 
Scanning: AP327, AP328 
Conical, X-Band Antenna; AP328 
Errors of Antenna Arrays: AP327 
Scattering: AP194, AP209, AP220, AP227, 

AP229, AP230, AP235, AP236, 
AP239, AP243, AP244, AP264, 
AP265, AP266 

Approximate Method for 
A227 

Backward, Schwinger Variational Prin- 
ciple: AP236 

Boundary Value Problems: AP264—-265 

Diffraction Theory for Cylinders and 
Spheres: AP220 

Forward: AP266 

Forward and Back, from Certain Rough 
Surfaces: AP244 

by an Infinite Grating: AP239 

Ionospheric, Dependence on Angle of 
Scatter: AP243 

Multiple: AP209, AP266, 

Panel Discussion on: AP266 


Problems: 


Chapter News 


Chicago, IIl.: October, p. 588 : 

Columbus, Ohio, Chapter (to be organized): 
April, p. 101 

Denver-Boulder, Colo.: October, p. 588 

Los Angeles, Calif.: October, p. 588 

Philadelphia, Pa.: October, p. 588 

San Diego, California: April, p. 101, Octo- 
ber, p. 588 

Washington, D. C.: April, p. 102, October, 
p. 588 


Meetings 


Calendar of Events: January, p. 4; April, 
p. 102 


by Randomly Distributed Obstacles: 


AP209 

by Prolate Spheroid, Radar Cross Sec- 
tion: AP230 

Scattered Fields and Scattered Power: 
AP194 


Theory, Vector Combining Electric and 
Magnetic Fields: AP229 
Variational Method in Problems: AP235 
Screens, Absorbing Plane, Diffraction of 
Electromagnetic Waves Caused by 
Apertures: AP203 
Seismic Pulse: AP260 
Sheet Arrays, Partly Reflecting: AP335 
Signals, Microwave Reflection from the 
Ocean: AP212 
Spin Wave Equations and Permeability 
Tensors: AP261 
Superrefractive Layers, Effect on Nonopti- 
cal Fields: AP214 
Surface Currents: AP189, AP197 
Controlled by Use of Channels: AP197 
Excited by Half-Planes: AP189 
Surface Wave: AP233, AP250 
Antennas, Spherical: AP250 
Occurrence of : AP233 


At 


Telescopes, Radio, Resolution, Pattern Ef- 
fects and Range: AP253 
Transformers, Simulated, Lens Matching 

by: AP187 
Transmission Characteristics of Inclined 
Wire Gratings: AP332 
Tropospheric Propagation, New York Uni- 
versity Research: AP228 


Nontechnical Index 


Communication by Scatter 
Symposium: January, p. 2 

Communications Theory and Antenna De- 
sign Symposium: October, p. 591 

International Radio Consultative Commit- 
tee: October, p. 589 

IRE National Convention, 1956, January, 


Techniques 


p. 4 

IRE National Convention, 1957: October, 
p. 591 

Microwave Techniques Symposium: Janu- 
ary, p. 3 

Optics and Microwave Symposium, October, 
p. 588 

URSI: October, p. 588 

Very Low Frequency Propagation Symposi- 
um: October, p. 589 


Ce) 


AW 


Waveguides: AP201, AP211, AP228, AP257, 
AP263 
Amplitude Concept Applied to Junction 
Problems: AP211 
Circular, with Gyromagnetic Material, 
Propagation in: AP257 
Electromagnetic, Variational Principles: 
AP201 
New York University Research: AP228 
Rectangular, Ferrites in: AP263 
Wave Propagation: AP196, AP214, AP228, 
AP238, AP240, AP241, AP242, 
AP257, AP260, AP268, AP334 
in Circular Waveguides Filled with Gyro- 
magnetic Material: AP257 
in Doubly-Refracting Media: AP268 
Effect of Superrefractive Layers on 
Nonoptical Fields: AP214 
Height-Gain Curves at 8.6 MM.: AP240 
Ionospheric: AP196, AP228 
New York University Research: AP228 
over Irregular Terrain: AP196 
Line-of-Sight: AP238, AP334 
Near-Field Corrections to: AP238 
in Randomly Inhomogeneous Medium: 
AP334 
Particle-Energy Relations: AP242 
Phase Measurements of Varying Signals 
over Turbulent Paths: AP241 
Seismic Pulse: AP260 
Tropospheric, New York University Re- 
search: AP228 
WKB Method, Refinement of and Ap- 
plication to Electromagnetic Wave 
Theory: AP226 
Yagi-Uda Antenna: AP252 


WESCON Announced: April, p. 102 


Miscellaneous 


Administrative Committee Responsibilities 
and Activities: January, p. 1 

Communications Section of PGAP Trans- 
actions: January, p. 4; April, p. 102 

Fringe Interest Groups in PGAP: January, 
paz 

Membership Survey: January, p. 2, April, 
p. 101 

News and Views Department Needs Sub- 
Editors: January, p. 4 

Transactions Policy: January, p. 2. 
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